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1 .  INTRODUCTION  AND  REVIEW  OF  PREVIOUS  LITERATURE 


The  purpose  of  this  report  Is  to  provide  students  with  a  means  to  calcu¬ 
late  the  energy  levels  of  Impurity  ions  in  solids  and  analyze  experimental 
results.  In  order  to  achieve  this  goal  it  is  necessary  to  review  material 
which  many  readers  may  have  been  exposed  to  elsewhere.  This  review  is  rather 
brief,  and  the  bibliography  (sect.  1.4)  in'cludes  references  to  specific  sec¬ 
tions  of  textbooks,  monographs,  or  papers  where  the  material  is  discussed  in 
more  detail.  Unfortunately,  as  in  most  fields  of  science,  a  number  of  dif¬ 
ferent  notations  (in  various  alphabets)  are  used  in  the  literature;  it  is 
hoped  that  this  review  will  help  overcome  some  of  this  difficulty. 

Many  mathematical  results  are  presented  without  proof  with  a  view  towards 
brevity.  The  review  of  angular  momentum  is  directed  towards  the  use  of  the 
various  tabulations  of  matrix  elements,  n-j  symbols,  and  group  tables,  and  our 
study  of  group  theory  is  simply  directed  toward  the  use  of  the  abundant  tables 
of  the  32  point  groups. 

Most  of  the  discussion  of  crystal-field  theory  is  devoted  to  the  point- 
charge  or  point-multipole  model.  This  is  quite  natural  since  most  of  its 
recent  development  has  been  done  by  research  workers  at  Harry  Diamond  Labora¬ 
tories  (HDL) .  Very  recently,  research  workers  under  the  direction  of  Paul 
Caro  in  France  have  significantly  extended  the  point-multipole  crystal-field 
model  and  their  contributions  are  also  discussed. 

Finally,  the  ongoing  work  at  HDL  on  the  transition-metal  ions  is  dis¬ 
cussed.  Most  of  these  results  are  unpublished  at  present,  but  computer  pro¬ 
grams  are  used  to  illustrate  directly  the  effect  of  various  free-ion  and 
crystal-field  interactions. 

1 . 1  The  Hydrogen  Atom 


We  assume  that  we  have  an  electron  of  charge  -e,  mass  m,  and  an 
infinitely  massive  nucleus  of  charge  Ze.  The  nonrelativistic  Hamiltonian  for 
this  system  is 


2 

EL  . 
2m 


(1.1) 


where  p  is  the  momentum  of  the  electron.  The  time-independent  Schrodinger 
equation  for  this  system  is  given  by 

Hi|>  ■  , 
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where  we  have  used  the  relation  ?  -  -ihV.  The  bound  solutions  to  equation 
(1.2)  satisfying  all  the  boundary  conditions  are 


^nim 


Vr)Vs"*) 


(1 .3) 


where  the  R  »(r)  are  the  associated  Laguerre  polynomials  and  the  (0»<|O  are 
spherical  harmonics.  The  radial  functions  R^tr)  are  of  little  interest  here 
and  are  not  discussed  further.  The  energy  in  equation  (1.2)  is 


E  .  -Z2me4 
n  2n2*2 


(1.4) 


where  n  £  1.  The  quantum  number  l  (angular  momentum)  is  restricted  to  the 
values 


&<n,  l  «  0,  1 ,  2  . . .  ,  (1.5) 

and  m,  frequently  referred  to  as  the  magnetic  quantum  number,  is  restricted  to 

-nun  .  (1.6) 

The  spectroscopic  notation  for  a  sequence  of  l  values  is 

i  -  o,  1,  2,  3,  4,  5,  6,  7  ...  (1.7) 

a,  p,  d,  f,  g,  h,  i,  k  ...  ; 

along  with  the  value  of  n  the  states  are  referred  to  as  Is,  2s,  2p,  3a  ... 

with  the  restriction  given  in  equation  (1.5). 

The  spherical  harmonics  are  given  by 


-  Nln,V003  9)e1”* 


where 


The  ( cos  9)  are  associated  Legendre  polynomials  and  are  defined  as 

vz>  -  (i  -  z2)"/2  (y«  v*>  • 

and  the  Legendre  polynomials  are 


(1.8) 


(1.9) 


(1.10) 
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The  definitions  given  in  (1.9)  and  (1.10)  are  restricted  to  m  i  0.  For  m  <  0 
we  have 


Y 


i-m 


(1.11) 


where  Y£m  is  the  complex  conjugate  of  Y^. 

The  spherical  harmonics  are  normalized  so  that 


Y# 

Vi 


,(9,*)Y4m(0,*)  d a 


VfcVm 


(1.12) 


where  dfl  =  sin  0  d9  d<j>,  and  the  integration  covers  the  range  OSeSir,  0  £  <|> 
S  2 it. 

1 . 2  Angular  Momentum  Algebra 

In  classical  mechanics,  the  angular  momentum  of  a  particle  is  defined 
by 

I  *  r  *  p  .  (1.13) 

Actually,  we  should  specify  that  the  angular  momentum  so  defined  is  about  a 
particular  origin,  and  r  i3  the  vector  distance  from  this  origin  to  the  par¬ 
ticle  with  momentum  p. 

If  we  use  the  commutation  relations 


[xifPj]  -  itfS^ 


(1.14) 


with  x^  »  x,  y,  or  z,  then  we  can  obtain  the  commutation  rules  for  angular 
momentum, 


[l  A  ]  -  Ifil 
L  x  y  z 


[£,«,]-  IfiH  ,  arid  [«,  ,i  ]  -  i-TU 
*-y  zJ  X  UZ  XJ  y 


(1.15) 


which  are  the  basic  commutation  rules  for  the  Cartesian  components  of  the 
angular  momentum.  For  convenience  here  we  shall  drop  the  1i  in  the  commutation 
relations.  This  does  not  mean  that  we  drop  ft  throughout;  we  restore  Ti  simply 
by  writing  the  interactions  involving  the  angular  momentum  so  that  the  -ti  is 
contained  in  tne  constants.  As  an  example  of  this,  consider  the  spin-orbit 
Hamiltonian 


(1.16) 


with  t  and  8  having  units  of  angular  momentum  (the  spin  angular  momentum, 
a,  we  v/i  1 1  diaouss  later).  When  these  are  written  in  terms  of  dimensionless 
t  and  s,  we  have 


w  i  au  *  5 

"2  *  .  2  2  r  3r  *’s  ’ 

2m  c 


(1.17) 


where  ?  and  s  obey  the  commutation  rules  in  equation  (1.15)  but  -  1 . 


For  our  purposes  here,  we  want  to  use  the  spherical  representation 
of  t,  which  is  given  by 

‘♦1  ■  ■  •=  (1X  *  Uy>  • 


*0  •  *z  ’ 

*1  -  4  u  -  iO  . 


(1.18) 


and  the  commutation  rules  are 


l^o’^+I  ^  “  *+l  ’ 

’  *-l  » 

[*-+1  »Li ]  “  _)lo  • 


(1.19) 


The  eigenfunctions  of  the  angular  momentum  are  the  spherical  har¬ 


monics,  and 


where 


Aq|  im>  -  m|im>  , 


i,|i.m>  -  ?  —  [(£?m)(J(,±m+1  i]1 /2|)l,m±1>  , 


|tm>  -  Y.m(0,«)  . 


(1.20) 


Frequently,  we  shall  use  the  unit  vector  P  to  indicate  the  argument  of  Y^, 
thus: 


When  the  Yn  are  wave  funotiona  suoh  aa  in  equation  (1.20),  we  have 
sun 


Yim(r)  -  \lm>  . 


la 


The  orthogonality  of  the  wave  funotiona  aa  given  by  equation  (1.12) 


<9,  'm'  |  i,m> 


{U,6mm'  * 


Further,  we  3hall  aaaume  that  the  spin  angular  momentum,  a,  obeya  the  aame 
oommutatlon  relationa  aa  given  in  equation  (1.19);  the  two-component  apinor 
wave  funotiona  are  repreaented  by  the  wave  function  |sm  >,  ao  that  the  aingle- 
electron  wave  function  for  orbital  and  spin  angular  momentum  ia 

|£m0>|am  >  .  (1.21) 

M  3 

The  wave  funotiona  given  by  equation  (1.21)  then  obey  the  following: 


(t)2|2m^>|ama>  -  i(4+1 )|£ml>|sm>  , 

a0|im!l>|sma>  -  mg|  lmft>  j  smg>  , 


(1.22) 


(a)  |i,ml>jama>  -  s(a+'i )  |  im^  |sma>  , 

where,  of  course,  s  -  1/2.  A  further  property  of  the  spherical  harmonics  is 
given  by 


IY£B<r)  -  (-D^m(r)  ,  (1-23) 


where  the  Inversion  operator  is  IP  -  -f,  a  property  that  will  be  used  fre¬ 
quently  in  our  analysis.  For  other  symmetry  operations,  an  explicit  expres¬ 
sion  for  the  spherical  harmonics  is  convenient;  table  1.1  is  included  for  this 
purpose.  While  many  of  the  Interaction  terms  of  the  Hamiltonian  were  derived 
by  using  spherical  harmonics,  it  is  convenient  to  introduce  the  tensor 
operators 


<Vr) 


\2Jl+1  / 


1/2 


Vr) 


(1.2l») 
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TABLE  1.1.  SPHERICAL  TEN80R3 ,  Cnm,  IN  RECTANGULAR  COOPDIHATF.S 


n,a 

M4 

^nm 

n,m 

M4 

^nm 

0  0 

1 

1 

5  5 

-h« 

(x  ♦  ly)1’ 

1  1 

-  J. 

’  ✓? 

x  ♦  ly 

5  H 

hr™ 

z(x  ♦  ly)4 

1  0 

1 

z 

5  3 

(9z2  -  r2)(x  ♦  iy)3 

5  2 

✓210 

5 

z(3*2  -  r2)(x  ♦  iy)2 

2  2 

✓37? 

(x  ♦  iy)2 

5  1 

(21z4  -  1i)z2r2  ♦  r4)(x  ♦  ly) 

2  1 

-✓372 

z(x  ♦  iy) 

5  0 

1 

5 

63z4  -  70z2r2  ♦  15r4 

2  0 

1 

2 

3z2  -  r2 

6  6 

✓23T 

32 

(x  ♦  ly)6 

3  3 

(x  ♦  ly)3 

6  5 

-*g 

z(x  +  iy)5 

3  2 

/tv? 

z ( x  ♦  ly)2 

G  it 

( 1 1 z2  -  r2)(x  +  iy)4 

3  1 

-I 

(5z2  -  r2) (x  +  iy) 

6  3 

✓T05 
"  "TS" 

zfllz2  -  3r2)(x  ♦  ly)3 

3  0 

1 

2 

5z2  -  3r2 

6  2 

✓T05 

32 

(33z4  -  I8z2r2  ♦  r4)(x  ♦  iyi2 

14  It 

✓70 

“TS 

(x  ♦  iy)4 

6  1 

iKz 
’  T? 

z(33z4  -  30z2r2  ♦  5r4)(x  +  ly) 

«  3 

z(x  +  ly)3 

6  0 

1 

r? 

231z6  -  315z4  +  105z2  -  5 

A  2 

✓To 

(7z2  -  r2)(x  ♦  iy)2 

H  1 

✓5 

■  nr 

z(7z2  ••  3r2)(x  +  iy) 

14  0 

i 

? 

35z4  -  3°z2r2  +  3r4 

a Multiplier  to  entry  on  the  right. 


Since  Y*  (r)  -  (-1 )mY0  (r),  we  have 
Urn  Jt.-m 

C •  (?)  -  (-1)%  (?)  .  0.25) 

dm  x.,-m 

m 


xn  ic^LiwwwvArjsf .\v.vais  vwx  wv.  u>  v«  i-h  itk  w>  l>  ,-xj.jum m  ■'*r*r*r*j  awa  a 


The  use  of  Cj,m  rather  than  Yj^  in  the  interaction  terms  eliminates  almost  all 
the  faotors  of  An  example  of  this  is  the  coupling  rule  for  spherioal 
harmonics  (Rose,  1957,  p  61):# 


1/2 

Vn»  ■  l  [tjHfTT1]  <*«»"«»  I ‘<°»  <k(q)n(.)|«q*.)>  , 

(1.26) 


but 

CkqCnm  "  ^  <><(0)n(0)  |  4(  0)>  <k(q)n(m)  |  i(q*m)>  Cfc  .  (1.27) 

fir 

In  equations  (1.26)'  and  (1.27),  all  the  tensor  operators  have  the  same  argu¬ 
ment.  The  quantities  in  angular  brackets  in  equations  (1.26)  and  (1.27)  are 
Clebsoh-Gordan  (C-G)  coefficients,  which  we  oover  in  seotion  2. 

1 .3  Problems 

1.  The  inversion  operator  I  converts  the  vector  r  to  -r.  Show  that 


IYom(r)  -  (-1)*Yim(r)  . 


2.  If 

C2(x) (x, 
C2(y) (x , 

y*’ 

C2(z)(x, 

C2(1)(x, 

C2(2)(x, 

o(]_x)(x, 


y,  z)  -*•  (x,  -y,  -z)  , 

y,  z)  -*•  (-x,  y,  -z)  , 

y,  z)  +  (x,  v.  -z)  , 

y,  z)  •>  (-x,  -y,  z)  , 

y,  z)  ■*  (y,  x,  -z)  , 

y,  z)  -*•  (-y,  -x,  -z)  , 

y,  z)  -*•  (-X ,  y,  z)  , 


o(]_y)(x,  y,  z)  -►  (x,  -y,  z) 


•References  are  listed,  alphabetically  by  author,  at  the  end  of  each  section,  along  with 
uncited  bibliographic  entries  pertinent  to  each  topic. 


o(Jj)(x,  y ,  z)  +  (-y,  -x,  z)  , 


0(^2) (x,  y,  z)  +  (y,  x,  z)  , 


show  that 


C2<*)Yim  -  <-1)*Y*-m(r)  , 


C2(y)Y4m(r)  -  (-1)^Ya>m(r)  , 


°h**mtf>  "  <-1>*+mYim(r)  , 


C2(‘>Ck,  ■  •1'\,  ’ 


C2(,)Ckq  '  (-,)kel("/2>\-q  • 


C2(2)Ckq  ■  <-'>"o‘('/2)Vq  - 


"(lx,Ckq  -  IC2(k)Ckq  -  Ck.Q  , 


'<ly>cka  '  lcp<y)^„  -  <-’)qS-q  • 


»(ll)C.  -  e 


kq  2  w  kq 
l(Tr/2)q, 


C. 


kq  *  v'k-q  ' 


-i  (ir/2)qr 

kq  w  Ck-q  ’ 


0<l2)C„n  -  e  . ~''C 


by  using  equations  (1.8),  (1.9),  and  (1.10),  or  by  using  table  1.1. 


3.  Using  the  generating  function  for  P^(cos  6),  obtain  the  result 


rrr2 


-  I 


r< 


r> 


1+1  l 


p.(cos  e)  , 


where  cos  0  -  ?i*?2  and  r<  is  the  smallest  of  the  vectors  r^  and  r 2- 
addition  theory  of  spherical  harmonics,  we  have 


From  the 


P,(cos  6)  -  I  c;B(i;,)cJ?2)  . 


m--i, 
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j.-vt  miiru  aru  srv.  *ru  w  1 r  j  *n*w\f*nj - 1— n  *rj  w  f  -j,  »r^ 


rj  wu  iru  wv. 


The  generating  funotion  for  the  Legendre  polynomials  is  given  in  Rainville 
(I960)  as 


o  TTo  "  1  n  P0(z)  *  with  |h|  <  1  and  z  -  oos  9 

2  Ml _ 1  '  f  "  MM  ™ 


[l  +  h  -  2hz] 


a-o 


4.  From  the  generating  funotion  for  P^(z),  show  that  P^(1)  -  1. 


P#(-1)  -  (-1)*'*  and 


i ,  ■ 1  • 

m— & 

it  4  4 

Show,  in  two  ways,  that  the  sum  in  problem  3  is  M)  when  rg  -  -r1 ,  for 
arbitrary  i. 

5.  By  expanding  r  *  t  show  that 


r,  .t  3  ir  *  i 
v  •  r  —  1 1  ■  — * — 
3r  r 


which  is  a  convenient  form  of  the  V  operator  for  spherical  problems. 

1 .4  Annotated  Bibliography  and  References 

Condon,  E.  U.,  and  H.  Odabasi  (I960),  Atomic  Structure,  Cambridge  University 
Press,  Cambridge,  U.K.  Chapters  3  and  4  give  a  thorough  discussion  of  the 
hydrogenic  wave  function.  Pages  190  and  191  tabulate  the  radial  wave 
functions. 

Judd,  B.  R.  (1963),  Operator  Techniques  in  Atomic  Spectroscopy,  McGraw-Hill, 
New  York,  NY.  This  outstanding  textbook  is  frequently  referred  to  in  this 
report;  the  book  is  almost  completely  free  of  typographical  errors.  The 
problems  in  general  are  very  difficult  and  require  considerable  time. 
Only  the  preface  and  the  first  two  pages  are  pertinent  here. 

Leighton,  R.  (1959),  Principles  of  Modern  Physics,  McGraw-Hill,  New  York,  NY, 
chapter  5,  The  One  Electron  Atom. 

Polo,  S.  R.  (1961,  June  1),  Studies  on  Crystal  Field  Theory,  Volume  I— Text, 
Volume  II— Tables,  RCA  Laboratories ,  under  contract  to  Electronics  Re¬ 
search  Directorate,  Air  Force  Cambridge  Research  Laboratories,  Office  of 
Aerospace  Research,  contract  No.  AF  1 9(604)— 5541 .  [Volume  II  gives  date 
as  June  1,  1961.]  The  symmetry  operations  are  discussed  on  pp  1-4ff; 
Clebsoh-Gordon  coefficients  on  pp  8. Iff;  and  excellent  tables  of  P^(z), 
Pim(z),  and  Yim  are  given  in  the  appendix,  all  in  Vol.  I. 

Rainville,  E.  D.  (I960),  Special  Functions,  Macmillan,  New  York,  NY.  This  is 
just  one  of  the  numerous  texts  written  by  this  outstanding  teacher. 
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Rose,  M.  E.  (1957),  Elementary  Theory  of  Angular  Momentum,  Wiley,  New  York, 
NY,  chapter  II. 

Sobelman,  I.  I.  (1979),  Atomic  Spectra  and  Radiative  Transitions,  Springer- 
Verlag,  New  York,  NY.  This  is  a  very  excellent  monograph  and  we  fre¬ 
quently  refer  to  Sobelman* s  derivations.  Beware  of  typographical 
errors  1!  Pages  1  through  12  are  pertinent. 

Watanabe,  H.  (1966),  Operator  Methods  in  Ligand  Field  Theory,  Prentice-Hall, 
Englewood  Cliffs,  NJ.  Pages  11  and  12  are  applicable  here;  also,  the 
Introduction  is  interesting.  We  frequently  refer  to  this  monograph.  Most 
of  the  equations  are  free  of  typographical  errors.  Tables  of  explicit 
expressions  for  Yj^m  for  0  S  i,  S  6  are  given  in  appendix  1.2. 


2.  CLEBSCH-GORDAN  COEFFICIENTS 


For  our  purpose,  it  is  convenient  to  define  the  Clebsoh-Gordan  (C-G) 
coefficients  as  the  coefficients  in  the  transformation  from  two  angular  momen¬ 
tum  spaces,  say,  I  and  s,  to  form  the  composite  space  j.  That  is, 


|  Jm>  -  l  <£(y)s(m-y)  |j  (m)>  |£.u>|s,m-y>  ,  (2.1) 

y 

where  the  quantity  <£(y)s(m-y) |j (m)>  is  a  C-G  coefficient.  The  limits  on  the 
sum  in  equation  (2.1)  are  not  given,  as  it  is  assumed  (and  will  be  assumed  in 
the  following)  that  the  sums  oover  all  values  for  which  the  C-G  coefficient 
does  not  vanish.  Since  we  wish  to  have  an  orthonormal  basis,  we  have 


since 


<j 'm'  |  Jm>  -  Sjj,V 

-  I  <£(u)s(m-y)  |  J  (m)>  <£(y')s(m,-y,)|j '(m')> 
yy’ 

<fcy'Uy>  -  6  ,  and  <s  ,mf  -y '  Is  ,m-y>  -  • 

1  yy'  1  m-y^'-p' 


(2.2) 


(2.3) 


Thus ,  we  have 


l  <£(y)s(m-y) |j (m)>  <£( y )s (m-y ) | J ' (m) > 
y 


(2.4) 


an  important  arid  very  useful  result.  If  we  ajsume  (correctly)  that  the  same 
coefficients  connect  the  j  space  to  the  i  and  s  spaces,  we  can  obtain  another 
condition  on  the  C-G  coefficients,  namely, 


’ll 


'®ss' 


I  <£(mjt)s(ms)  |j(raJl+ms)> 


<i'(m£)s'( 


<ng)  |j( 


(2.5) 
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Some  other  relations  among  C-G  ooeffioients  are 


<a(a)b(8)  |o(Y)>  -  0 

if  | a |  >  a,  or  |g|  >  b,  or  |y|  >  c,  (2.6) 

and  if  Y  «•  o  +  6  . 

The  C-G  coefficients  vanish  unless  the  three  angular  momenta  obey  the  triangle 
condition,  that  is,  |a-b|Sc£a  +  b  for  any  permutation  of  a,  b,  or  c. 
Three  of  the  most  important  symmetry  relations  of  C-G  coefficients  are 


<a(a)b(8)|o(Y)>  -  (-1)a+b“c  <a(-a)b(-8)|o(-Y)> 

-  (-1)a+b'°  <b( 8 )a(a) |c( Y) > 

1  /2 

-  <-1)a"“  (2g— -j)  <a(a)c(-Y)|b(-8)> 

where  a  +  8  -  Y  . 


(2.7) 


The  coupling  procedure  given  in  equation  (2.1)  can  be  deceptive.  If  we  had 
formed  the  wave  function  jjm>  by  coupling  s  to  i  as 


|jm>  -  l  <s(y)a(m-y)  Jj  (m)>  |sy>|il,m-y>  ,  (2.8) 

p 

the  relation  to  equation  (2.1)  would  be,  from  equation  (2.7), 


|jm>  -  l  (-1  )s+i  ^<i,(y)s(m-y)  |  Jm>  |fcy>|s,m-y>  , 

P 

in  which  we  have  changed  the  summing  index  y  m-y.  This  overall  phase  factor 
may  have  no  effect,  but  in  a  long  involved  problem  a  switch  from  i-s  to  s-2. 
coupling  can  cause  errors.  Algebraic  expressions  for  <a(a)b(8) | c(a+8) >  are 
given  by  Rose  (1  957,  pp  224-225),  for  c  -  1/2,  1;  algebraic  expressions  for 
the  related  3“j  symbols  for  1/2  £  c  £  2  are  given  in  Brink  and  Satchler  (1962, 
p  36).  A  few  special  C-G  coefficients  of  interest  are 


,  _/?  /2c  +  1\1/2  T(s) 

<a(0)b(0)  c(0)>  -  (-1)s/2  ( - )  — r-r - r-r — r 

\s  +  1/  T9,  Ts2Ts3 


(2.9) 


where 


-a+b+o,  Sg-a-b+c,  s^-a+b-c,  s-a+b+c  , 
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and 


T(s) 


<a(a)b(|3) |a+b(a+8)> 


9 


(bj-birrHrr  •  <2->°> 

<a(a)0(0)|b(8)>  -  6a  >6  0  , 

ct  |  D  Ot  f  p 

a(o)k+2(o)|Ko)>  -  -  {£i  [rlf^rffrH^k]1/2  <*(°)k<0>l*(0>>  • 

The  commutation  relations  for  the  spherical  components  of  the  angular 
momentum  of  a  single  electron  given  in  equation  (1.19)  can  be  written  com¬ 
pactly  in  terms  of  C-G  coefficients  as 

[*y,*v]  -  V2  <1(v)1(y)|l(n+v)>  and  (2.11) 

[sy,sv]  -  S2  <1(v)1(u)|1(m+v)>  sy+v  .  (2.12) 


The  total  angular  orbital  momentum  operator  for  a  system  of  N  electrons  is 


N 

2  - 1  i(  i )  , 

i 


(2.13) 


and  the  total  spin  angular  momentum  operator  is 


N 

§  -  I  s(i )  . 

i 


(2. 14) 


The  spherical  components  of  these  operators  obey  the  same  commutation  rela¬ 
tions  as  equations  (2.11)  and  (2.12),  or 


[L.Lv]  -  /2  <1(v)1(u)|1(m+v)>  Ly+v  and 


[Sy,Sj  -  /2  <1(v)1(p)|Ky+v)>  Sy+v 


(2.15) 

(2.16) 
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Also,  it  should  be  noted  that 


[L^,  S-v(  i  )  ]  -  /2  <1(v)1(y)  |Ky+v)>  £y+v(i)  and  (2.17) 

[S  ,av<i)]  -  /2  <1(\>)1(y)|l( y+v)>  3y+v(i)  .  (2.18) 

Consequently,  from  using  equations  (2.15)  through  (2.18),  we  have 

[L  ,C.(i)]  -  /k(k+1  )"  <k(q)1  (y)  |k(q+y)>  C.  .  (i)  ,  (2.19) 

y  Kq  k , qTy 

[Ju.ckqa)l  -  [ivVn]  (2-20) 

with 

3  -  L  +  $  . 

The  Clebsch-Gordon  coefficients  used  here  are  related  to  the  3-J  symbols  by 

<a(a)b(0)]c(Y)>  -  /2?H  (-1)“a+t>"Y  £  _°)  (2,21) 

and  the  symmetry  conditions  on  the  3~j  symbol  can  be  obtained  from  equation 
(2.7).  The  3-J  symbols  are  extensively  tabulated  by  Rotenberg  et  al  (1969). 

2. 1  Problems 

1 .  Show  that 

t2|im>  -  l(t+1 )|tm>  , 

2 

s  |sm>  -  s(s+1 ) | sm> 

2.  Consider  the  interaction  (spin-orbit  interaction) 


Show  that  the  matrix  elements  of  this  interaction  using  the  states  given  in 
equation  (2.1)  are 

<j'm'|Hg_o|jm>  -  |  Cj(J  +  1)  -  Ml+1)  -  +  • 

Hint:  J  -  t  +  s;  consider  J2. 
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2.2  Annotated  Bibliography  and  References 


Note:  Check  your  favorite  quantum  mechanics  text;  there  may  be  a 

section  on  C-G  coefficients  or  3-J  symbols. 

Brink,  D,  M.,  and  G.  R.  Satohler  (1962),  Angular  Momentum,  Clarendon  Press, 
Oxford,  U.K. 

Condon,  E.  U.,  and  H.  Odabasi  (1980),  Atomic  Structure,  Cambridge  University 
Press,  Cambridge,  U.K.,  appendices  2  and  3« 

Condon,  E.  U.,  and  G,  H.  Shortley  (1959),  The  Theory  of  Atomic  Spectra,  Cam¬ 
bridge  University  Press,  Cambridge,  U.K.  The  C-G  notation  used  by  Condon 
and  Shortley  is  related  to  that  used  in  this  report  as  follows: 

^ J 1 J  2m  1  I J  1 J  2 J m ^  ’  <J  1  (m1  )J2 Cm2)  1 J  <m)  >* 

Edmonds,  A.  R.  (1957),  Angular  Momentum  in  Quantum  Mechanics,  Princeton  Uni¬ 
versity  Press,  Princeton,  NJ.  Relationship  of  the  Y^m  of  various  authors 
is  given  on  page  21.  His  c(k)  are  the  same  as  those  of  Judd  and  are  the 
same  as  our  Ck  .  The  relation  of  the  C-G  coefficients  to  other  notations 
is  given  on  page  52.  This  is  a  good  book — but  look  out!  It's  loaded  with 
typographical  errors. 

Rose,  M.  E.  (1957),  Elementary  Theory  of  Angular  Momentum,  Wiley,  New  York, 
NY,  chapter  III.  The  relation  of  C-G  coefficients  to  other  symbols  is 
given  on  page  Hi.  The  commutation  rules  for  in  fcer’m3  of  C-G  are 

given  on  pages  8H  and  85. 

Rotenberg,  M.,  R.  Bevins,  N.  Metropolis,  and  J.  K.  Wooten,  Jr.  (1969),  The  3-J 
and  6-j  Symbols,  MIT  Press,  Cambridge,  MA. 
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3. 


WIGNER-ECKART  THEOREM 


The  Wigner-Eckart  theorem  states  that  if  we  have  a  spherical  tensor  Tkq  in 
the  space  spanned  by  the  wave  functions  |JM>,  then  the  matrix  elements  are 

CJ’M'lT,  |JM>  -  <J(M)k(q)| J'(M')>  <J*|T.  |j>  (Rose,  1957)  or 

kq*  k'  (3.D 

<J'M'  |T  I  JM>  -  (_Jj»  q  m)C J*  lTkl J)  (Judd.  1 963 ;  Wybourne,  1965). 

The  projection  (q)  dependence  is  contained  in  the  C-G  coefficients,  and  the 
factors  <J'|Tk|j>  are  called  the  reduced  matrix  elements. 

If  we  have  a  mixed  spherical  tensor,  rank  k  and  projection  X  in  spin 
space,  and  rank  k  and  projection  q  in  orbital  space,  the  Wigner-Eckart  theorem 
then  is 


<L'M£S'M£|Tjq|LMLSMs>  -  <L( ML ) k ( q )  |  L *  ( ) >  <s(Mg)ic( X )  | S » (M£)>  <L'S'  |TKk|LS>  . 

(3.2) 

Since  the  C-G  coefficient  is  purely  a  geometrical  factor,  all  the  physics  is 
contained  in  the  reduced  matrix  element.  The  Wigner-Eckart  theorem  allows  the 
extraction  of  the  geometrical  factors  from  many  complicated  matrix  elements; 
it  also  serves  as  perhaps  the  main  motivation  for  the  development  or  Racah 
algebra  in  dealing  with  angular  momentum  states. 
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3. 1  A  Single  d  Electron  In  a  Crystal  Field 


As  an  example  of  the  use  of  the  Wigner-Eokart  theorem,  wo  consider  a 
problem  that  is  simple  in  tensor  algebra  but  rather  important  in  the  spectra 
of  impurity  ions  in  crystals:  the  case  of  a  single  3d  electron  in  an  axial 
crystal  field.  Such  a  system  could  be  the  doubly  ionized  scandium  ion,  Sc  , 
substituted  for  a  doubly  ionized  constituent  ion  of  approximately  the  same 
ionic  radius.  The  solid  could  be  hexagonal  with  nearest  neighbor  ions  located 
along  ,£+he  c-axis  and  the  off-axis  ions  too  distant  to  have  an  effect  on 
the  Sc  ion. 

2  + 

We  assume  that  the  remainder  of  the  electrons  on  Sc  are  replaced  by 
an  appropriate  spherical  potential.  The  wave  functions  for  the  system  are 
taken  as 


R3dY2m 


(3.3) 


and  as  indicated  we  ignore  any  effects  of  the  spin  in  our  approximation.  In 
general,  the  radial  function,  R^,  oan  be  calculated  by  a  numerical  technique 
such  as  Hartree-Fock.  The  Hamiltonian  for  the  problem  we  consider  is 


hCEF  "  A20r2c20(F)  +  Ai|or'J4Ci4oCr* ) 


(3.M 


(an  axial  crystal  field  is  defined  as  Ak(J  ■  0,  q  *  0)  and  we  shall  assume  that 


kO 


\o<rl<> 


with 


.  kN 
<r  > 


f 

Jo 


k  2  ,  .  2  . 
r  R-.(r)r  dr 
id 


(3.5) 


The  series  in  (3.*0  is  terminated  through  four-fold  fields  (C^0);  odd-k  terms, 
if  present,  are  omitted  from  the  problem.  The  matrix  elements  of  HCEp  are 
given  by 


am-  |HCEF|  Jlm> 


B20  <lm’lC2ol*m>  +  B4o  <l[n'  lc40l*m> 


(3.6) 
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Fop  the  S4  site  in  yttrium  aluminum  garnet  (YAG),  the  axial  point  oharge 
lattioe  sums  are 

*20  “  6355  . 

A4o  -  25,089  onf’/A4  . 


For  So2+ 


<r2>  -  1.372  A2, 
<r4>  -  4.053  A4, 

®20  *  871 9  om~  ^ , 
B4o  -  1 01 , 686  cm 


If 


the  site  were  cubic,  B20  “  0  and  B44  » 
By  the  Wlgner-Eckert  theorem  we  have, 


60,769. 


generally, 


am’  |ckq|£m>  -  <Jl(m)k(q) | Z(m»)>  <MCJ*>  (3.7) 


and  from  (1.26)  we  have 

j  Y£m» YkqYg.ra  dfl  “  [nf~]  <t(m)k(q)|A(m')>  a(0)k(0)|i(0)>  (3.8) 

where  we  have  rearranged  the  order  of  l  and  k  in  the  C-G  coefficients.  By 

using  the  relation  of  Ckq  and  Ykq  in  equation  (3.8),  we  have 

<Jlm'  |ckq | «,m>  -  a(ra)k(q)|  £(m')>  <£.(0)k(0)  |£(0)>  (3.9) 

and  from  (3.7)  we  have  the  important  result 


<H|Ck|)l>  -  <£(0)k(0)|H(0)>  .  (3.10) 


In  general 

rpo +1  T1/2 

<£'|ck|)i>  -  a(o)k (0) | n ' ( 0) >  . 
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Thus  using  equation  (3.9)  in  (3.6)  we  have 


<2m* |HCEFj2m>  -  B20  <2(m)2(0) j  2(m)>  <2(0)2( 0) | 2(0) > 

+  <2(m)4(0) | 2(m)>  <2(0)4(0) | 2(0)> 


(3.11) 


and  we  notice  from  the  symmetry  properties  of  the  C-G  coefficients  (eq  (2.7)) 
that  <2(-m)k(0) | 2(-m)>  -  <2(m)k(0) j 2(m)>,  so  that  the  states  with  negative 
projection  (m  <  0)  have  the  same  matrix  elements  as  those  with  positive  m. 
Notice,  also,  that  had  we  considered  terms  in  the  potential,  Ckq,  with  k  >  4, 
they  would  not  contribute  since  <2(m)k(0) | 2(m)>  »  0  for  k  >  4. 

The  C-G  coefficients  in  equation  (3.11)  can  be  found  in  Rotenberg  et 
al  (1969)  (the  relation  of  the  C-G  coefficients  and  3“j  symbols  is  given 
therein)  and  are 


<2(0)2(0)  | 2(0)>  -  -(2/7)1/2 
<2(0)4(0) | 2(0)>  -  (2/7)1/2  , 


<2(1)2(0)|2(1)> 
<2(  1  )4(0)| 2(1 )> 
<2( 2) 2( 0) | 2( 2) > 
<2( 2) 4( 0) | 2(2) > 


-  -(1/l4)1/2  , 

-2(2/7) 1 /2 
3 

-  -(2/7)1/2  , 

-(1/14)1/2 

- 3 -  ' 


It  is  perhaps  easier  to  obtain  the  above  C-G  coefficients  by  using  equation 

(2.10)  . 


Thus  the  energy  is  given  by 


Eo 

"  2B20/7  +  2B40/7  • 

m  »  0  , 

E*1 

"  B20/7  -  4B40/21  • 

m  »  ±1  , 

(3.12) 

E±2 

“  -2B20/7  +  B40/21  * 

m  "  ±2  . 
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It  should  be  noted  that  the  trace  Eq  +  2E,  +  2E2  vanishes.  This  is  a  general 
feature  of  spherioal  tensors.  Sinoe  the  diagonal  matrix  elements  are  given  by 


<Hm|Ckq|Ara>  -  a(m)k(0)|£(m)>  <£|Ck|£> 


then 


l 

l  <Hm|Ck0|im>  -  <*|ck|)l>  l  <H(m)k(0)  |  i,(m)>  . 

m--£,  m 


The  C-G  coefficient 

< 8. (m ) 0(  0)  |  £(m)>  -  1 

-  (-1)a‘ra  /2UT  <£(mK(-m)  |  0(  0)  > 
and 

1  /2 

<P,(m)k(0)  |  Jt(m)>  -  (-1 ) "*+m(§£^-)  <«.(m))t(-m)  |k(0)>  . 


Therefore, 


l  <£(m)k(0)  j  S,(m)>  -  l  |k(0)>  <£(mH(-m)  |  0(0)  > 

m«-£.  /2k +  1  m 

-  (2£+1)5k0  (3.13) 


»  0  (k  >0,  the  only  values  of  interest  here)  , 


I 


where  we  have  used  the  orthogonality  condition  given  in  equation  (2.4). 

The  result  given  in  equation  (3.13)  is  a  very  useful  check  on  the 
calculation  of  the  energy  matrices,  since  it  is  very  easy  to  make  a  mistake  in 
the  evaluation  of  the  C-G  coefficients. 

In  our  later  work  we  will  encounter  problems  where  the  matrix  ele¬ 
ments  of  the  Hamiltonian  <£m’|H|£m>  <•  0.  In  these  cases  we  have  a  set  of 

basis  functions  <6  (such  as  1 2m>  above)  and  we  assume 
m  1 
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From  Sohrodinger’s  equation  we  have 


HiJ/  ■  E\|i  or 

I  £*  H$  -  I  aE6 
m  m  “mm 

and  we  multiply  through  by  $*,  and  integrate  to  obtain 

T  a  H  ,  ■  a  ,E  or 

L  m  m’m  m* 
n 

I  am[H  .  -  E6,]  -  0  , 

tn  mm  nrm 
in 

which  is  the  secular  equation  for  determining  the  energy  levels, 
system. 


(3.1*0 


E,  of  a 


All  our  efforts  will  be  directed  toward  obtaining  equation  (3. 11*)  for 
many  electron  systems.  We  shall  use  the  methods  of  group  theory  and  other 
techniques  to  reduce  the  number  of  components  in  equation  (3.1*0  to  a  minimun. 

In  the  previous  example  the  matrix  elements  Hm,m  *  0  (m'  #  m)  and  the 
energy  levels  were  given  by  Hmm.  If  we  consider  a  tetragonal  crystal  field 
given  by 

HCEF  “  B20C20  +  BU0C40  +  B44^C44  +  C4-4^  (3.1f>) 

where  is  real,  and  -  B^,  then  the  only  matrix  element  different  from 
those  of  equation  (3.12)  is 


<2-2|hcef|22> 


/70 

21  °44 


(3.16) 


H-22  * 

and  the  secular  equation  for  these  states  (|2±2>)  is 


-  E 


with  H_2-2  m  ^22  Hg_2  *  *^~22 


[22  "  E 

H-22 

!2-2 

H-2-2 

l  ”  H-22’ 
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(3.17) 


Expanding  the  determinant  in  equation  (3.17)  gives 


(E 


E 


f 


and  from  equations  (3.12)  and  (3.16)  we  have  (assuming  >  0) 

E  _  -  2  .  i_  D 

*+  7  20  21  BH0  21  B44  ' 

E  .  _  2  B  _l_n  M  D 

E-  7  B20  .21  B40  21  B44  • 


(3.18) 


and  with  the  energy  levels  E0  and  E±1  given  in  equation  (3.12),  we  have  all 
the  energy  levels. 

The  wave  functions  corresponding  to  E±  are 

*+  -  —  [|22>  ♦  | 2-2>] 

/2 

<Ji  -  [  |  22>  -  |  2~2>]  . 

/2 

We  then  have  all  the  energy  levels  of  a  single  electron  in  a  tetragonal 
field.  In  all  cases  we  have  found  five  energy  levels,  which  is  the  number  of 
states  of  the  free  ion  (28,  +  1  -  5). 

An  important  result  can  be  obtained  for  the  tetragonal  crystal  field 
if  we  let  B20  -  0  and  -  5Bi(0//70.  This  is  the  limit-  of  a  cubic  field;  in 
this  limit,  from  equations  (3.16)  and  (3.10),  we  have 


E0  "  E+  “  2W7 


/21  . 


(3.19) 


The  doubly  degenerate  level  (E^.E*)  is  denoted  E  and  the  triply  degenerate 
level  (E±i,E_)  is  denoted  T2.  These  labels  are  for  the  irreducible  represen¬ 
tation  of  the  cubic  group  in  the  Mulligan  notation;  E  is  referred  to  as  T 3  and 


30 


Tj  as  T5  in  the  Bethe  notation.  The  differenoe  between  the  two  sets  of  energy 
levels  given  in  equation  (3.17)  is  frequently  referred  to  as  lODq  and  is 


lODq  -  1  OBjjq/21  or 
B^q  ■  21Dq  . 


(3.20) 


The  quantity  Dq  or  lODq  is  frequently  reported  as  an  experimentally  determined 
parameter  in  papers  on  optical  data  taken  on  transition-metal  ions  containing 
d  electrons  whether  in  cubic  sites  or  not.  The  relation  given  in  equation 
(3.20)  holds  for  the  many-eleotron  configuration  nd  .  The  various  coeffi¬ 
cients  such  as  Dq  which  are  commonly  used  for  other  symmetries  are  given  by 
Konig  and  Kremer  (1977). 

3. 2  Problems 


1.  A  spherical  component  of  the  angular  momentum  of  an  eleotron, 
£  ,  has  the  matrix  elements 

ry  9 


<JLmf  1 8.^1  JLm>  -  <£(m)  1  (a)  |  Urn' ) >  <£|£|£> 


from  the  Wigner-Eckart  theorem.  Also  we  have  Hq |  Jtm>  -  m|£m>.  Using  a  table 
of  C-G  coefficients  (see  Rose,  1957,  appendix) ,  evaluate  <£(m)1  (0)  |  £(m)>  and 
then  show  that 

<£|£|£>  -  /FTThT  . 


2.  The  tensor  Tkq  has  the  property  T* 
element  q 


(-1) 


Consider  the  matrix 


<rm'|Tkq|lm> 


and  its  hermitian  conjugate 

(<£'m' |Tkq|Hm>)t  -  <lm|T»  |£’m'>  . 


By  using  the  Wigner-Eckart  theorem  on  the  above  matrix  element  and  its  hermi¬ 
tian  conjugate  show  that 


<£|Tk|d'>  -  (-I)8'-*'' 
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In  many  books  the  reduoed  matrix  elements  are  written  (A'jT^  |*)  with 


U|T  |4* )  -  /2TH  <A|T  \V> 


so  that 


(a|t  Jam  -  (-i)1"£,(a'|t  |A) 


solid. 


3.  An  electron  is  trapped  at  a  negative  ion  vacancy  site  in  a 
Taking  the  effective  potential  the  eleotron  sees  as 


B-Ac  +  co_?)  » 


H  “  B20C20  +  “22  ^“22  '  “2-2' 


calculate  the  energy  levels  of  the  p  state  (Yim(r))  of  the  electron.  (The  C-G 
coefficients  can  be  found  on  page  225  of  Rose,  1957).  You  will  need 


am|c20|im>  -  <A(m)2(0)|Jt(m)>  <A(0)2( 0)  | 1(0) >  , 


<"|c22h-i> 
<1  (  0) 2( 0) | 1 (0)> 
<1  (±1 ) 2 ( 0 ) 1 1 (±1 )> 


<1 (-0)2(2) | 1 ( 1  )>  <1  (0)2(0) | 1 ( 0) >  , 


/2 

5 


1 


<1  ( 1)2(2)  |  1  ( 1 )  >  - 


Ao 

£ 

5 


The  answer  can  be  obtained  from 

<io|h|io> 


2  B 

5  B20  ’ 


<1  ±1  |  H 1 1  ±1 > 
<11 | H | 1 -1 > 


5  B20  * 


5  B22 
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*4.  By  using  the  results  given  on  page  101  of  Ballhausen  (1962)  in 
equations  (3.12)  and  (3.16),  show  that 

B20  -  -70s  , 

Bjj0  -  21  (Dq  -  Dt)  , 

B|jij  “  \  /TODq  . 


Show  that  for  a  crystal -field  interaction  of  Cg  symmetry 


B^q  -  -l4Dq  -  21 Dt  , 

B43  -  2/70Dq  , 

as  given  on  page  104  of  Ballhausen  (note  that  the  term  involving  in  eq 
(3.15)  is  replaced  by  -  C^)). 
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4.  UNIT  SPHERICAL  TENSORS 
4.1  Discussion 


Because  of  the  power  of  the  Wlgner-Ekhart  theorem,  it  occurred  to 
Racah  to  cast  the  various  operators  representing  the  interactions  in  terms  of 
universal  quantities  that  could  be  tabulated  for  a  frequently  used  many- 
particle  system.  Toward  this  end,  Racah  introduced  the  unit  spherical  tensors 
for  the  electronic  configuration  n2,M,  which  we  define  as 

CA’m*  lukq| Am>  -  a(u)k(q)|jUm»  )>  «u, 


for  the  orbital  space  and 

<£'m's'm^|v^|amsma>  -  a(m)k(q)|  fc(m*  )>  <s(m3)icU)  |s(m^)>  6U,6SS,  , 


(4.1) 


for  orbital  and  spin  space. 


The  generalization  to  an  N-electron  system  is  simply 

Ukq  '  £  ukq(I)  ana 


(4.2) 


V<k  -  T  vKk(i) 

*q  f  AqU;  * 
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A  simple  and  often  used  example  of  these  tensors  in  orbital  space  is 


l  VU  '  \  <l|CK|l>  V0 

(4.3) 

-  <A|C|c|ft>  UkQ  . 

where 

<9.|Ck|il>  -  <A(0)k(0)  1 1(0) > 


(we  omit  the  upper  limit  on  the  i  sum  in  the  remainder  of  the  discussion). 
The  angular  momentum  is  simply  related  to  unit  tensors  by 


- I  yn  -  l  u1mu) 

Ly  -  U1y  ,  and  (4.4) 

-  AU+1 )  . 


An  example  of  a  tensor  in  a  mixed  spin  and  orbital  space  occurs  in 
the  hyper-fine  interaction  H5,  given  by 


-  (20BNuN/l)  l  yf/r^  ,  (4.5) 

where  0  is  the  Bohr  magneton,  0N  is  the  nuclear  magneton,  uN  is  the  nuclear 
moment ,  and  I  is  the  nuclear  spin.  Now 


1 


*1  '  Si  *  3?l(VJl)/r? 


(4.6) 


or 


N  (i)  -  l  (i)  -  /TO  l  <1  ( v)2(q-v) j 1  (q)>  Sv(i)C2  (i)  (4.7) 

M  M  ^  » W 

(we  show  in  sect.  6  how  eq  (4.7)  is  obtained  from  eq  (4.6)). 


36 


as 


The  part  of  Nq(i)  containing  4  (i)  can  be  written  in  terms  of  U1c[, 
in  the  second  part  of  equation  (4.4):  ** 


l  s  (i  )C0  n  (i)  -  <s  |s  |s><fi.  |C9|4>  V1  2  . 

j*  v  2,q-v  11  *2'  v,q-v 

A  component  of  i5  -  \  5  can  be  written 

i  1 

Nq  -  /4(4+1 )  U1q  -  /To  /s(s+1 )  <«,( 0) 2( 0)  | &( 0) > 

X  X  <1(v)2(q-v)|Kq)>  V^qfv  . 

Thus,  equation  (4.5)  can  be  written 


(4.8) 


(4.9) 


h5  -  (2eeNuN/i)  <i/r3>  l  n  i;  ,  (4.10) 

q  4  M 

with  Nq  given  by  equation  (4.9). 

Various  authors  use  different  normalizations  of  the  unit  spherical 
tensors.  The  relation  of  the  spherical  tensors  used  here  to  those  tabulated 
by  Nielson  and  Koster  (1963)  and  by  Polo  (1961)  are 


<L'S,a' |Uk|LSa> 


(L'Sa'|U(k)|LSa)  ~+1  - 
/2L' +1 


(4.11) 


<L*S’a' |V<k|LSa>  -  (L'S’a' IV^jLSa) 


<k  i 


✓4(24+1) 


C3(2L,+1)(2S'+1)] 


1/2 


(4.12) 


In  addition  Polo  tabulates 

(L'Sa* |C<k>|LSa) 


where 


CqK>  '  I  CKqU)  ’  (“-'3) 

Nielson  and  Koster  (1963)  tabulate  the  reduced  matrix  elements  of  V11  only; 
the  reduced  matrix  elements  for  the  electronic  nd”  configuration  of  V  ,  V1^, 
and  V ' 4  are  calculated  by  Wai-Kee  Li  (1971). 


37 


vmw\M  \JMif*  W*  V*  iTa  V*  k" M  ILM  1Ui  FJUSJi T A/1AA 


H.2  Bibliography  and  References 


Li,  Wai-Kee  (1971),  Reduced  Matrix  Elements  of  V^1^,  and  ^  for  dn 

Configurations,  Atomic  Data  £,  263. 

Nielson,  C.  W.,  and  G.  F.  Koster  ( 1  96 3 ) ,  Spectroscopic  Coefficients  for  the 
pn,  dn,  and  fn  Configurations,  MIT  Press,  Cambridge,  MA. 

Polo,  S.  R.  (1961,  June  1),  Studies  on  Crystal  Field  Theory,  Volume  I — Text, 
Volume  II — Tables,  RCA  Laboratories,  under  contract  to  Electronics  Re¬ 
search  Directorate,  Air  Force  Cambridge  Research  Laboratories ,  Office  of 
Aerospace  Research,  contract  No.  AF  1 9( 604)-551<1 .  [Volume  II  gives  the 
date  as  June  1,  1961.] 

Slater,  J.  C.  (I960),  Quantum  Theory  of  Atomic  Structure,  Volume  II,  McGraw- 
Hill,  New  York,  NY,  appendix  26  and  chapter  22. 


5.  RACAH  COEFFICIENTS 

The  Raoah  coefficients  arise  in  the  coupling  of  three  angular  momenta 


(Rose,  1957,  p  107)  to  form  a  final  resultant, 
momenta,  we  consider  two  coupling  schemes: 

In  the 

coupling  of 

the  angular 

scheme  A: 

J, 

+  h  "  ^12  ’ 

Jl2  4  I3 

-  1  , 

(5.1 ) 

scheme  B: 

+  ^3  “  ^13  ’ 

J13  4  ?2 

• 

(5.2) 

Coupling  scheme  A  is  represented  by  the  wave  function 
|A>  -  l  <J1(m1)j2(m2)|j12(m1+m2)>  <j12(m1+m2)j3(m  ) |j(m)> 

m  rn  m 


m^  m2m3 


*  |  J  -j  m^  j 2^2  J  ^  * 


scheme  B  is  represented  by  the  wave  function 
|B>  -  l  <J1(m1)j3(m3) ij13(m1+ra3)>  <J13(m1+m3)J2(m2) |j(m)> 


m^  m2m3 


*  I J  ^  m  ^  j  2*^2  ^  3^3  ^  * 
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(5.3) 


(5.4) 


The  coupling  schemes  A  and  B  are  oonneoted  by  a  unitary  transformation 

|B>  -  l  <A|B>|A>  5  (5.5) 

A 

the  coefficients  of  the  unitary  transformation  are  determined  by  taking  the 
inner  produot  of  equation  (5.3)  with  equation  (5.4). 

We  define  the  Raoah  coefficients  as  follows: 


W  { J  y  j  ^  2J  i  3*^  3  ’  ^  J ) 


_ 1 _ 

[(2.j12+l)(2j13-vl)]1/2 


<a|b> 


Thus, 

[(2j12>l)(2j13+l)]1/2w(j2j12j13j3!j1j) 


(5.6) 


-  I  <J1  (m1  )j2(m2)|j12(tt.1+m2)>  <j12(m1+m2)j3(m-in..-m2)|j(m)> 
m1m2 

(5.7) 

*  <J1 (m1 )j3(m-m1-m2) |j13(m-m2)>  <J1 3(m-m2)j2(m2) |j (m) >  . 


The  following  equation  can  be  obtained  from  equation  (5.7): 

<J2(m2)j1  (m1  )|J12(m1+m2)>  <J1 2[m}  +m2)j3(m-m1  -m,,)  |  j  (m)> 

-  I  [(2j12+l)(2j13^l)]1/2w(j2j1JJ3!j12J13)  (5.8) 

J13 

*  <J1  (m1  )j3(m-m1 -m2) |J1 3(m-m2)>  <J2(m2)j1 3(m-m2) | j (m)>  , 
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which  is  a  relationship  used  often  in  our  analysis.  For  clarity  we  rewrite 
(5.8)  as 

<a(a)b(e)|e(«+3)>  <e(a+8)d(6-a-B)|c(6)> 


(5.9) 


l  /( 2f+1 ) (2e+1 )  W(abcdjef)  <b(8)d(6-a-B) |f (6-a)>  <a(a)f (6-«) j o (6) > 
f 


(Rose,  1957).  The  Racah  coefficient  is  related  to  the  symmetrized  M6-J" 
symbol  by  the  following  equation: 


Wtabodief)  -  (-)a*b*0*a{a  . 


(5.10) 


The  symmetry  of  the  "6-J"  symbol  is 


SJ1  J2  J3 
Vi  *2  V 


(J2  J1  J3)  ^  pi  J3  J; 
V2  41  4)  Vi  i3  li 


1 J 1  *2  *3 

;4i  J2  J3. 


(5.11) 


and  all  combinations  of  the  relations  in  equation  (5.7).  The  four  triads 

(jl  J2  J3)*  ( J 1  a2  a3)»  Ul  J2  *3) »  and  Ui  *-2  J3)  must  b  able  to  form  a 

triangle.  That  is, 


I  J,  -  J2|  s  j3sj1  ♦  j2 


(5.12) 


with  similar  relations  for  the  other  triads. 

An  example  of  the  use  of  Racah  coefficients  is  in  the  calculation  of 
single-electron  matrix  elements  of  the  operator 


Vq 


l  <k(q-A)  1  ( A )  | k *  (q)>  , 


(5.13) 


which  arises  in  numerous  applications.  We  consider  the  matrix  element 


<rm'|Ek,q|am>  -  <l(m)k'  (q )  |  JL f  (m')>  <£'|Ek,|Jt> 


(5.1«0 
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by  applying  the  tfigner-Eokart  theorem,  equation  (3.1).  Also,  by  taking  the 
same  matrix  element  of  equation  (5.12)  we  have 


a'm’|Ek,q|Jlm>  -  I  <k  ( q-X)  1  ( A )  |  k '  (q)>  a’m’  Ick  JHm>  .  (5.15) 

Now  we  further  consider  the  matrix  element  in  equation  (5. 11*)  to  obtain 


<rm*|C.  .£.|4m>  -  I  a'm’jc.  . |£"m»>  <£"m" ] l.  | £m>  ,  (5.16) 
A  K,q"A 

where  we  have  used  matrix  algebra  on  the  product  of  two  operators.  If  we 
apply  the  Wigner-Eckart  theorem  to  the  last  matrix  element  in  equation  (5.16), 
we  obtain 


a"m"|fcJUm>  -  a(m)lU)|£'(m!')>  fi  <£|l|fc>  ; 


also,  m"  “  m  +  \  as  required  by  the  C-G  coefficient.  We  have  previously  3hown 
that 

a"|£|£>  -  /£<£+l)  Sln  .  (5.17) 

Therefore, 

<£"m"  |  |  £m>  -  <9.(m)  1  (>. )  |£"  (m-f-A)  >  6u,/i(l+l )  .  (5.18) 

Using  these  results  in  the  remaining  C-G  coefficient  in  equation  (5.16),  we 
have 

<£’m' |Ck  |£(m+A)>  -  <£(m+A)k(q-A)  |  £  •  (in' )>  <£’  |Ck|£>  .  (5.19) 

Substituting  the  result  of  equations  (5.19)  and  (5.18)  into  equation 
(5.16),  we  have 


<£ '  m*  |  C.  .  £.  |  £m>  -  A  (JM7<£  *  |C.  |£>  <£(m)1  (A)  I  £(m+A)> 

K,q"A  A  k  (5.20) 

x  <£(m+A)k(q-X)j£,(m')>  , 


giving  the  matrix  element  in  equation  (5.15).  If  we  substitute  the  result  of 
equation  (5.20)  into  equation  (5.15),  then  we  have 


<£'m'  |Ek,  |£m>  -  /£(£+1 )  <l'|Ck|£>S 
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(5.21 ) 


where 


S  -  I  <k(q-A)  1  ( A) |k*  (q)>  <2,(m)  1  ( A)  |  fc(m+A ) >  <Jl(ra+A)k(q-A)  |  V  (m* )  >  . 

A 

(5.22) 

The  last  two  C-G  coefficients  in  equation  (5.22)  can  be  recoupled  by  using 
equation  (5.8)  or 


<l(m)1  (  A)  j  H(m+A)>  a(m+A)k(q-A)|fc'(m'  )> 


(5.23) 


-■  l  J  ( 2f +1 )  ( 24+1 )  •  W  U  U  •  k }  l  f )  <l(A)k(q-A)|f(q)>  <l(m)f  (q)  |  T  (m* )  >  . 

f 

The  C-G  coefficients  in  equation  (5.23)  can  be  rearranged  by  using  the  sym¬ 
metry  rules  of  equation  (2.7)  to  give 


<1( A)k(q-A)|f(q)>  -  (-1)1+k“f  <k(q-A) 1 ( A) | f (q) >  .  (5.24) 


This  C-G  coefficient  and  the  first  C.-G  coefficient  in  equation  (5.22)  are  the 
only  two  C-G  coefficients  containing  A,  so  that 


l  <k (q-A) 1 ( A ) | k’ (q) >  <k(q-A) 1 ( A) | f (q)>  -  6  (5.25) 

A  IK 

because  of  the  orthogonality,  as  shown  in  equation  (2.4),  of  the  C-G  coeffi¬ 
cients.  Thus,  we  get 


S  -  (-i)1+k“kV(2k’+1)(2iM)  WUU'kjik')  a(m)k '  (q)  j *  <n?' )  >  ,  (5.26) 


which  when  substituted  into  equation  (5.21)  gives 


<rm'|Ek,q|£m> 


(-1)1+k“k,A^l)(2i,Tl)(2k’+1)  WdU’kjftk*)  <2,' 

x  <f.(m)k '  (q)  1 2, 1  (raf )  >  . 


(5.27) 


Upon  comparing  the  result  given  in  equation  (5.14)  with  equation  (5.27),  we 
have 


<t'|Ek,|l> 


<t’|Ck|»>  WUU’kilk’)  , 

(5.28) 
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which  is  a  useful  relation  if  we  wished  to  express  the  tensor  Ek)(.  in  terms  of 
unit  spherioal  tensors;  in  that  case  we  would  specialize  equation  (5.28) 
to  i'  -a  and  simply  replace  Ckg  in  equation  (M.3)  by  Ej.g  with  the  reduced 
matrix  element  given  by  equation  (5.28).  We  shall  have  frequent  ocoasion  to 
express  our  results  in  terms  of  Raoah  coefficients  by  using  equation  (5.8), 

5. 1  Problems 

It  is  frequently  convenient  to  build  spherioal  tensors  of  higher  rank 
by  coupling  produots  of  angular  momentum  operators.  One  such  tensor  is 


Tkq  "  2  <KcOl(q-a)|k(q)>  *a*q_o  ,  (a) 


where 


(a  a 

1  a  q* 


r 


-q+a 


(-!> 


> 


in  which,  by  the  properties  of  the  C-G  coefficient,  k  is  restricted  to  0  £  k  S 
2.  The  Tk(,  thus  constructed  is  patently  a  spherical  tensor  of  rank  k,  projec¬ 
tion  q.  Tne  application  of  the  Wigner-Eckart  theorem  (eq  (3.1))  gives 


<am’|Tkq|tm>  »  <a(m)k(q)|a(m,)>  <«. |tr |a>  (b) 


with  m’  -  m  +  q. 

Direct  calculation  of  the  matrix  elements  gives 


<am'|V.  |am>  -  l  <1(o)1(q-o)|k(q)>  <am'|a  an_  |am>  . 

KQ  a  Q“0t 

a  ^ 


(o) 


By  using  the  rules  for  matrix  multiplication  we  have 


<am’U  a  |arn>  -  T  <am’|a  |a"m">  <a"m"|a„  |am>  .  (d) 

04  q_a  a'^m"  °  q"a 

Using  the  Wigner-Eckert  theorem  on  the  first  matrix  element  in  (d)  gives 


<am,|aa|a',m">  -  <a"(m")i(a)|a(m,)>  <a"|a|a> 


(e) 


swv/rwuwwuwujft/  «tu 


But,  as  we  have  shown  earlier,  <8,"|l|i,>  -  /ITT+TT  and  the  C-G  ooeffioient 
requires  m"  -  m*  -  a.  Substituting  these  results  into  (d)  gives 

CAm’lfc  A  Ura>  -  <A(m' -a)  1  (a)  I  i.(m*  )>  <A(ra)  1  (q-a)  I  i.(m'  -a)  >  A(A+I)  (f) 

i  qj  Q—  Qt  1  1  ^  1 

where  we  have  used  the  same  technique  used  in  (e)  on  the  second  matrix  element 
in  (d). 


The  two  C-G  coefficients  in  (f)  oan  be  reooupled  using  equation  (5.8) 
in  the  text  to  give 

a(m)1  (q-a) |  Mm' -a) >  <t(m* -a)  1  (a)  |  J.(m'  )>  (g) 

-  I  /(2f+1  W21  +  1  )  WUUIjAf)  <1  (q-a)1  (a)  |?(q)>  <A(m)f  (q )  |  i(m' )  >  . 

f 

We  now  have  the  C-G  coefficients  uncoupled  so  that  if  we  consider  the  sum 
on  o  given  in  (a)  and  the  C-G  coefficient  in  (g) ,  we  have 


I  <1  ( a)  1  (q-a)  |k(q)>  <1  (q-a)  1  (a)  jf  (q)>  -  (-1)k5fk 
a 

where  we  have  used 

<1 (a) 1 (q-a) | k (q) >  -  (-1)k  <1 (q-a) 1 (a) |k (q) >  . 


Collecting  all  these  parts  together — (g),  (f),  and  (h) — and  substituting  into 
(c)  we  obtain 

<Am ’  |  T.  |  Am>  -  (-1  )kJl(  1+1  ) [( 2k  +  1 ) (2)1+1 ) l1  /2W( IU1  ; Ik)  <A(m)k (q)  j  £(m’ )  >  . 

Kq 

Comparing  this  result  to  (b)  we  obtain 


a|T  |A>  -  (-1)kil(a+1)C(2k+1)(2£+1)]1/2  W(tui|lk) 


(h) 


1.  From  the  results  obtained  above,  find  the  following: 


for  k  »  0 


<l|T0|l>  -  AU+1)  /2A+1  W(mi;ftO) 

-1 


WUU1  ;  S.0 ) 


✓3C2A+1 ) 


(Rose,  1957,  p  113) , 


<*|T0|!l>  - 


-&U  +  1  ) 
/3 
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But,  from  (a) 


T0  -  l  <1 (o)1 (-a) |0(0)>  l  1 
a  a  -a 


and 


< 1 (o) 1 (-o) | 0(0) >  -  (-1)1-0  —  <1(a)0(0)|l(o)> 

/3 


Substituting  the  above, 


(-1) 


1-ot 


✓3 


—  I  M)a!  l 

/3  a 


a  a 


also, 


Show  also 


--l(l)2  , 

/J 


(J)2|lm>  -  1(1+1 )|lm>  , 

<£m'  |  TQ j  im>  -  <!(m)0(0)|l(m)>  < ft | TQ 1 2, > 


-  -  —  ft( ft  +  1 )  . 

/3 


W(Jt111  ;11)  - 


[ft(ft+1 )  (2ft+1 )  j1  /2 


W(!1!1 } 12) 


✓(21*3) (21-1 ) 
3/30  /lU  +  1  )(21+1  ) 


1/2 


<l|T,  H>  .  .  /EEH  .  ,.|t  l.s  -  Ct(U1)(2>*3)(2t-1)] 

'  S2  2  3/5- 

2,  Given  the  results  of  problem  1,  calculate  the  reduced  matrix 
elements  of  the  tensor 


WKQ  “  E  <1(B  )k  (Q-8 )  |  K(Q  )>  laT, 


B  k ,Q-S 
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where  Tk^  is  given  in  the  discussion  at  the  beginning  of  seotion  5.1. 

The  reduoed  matrix  elements  of  W^q  are  given  by  the  Wigner-Eokart 
theorem  as 


am'  | WKQ | £m>  -  a(m)K(Q)|H(m')>  a|WK|H>  . 

Then,  by  computing  the  matrix  elements  of  WKq  by  the  methods  of  problem  1, 
show  that 

a|wK|4>  -  -(-i)k“Kuu+i)(?ui)(2K*i)]1/2  a|Tk|&>  wouaisno  . 


Evaluate  a|wk|)l>  for  K  -  0,  1,  2,  3.  (The  latter  Raoah  coefficient  is  given 
in  the  appendix  in  Rose,  1957.) 


Show  the  following: 


WUKJ.1  ;4K) 


K(K+1 ) 


4K(K+1  )  (2K+1  )i(t+1  H2J.+1 


'  1 
T 


/2 


k  -  K  , 


WUK+U1  ;ll,) 


WUK-U1  ;  IK) 


O  1  1  /? 

(K+2&+2)(K+1)  (2fc-K) 

4 (  2K+1  )  ( K+ 1 )  l ( I+T7 (2H+1 )_ 


?  T 1  /? 

(K+2&+1  )K  (2&+1-K) 

4(2K+1)K(2K-1)i(2t  +  1)U+1)_ 


k  -  K  +  1 


k  -  K  -  1  . 


We  write 


a|wK|i>  -  a|wK(k) |«>  . 


Then,  when  K  -  0,  we  have  k  -  1  and 


a|w0(i  )|£>  -  . 

°  /S" 

For  K  -  1 ,  we  have  k  -  0,  1,  2,  and 

<t|w.(l)|8>  -  -  [KKOUt-IHWjl 

2/3 

<Mw2(2)|»>--[tu,1)(2V)(2U3)]1/2  • 


47 


For  K  »  3#  we  have  k  -  2  and 


alw  (2>u>  -  [(t-1)t(t<-1)(H*2)(2t-1)(2t»3)3’/2 

3  3/re 
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6 


RACAH  ALGEBRA 


It  Is  convenient  in  many  vector  problems  to  express  the  vectors  in  terms 
of  spherical  bases  given  by 


Then 


Q±1  “  *(®x  ±  *®y )/&  » 


®0 


S*  -  C-1  )**•-„  . 


(6.1) 


®ux®v  *  “1^?  <1 (v) 1 (p) | 1 (p+v) >  ey+v  , 


(6.2) 


A  u  A 

®J*ev 


Jpv 


The  vector  t  can  be  written 


t  -  I 


e*A 

M  P 


t  -  l  e  A* 

;  m  »* 


(6.3) 


and 


A  -  e  , 
P  P 


x-s  -  z  a*b 

“  n 


p  P 


(6.4) 


I  A  B* 
P  P 
P 


Z  (-1 ) yA  B 
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Thus,  it*s  in  equation  (1.16)  oan  be  written 


bt  -  l  (6.5) 

so  that  the  spin-orbit  interaction  given  in  equation  (1.16)  is  immediately  in 
spherioal  tensors,  since  l  and  s  are  spherical  tensors.  That  is,  for  the 
manyeleotron  configuration1!  y 

H  -  C (r)  l  /aU+1  TsTsTmT  (-Dp  v”  .  (6.6) 

y  u 

An  example  of  Racah  algebra  is  the  reduction  of  the  latter  part  of  equation 
(M.6).  That  is,  we  wish  to  cast  the  interaction 


»  i-C-s  +  3n(r*s)] 

M  -  l  1*T  , 

q  q 


(6.7) 


where  we  have  dropped  the  subscript  i  on  the  components  and  introduced  the 
unit  vectors  r  -  r/r.  First  we  note  that 


f  -  i 


(6.8) 


as  in  equation  (6.3),  and  we  oan  write 


*  -+ 
r 


(6.9) 


Then  from  (6.7)  and  (6.8) 


T,  '  \  <-1>BC,_B(r,s8  . 

0 


(6.10) 


The  recoupling  given  in  equation  (1.27)  can  be  used  to  get 


C1qC1-B  "  E  <1  (0)  1  (0)  |k(0)>  <1(q)1(-®)|k(q-0)>  ,  (6.11) 
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and  since  <1 (0)1 (0) | 1 ( 0) >  -  0,  the  terms  in  (6.11)  are  restricted  to  k  -  0  and 
k  -  2.  Then  we  have 

C)qC^  -  <1  (0)  1  c 0)  1 0( 0)>  <1(q)1(-q)|0(0)>  66q 

(6.12) 

+  <1(0)1 (0)| 2(0)>  <1(q)1(-e)|2(q-0)>  C2>q_g 

From  the  relation  equation  (2.10)  we  have  <1 (q)0(0) | 1 (q)>  -  1  and  by  symmetry 
(eq  (2.7)) 

<1(q)1(-q)|0(0)>  -  (-1)1“q//3  (6.13) 


and 


C 1 Q° 1 "6 


(-1)q 


e.q 


~  (-1)B  <1 (B)2(q-B) | 1 (q)>  C 


2,q-B 


(6.HO 


where  we  have  used  the  symmetry  relation  of  equation  (2.7)  on  the  C-G  coef¬ 
ficient  < 1 (q) 1 (“3 ) I 2(q-B)> ;  also  <1 (0) 1 (0) \ 2(0)>  -  -(2/5)1/2.  The  first  term 
in  equation  (6.14),  when  substituted  into  equation  (6.10),  cancels  the  -sq 
term,  and  the  remainder  gives 

T  -  -  /To  l  <1(B)2(q-B)|l(q)>  s  C-  ,  (6.15) 

q  g  B  2,q  B 

which  is  the  form  used  in  equation  (4.7)  where  this  interaction  was  cast  into 
V12  tensors. 


As  a  further  application  of  Racah  algebra  and  some  of  the  other  material 
discussed  above,  we  shall  derive  the  gradient  formula  (Rose,  1957,  p  120).  A 
convenient  form  of  the  gradient  operator  is 


7 


(6.16) 


and  we  would  like 


grad  $(r )Ckq(r)  ■  [7,<j»(r )Ckq] 


(6.17) 
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First  we  observe  that 


rCkq  "  E  (_1 )X®-AC1ACkq 

-  xIk,  C-1 ) Xe_x  <1 (0)k(0) | k* (0) >  <1 ( A)k(q) |k' (q+A)>  Ck,>q+X  ,  (6.18) 


where  we  have  used  the  coupling  rule  for  spherical  harmonics,  equation  (1.27) 
(Rose,  1957,  p  61).  Now  we  write 


r*t  -  l  (-1)a+8e«eC,i  .  ,  (6.19) 

a8  a  B  1  a  B 


and  we  use  equation  (6.2)  to  eliminate  the  cross  product  to  produce 

rx£  -  -i/2  l  (-1)*  <  1  (X-ot)  1  (a)  I  1  ( X) >  e.C.  I  .  ,  (6.20) 

X.a  X  1-0  °"X 

where  we  have  replaced  the  sum  on  8  by  letting  8  -  A  -  a.  Now  in  calculating 
the  commutation  we  need  only  consider  the  operators  in  equation  (6.20);  thus, 
we  need 


[C1-ala-rCkq]  •  (6-21) 

Since  <|>(r)  commutes  with  Ci_at^_ai  we  need  not  consider  it  at  present.  Fir3t 
we  expand  the  commutator  to  obtain 


[Ci- 


a^a-A ’^kq 


^  “  C1-a*a-ACkq 


CkqC1-a*a-A 


(6.22) 


we  then  use 


^a-A^kq 


[l 

1  0“ 


A  ,Ckq^  + 


C.  I  . 
kq  a- A 


(6.23) 


in  equation  (6.22)  to  obtain 


[C1-aaa-A,CW 


C1-cJla-A,Ckq] 


C1-aCkq*ct-A 


ukql'1 


l 

-a  or 


(6.24) 
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The  last  two  terms  cancel  since  C,  and  C.  „  commute.  Thus,  we  obtain 

1  -a  kq 


^1-a^a-X  ,Ckq^  "  C1-a^a-X,CkqJ 

[C1-aVx*CkqJ  -  Ci^/kTkTTT  <k(q)1 (a-X)  |k(q+a-X)>  C(<>q+o.x 


(6.25) 

(6.26) 


where  we  have  used  equation  (2.19)  with'La_^  -  8,a_^  (which  are  identical  in 
the  commutation  brackets).  The  result  in  equation  (6.26)  is  not  quite  in  the 
form  we  want,  but  by  using  the  coupling  rule  for  spherical  harmonics  given  in 
equation  (1.27),  we  finally  obtain 


tC1-a4a-X  'Ckql  "  )  <k(q) 1 (a-X) |k(q+a-X)> 


(6.27) 


«  l  <k(0) 1 (0) | k" (0) >  <k(q-a-X)1(-a)|k"(q-X)>  C.„  . 

k  ,q  a 


In  equations  (6.16),  (6.17),  and  (6.20),  we  need 

[?»t.Ckq]  .  (6.28) 

We  can  see  from  equations  (6.27)  and  (6.20)  that,  when  this  is  formed,  the 
terms  dependent  on  a  are 


S  -  I  < 1 ( X-a) 1 (a)  1 1 ( X)>  <k(q) 1 (a-X) |k(q+a-X)> 
a 


x  <k(q+X-a) 1 (-a) | k" (q-X ) >  ; 


(6.29) 


that  is, 


[rx$,C.  ]  -  i  /2  l  (-l)Je  <k(0) 1(0) |k"(0)>  A(k  +  l)  SC.  „ 

kq  Aomk"  a  k  ,q  a 


(6.30) 


with  S  given  by  equation  (6.29).  The  sum,  S,  given  by  equation  (6.29)  can  be 
reduced.  First  we  write 

<k(q ) 1 (a-X) | k (q+a-X ) >  <k(q+a-X) 1 (-a) J k" (q-X )> 


(6.31) 

-  l  /(2f+1  )~(2k+1 )  W(k1k«1;kf)  <1  ( a-X  )  1  (-a)  | -(-X  ) >  <k(q )f  (-X )  |  k"  (q-X ) >  , 


D AMXf. 


where  we  have  used  equation  (5.8).  Thus,  the  sum  over  n  contains  the  terms 


l  <1 (a-\)1 ( —a ) 1 1 < — X) >  <1(a-X)1(-a)|f(-X)>  -  6f1  (6.32) 

a 

by  the  orthogonality  of  the  C-G  coefficients.  We  can  use  equations  (6.32)  and 
(6.31)  in  equation  (6.29)  to  obtain 

S  -  /3(2k+1 )  W(k1 k"1 5 kl )  <k(q)1(-x)|k"(q-X) >  .  (6.33) 


Using  the  results  of  equation  (6.33)  in  equation  (6. 30)  gives 

[rx£,C.  ]  -  i/6k(k+1)(2k+1 )  l  (-1)Xe  .  I  <k(0)1(0) |k( 0)>  W  <k1k"1;k1> 
Kq  X  k" 


(6.34) 


X  <k(q)1(X)|k"(q+X)>  Ck,I  q+x  , 


where  we  have  changed  the  sign  of  X  in  the  sum.  Multiplying  the  results  given 
in  equation  (6.34)  by  -i$(r)/r  and  combining  them  with  equation  (6.28),  we 
have  (changing  k"  to  k') 


[V , <j>(r )Ckq]  -  l  [I?  +  p  /6k(k+1  W2lc+1  J  W(k1k'1  ;k1)] 

X  k 

(6.35) 

X  <k(0) 1(0) |k*(0)>  <k(q) 1 (X) |k* (q+X)>  Ck?>q+X  . 

The  Racah  coefficients  in  equation  (6.35)  are  of  simple  form  and  are  given  by 
Rose  (1957,  p  227).  These  are 

Wtklk1 1  ; kl )  ■  -[^k.,i^|<.1)]1/2  .  K'-K*’  . 


(6.36) 


k'  -  k  -  1 


which  are  the  only  values  of  k'  allowed.  These  results  used  in  equation 
(6.18)  can  be  written  as 


! 

I 


[V,,<r)Ckq]  -  I  (-1 ) 


I 


k' 


<k(q  )1  (  X)  |  k r  (q+X )  >  Ck,fqfXDk  <J»(r) 
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(6.37) 


where 


k  +  1  , 


(6.37a) 


k'  - 


(6.37b) 


and  we  have  used  the  result 

<k ( 0) 1 (0) |k+1 (0)>  -  1 /2  and 


<k(0)  1  (0)  |  k-1  (0)>  -  -J-JL.]172 

(from  Rose,  1957,  p  225).  The  two  most,  common  forms  of  4>(r)  that  we  will 
encounter  are  rK  and  1/rK  1 .  For  *(r)  -  r  ,  we  obtain 

[vp,<ckq]  "  "  /k(2k+1)  <X(q)1(u)|k-1(q+u)>  rk_1  ,  (6.38) 

and  for  <j>  -  1/rl<+1, 


[V1/rl<+1ckq]  “  "  *^X  +  1 )  (2k  +  1 )  <k(q)1(y)|k+1  (q+p)  >  1/rk+2Ck+1  . 

(6.39) 


The  results  given  in  equation  (6.39)  are  easily  checked  for  k  -  0,  since  for  k 
-  0  we  have 

[Vu,1/r]  -  -<0(q) 1 ( p) | 1 (q+p) >  1/r2C1>q+y  ,  (6.39a) 


and  from  the  properties  of  the  C-G  coefficients,  v;e  know  that  c  -  0 
and  <0(0) 1  ( p) 1 1 ( p) >  -  1 .  Then 

^V1/rC1qJ  “  "Cl/r2  ’  (6.39b) 

Also,  we  know  from  vector  analysis  that 


grad  1/r  -  -r/r3  (6.39c) 
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and 


r  -  r  l  •  (6. 39d > 

Then  we  substitute  equation  (6.39d)  in  (6.39o)  to  obtain 

(grad  1/r)y  “  -C^/r2  ,  (6. 39e ) 

which  is  identical  with  the  result  of  equation  (6.39b).  We  shall  use  the 
result  given  in  equation  (6.35)  frequently  later  on,  particularly  in  the  form 
given  in  equations  (6.38)  and  (6.39). 

6.1  Problems 

1.  In  section  3.1  we  gonsidered  the  crystal-field  Hamiltonian  in  the 
form  Hqef  “  B20c2o(r)  *  BHoCi»o(r )  and  obtained  the  energy  splitting*  of  a 
single  d  electron.  All  the  previous  calculation  was  for  the  orbital  states 
only,  neglecting  the  spin-orbit  coupling.  For  the  states  |jm>  with 
J  -  t  ♦  3  we  showed  (in  sect.  2.1,  Problems)  that 

<j'm'As|Hg_o|jmAs>  -  . j 6m»m  2  “  s(s+1)]  , 

where  the  states  |jmfl,s>  are 


|jmAs>  -  J  <S,(p)s(m-u)  |  j  (m)>  1 4u>| sm-y>  . 
u 

For  these  same  states  show  that 


< j  *8,3  |ck  |jas>  -  (-1  )J'“'V(2j  +  1)(2jl+1  )  W(k?,j * s;AJ  )  <l(0)k(0)  j  1(0)  >  , 


and  evaluate  this  quantity  for 

l  -  2  , 

1 


and  the  quantity  {  }  ia  a  6-j  symbol  sbioh  oan  bo  found  in  Rotenberg  et  al 
(1969).  Using  the  above  results  oaloulate 


<J*m|rk0|jm>  -  <J (w)k( 0) | j * (m)>  <J 'As |Ck |jAs> 


and  obtain  the  following  table  for  the  matrix  elements  of  the  orystal  field. 


y 

m* 

J 

m 

°20 

B40 

1 

2 

1 

2 

3 

2 

1 

2 

1/5 

0 

5 

2 

1 

2 

5 

2 

1 

2 

3/35 

2/21 

3 

2 

1 

2 

5 

2 

1 

2 

-/S/35 

-2/S/21 

3 

2 

3 

2 

3 

2 

5 

2 

-1/5 

0 

5 

2 

3 

2 

5 

2 

3 

2 

2/35 

-1/7 

5 

2 

5 

2 

5 

2 

5 

2 

-2/7 

1/21 

3 

2 

3 

2 

5 

2 

3 

2 

-6/35 

2/21 

To  obtain  the  energy  levels  we  lot 


<2  llHCEFli  2*  "  H11  ' 


<-2  2 1 HCEF 1 2  2>  "  H22  ’  and 


<2  2 lHCEF^ 2  2->  H12 


(H21  -  Hi2^'  Then»  for  fc>le  m 


^  levels,  we  have  the  secular  equation 


Hn  *  Vo(J  ■  V  - E 


‘1  2 


12 


H22  *  Vo 


-  I)  -  E 
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with  a  similar  result  for  m 


•#.  What  Is  the  result  for  m 


5 

2 


? 


2.  In  the  consideration  of  the  nonrelativlstlo  limit  of  the  Dirao 
equation,  a  correction  to  the  nonrelativlstlo  Hamiltonian  arises  of  the  form 


H’ 


2  2 

°*  ao  fav  i_ 
>T“  ^r  Or 


for  terms  involving  the  orbital  motion  only. 

e2 

a  «  the  fine  structure  oonstant  - 


In  equation  (a) , 


(a) 


a  u  the  first  Bohr  radius  -  — r, 
o  c 

me 

V  “  the  potential  energy  (not  neoessarily  spherloally  symmetric),  $  -  [J,V], 
Take  V  -  fk(r )Ckq(f*)  and  show  that 


fk(r)  /k(k+1 )  l  (-1)A  <k(q)1(A)|k(q+A}>  C 


kq+A^-A 


Using  this  result  show  that 


H' 


-a 


2a§ 


3fk(r) 


9r 


Or 


Wr) 


fk(r) 


/k(k+1 )  l  < 1 ( —  A ) k (q+A ) | k(q)>  Ckq+xfc_A 
A 


Further,  show  that 

<A * m*  |  l  <1(-A)k(q+,.)|k(q)>  Ckq+xfc_x | im> 

A 

-  [£U+1  )(2l+1  )(2k+1 )  ]1/2WUU’k;tk)  <t(m)k  (q)  |  V  (m* )  >  <r|Ck|fc> 
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and 


a  a 

<t,m'|H*|i>  - - jjS.  FkUf(r)  <l(m)k(q)|fc'(m')‘>  <*'|Ck|.<t>  , 

3  )  1/9 

FkU,(r)  -  — |jr  ♦  -=5—  Ct(t+D(2l+i)k(k+i)(2k+l)],'iw(»ii»kjlk)  , 

r  “ 

H(«1fk,lk)  ■  ♦..KM),-  »’<*'♦'>  , 

JUU*1)(2l*i)k(k*1)(2k«1)]  '2 

3f  f 

FU9#  ,(r)  -  -Js  t-  +  i  [£(£+1)  ♦  k(k+1 )  -  I’U'+I)]  . 

**  3r  3r  2r 

If  S. *  •  !.  and  f,  -  C  rK  (C  i3  a  constant),  show  that 
k  o  o 

J  r2  <*•  \e\aiir)h  ■  0  • 

that  is, 

k-'  3  k(k-1)  fc-2 

'Kr  3r  2 

3.  In  the  basis  of  states  given  by 

|jm>  -  I  <a.(u)s(m-p)  |  j  (m)>  |i,p>|sm~p>  , 

y 

calculate  the  matrix  elements  of  and  the  matrix  elements  of  s^.  Show  that 
the  following  statements  are  true: 

<J'|a|J>  *  ("1  )J  )(2t+1  )(2J  +  1  )]1/2W(Uj's;)lj)  , 

<J  ’  IsJJj  >  -  [3(3+1  )(2s+1  )(2j+1  )]1/2W(1sjU;sj)  . 

4.  We  write  a  generalized  spin-orbit  interaction  as  the  mixed  tensor 

T  -  l  <1  (a)1  Cq-o)  |k(q)>  l  . 
a 
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and  we  wish  to  oaloulate  matrix  elements  of  In  the  coupled  spaoe  (single 
eleotron)  |jm>  with  J  -  £  ♦  s.  By  the  Wigner-eoKart  theorem, 


<J,m'|T«q|jm>  -  <J (m)k(q) | j ' (m’»  <j|Tk|j>  ,  (a) 

whloh  is  the  easy  part.  We  now  write  the  state  |jm>  as 

|jm>  -  l  <i(y)s(m-y) |j (m) >  |fcy>|sm-y>  (b) 

y 


and  a  similar  expression  for  <j'm'|,  so  that  the  matrix  elements  given  in  (a) 
are 

CJ’m'lT  |Jm>  -  I  <1 (a)1 (q-a) |k(q)>  l  <£(y)s(m-y) |J (m)> 
q  a  yy ' 

*  <My')s(m,-y')|j,(m')>  <fty '  |  Jta|  Hy>  (o) 


x  <sm,-y'|s<,  |sm-y>  . 
q-a 

The  Wlgner-Eckart  theorem  can  be  used  on  the  last  two  matrix  elements  in  (o) 
to  give 


<iy '  |  fta|  £y>  -  <i(y)  1  (a)  |  #,(y  ’ ) >  /UTmT 

<sm,-y' |Sq_a|sm-y>  -  <s(m-y)1 (q-a) |s(m'-y,)>  /s(s+1 ) 
where  we  have  used  <fc|g,|l>  -  /TTT+TT  and  <s|s|s>  -  /s(s+1 )  . 


(d) 


The  first  C-G  coefficient  in  (d)  requires  y'+y+a,  and  the  second  C-G 
requires  m-y+q-oun’ -y ' ;  when  y'  is  used,  this  results  in  m'  -  m+q,  which 
agrees  with  the  restriction  on  m'  given  by  the  C-G  coefficient  in  (a).  Sub¬ 
stituting  the  results  of  (d)  into  (c)  gives 


<J,m'|Tk  |jm>  -  l  <1 (a)1 (q-a) |k (q) >  <t(y)s(m-y) |j (m)>  <fc(y+a)s(m+q+y-a) |j '(m')> 
4  ya 


<l(y)  1  (a)|t(ii+a)>  <s(m-y)  1  (y-a)  |s(m+q-y-a)  >  [s  (s+1  )t(t  *-1 )  j1  /2 


H 


5 


(e) 
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and  we  shall  refer  to  the  C-G  coefficients  by  a  number  referring  to  the  posi¬ 
tion  in  the  sum  in  equation  (e).  First  we  reoouple  H  and  3  using  equation 
(5.9)  to  give 

<t(y) 1 (a) | i(y+a)>  <l(y+a)s(m+q-y-a) | J  * (m* ) > 


-  (-1)  l  /Tsf +1 ) ( 2i+1 )  Wdij'sjtf)  <i(y)s(m+q-y-a) |f (m+q-a)> 
f 


x  <1(a)f (m+q-a)|j*(m')> 


(f) 


where  we  have  used  the  symmetry  of  the  C-G  on  M  to  reverse  the  order  of  the 
first  two  angular  momenta.  We  now  reoouple  5  in  equation  (e)  with  the  first 
C-G  in  equation  (f)  to  give 

<s(m-y)1 (q-o) |s(m+q-y-a)>  <i(y)s(m+q-y-a)|f (m+q-a)> 


(_1)t+s-f+1  £  /(2q+1  W'2s+lV  W( 1sfi;sg) 


(g) 


X  <s(ra-y)i(u) |g(m)>  <1 (q-a)g(m) | f  (m+q-a)> 


where  we  have  used  the  symmetry  of  the  C-G  to  rearrange  the  order  of  the 
angular  momentum  in  both  coefficients.  The  second  C-G  coefficient  in  (e)  and 
the  first  C-G  coefficient  in  (g)  are  the  only  coefficients  containing  y,  and 
the  sum  is 


l  <s(m-y)t(y)  |g(m)>  <i(y)s(ra-y) |j (m)>  -  (-1)^  ^  S6 


gj 


(h) 


Collecting  what  remains  of  (e)  using  the  results  of  (f),  (g),  and  (h),  we  have 


<J,m’ |Tkq|jm>  -  (-1)J[s(s+1)i(i+1)(23  +  1)(2j  +  1)(2i+1)] 


1/2 


*  l  (-1)f/2f+I  W( 1 ij ' s; if  )W( 1 sfi; sj ) 
f 


X  l  <1 (a>1 (q-a) | k (q) >  <1 (q-a)j (m) |f (m+q-a) > 


(e' ) 


x  <1 (a)f (m+q-a) |j ’ (m' ) >  . 
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We  now  reoouple  the  last  two  C-Q  ooeff Iolanta  in  (e»)  to  obtain 


<1 (q-a)j (m) | f (m+q-a)>  <1 (a)f (ra+q-u) | j • (m* ) > 


-  (-1)J  <j(m)1(q-a)|f  (<n+q-a)>  <f (m+q-a) 1 (a) | J ' (m* ) > 

(i) 

-  l  W2h+1)(2:M)  W(JIJ'ljfh)  <1  (q-a)l(a) jh(q)> 


X  <j(m)h(q)|j*(m')>  . 

Finally,  the  sun  over  a  with  the  remaining  C-G  in  (e’)  and  the  first  C-G  in 
(i)  gives 


l  <1 (a) 1 (q-a) |k(q)>  <1 (q-a) 1 (a) |h(q)>  -  (-1)k6 


hk 


(J) 


substituting  into  (5a)  produoes 


<  j 'm'  |  Tkq |  Jm>  -  (-1  )k+^ '  [s(s+1 )  (2s+1  H(  Jl+1 )  (2JI+1 )  (2J+1 )  (2k+1 )  D1 /2 


x  <j(m)k(q) )j’m,)>  S 


where 


(e") 


s  -  I  (-1)f(2f+1)  W(Hj's;lf>  WdsftjsJ)  W(J1JM;fk)  . 
f 

From  Rose  (1957,  p  191)  or  from  Judd  (1963,  p  64),  we  have 

S  -  (-1  )JX(1  Ik;  ftsj  ;Hsj  * )  (k) 

where  X  is  the  X-ooefficient  and  is  identical  to  a  9-j  symbol. 

By  comparing  (o”)  with  equation  (a)  we  obtain 

<j'lTk|j>  -  M)k+J '+J  [s(s+1  )(2s+1  HOl-H  ) (2)1+1  )(2j+1  )(2k+1)]1/2 


(e"») 


x  X( 1 1 k;tsj ;tsj  f )  . 
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Prom  the  rules  for  rearranging  the  arguments  of  an  X-ooeff ioient ,  we  have  the 
above  equal  to 

XUsJitsj'  jllK) 

and  when  k  -  0  we  get  (Judd,  1  963 »  p  67?  Rose,  1957,  p  192) 


(-1  )^+^  ’^""8 

X(isj;isj'?110)  -  - ,-jz  W(isls;J1)  6  . , , 

[s(2J+1)]1/Z  J  J 

(see  also  Brink  and  Satohler,  1962,  p  119).  The  latter  reference  gives 


(1) 


therefore 


XUaj;*sJ?111)  -  0  ? 

<J,|T1|j>-0  . 

From  equations  (1)  and  (e’’’),  we  have 


,1/2 


<j •  |t0|j>  -  ) <a«» >  j’  w(u».,hi )  . 

/3 


From  Rose  (1957,  p  227),  we  have 


2[ft(l+1 )(2A+1 )s(s+1  )(2s+1 )]  /2 


and  from  equation  (m) 


<J'|T0|j> 


(-1) 


2J 


2/3 


UU+1)  +  s(s+1)  -  J(j+1 )]  . 


Now  when  k  *  0  in  Tkq  as  defined  at  the  beginning, 

T00  .  I  <1  (a)1  (q-a) |0(0)>  Vq-aSq,0 


,1-0 


l  — -  <1  (a)O(O)  1 1  (a)>  )las_a 

a  /3  ! 


(m) 


(n) 


(o) 


(P) 


J-  i.s 

/3 
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In  seotion  6  we  have  shown  that 


therefore 


£. a  -  1/2  [J2  -  t2  -  a2]  » 


<j'm'|T0|jm>  -  -  ~  CJ<J+D  -  1U*1)  -  s(s*1>] 


«*> 


and  in  equation  (o),  J-Jt±s-fc±  1/2,  2j  -  21.  ±  1  s  therefore, 


(-1)2J  -  -1  , 


so  that  the  result  given  in  equation  (o)  agrees  with  equation  (q)  but  is  muoh 
more  diffioult  to  obtain. 

We  have  an  interaction  between  two  electrons  given  by 


Tkq  -  I  <a(o)b(q-a) |ktq)>  Caq<1>Cbq..a<2)  . 

By  the  above  techniques  obtain  the  reduoed  matrix  element  in  the  expression 


<L‘MsjTkq|LM>  -  <L(M)k(q)|L'(MV)>  <L * | | L> 
where  the  states  are  given  by 


|LM>  -  l  <l(u)t(M-|i)jL(M)>  j£u;1  >|4M-v;2> 

M 

where  the  last  index  in  the  angled  brackets  refers  to  the  individual  elec¬ 
trons.  Hint:  The  result  should  have  something  like 

X(abk;££L;££L) 


6U 
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7.  FREE-ION  HAMILTONIAN  UNFILLED  CONFIGURATION  nS,N  [N  <  2(24,  +  1)] 

7.1  Background  for  Free  Ions 

The  approximations  made  in  the  analysis  of  the  spectra  of  ions  are 
not  new.  In  faot,  they  go  back  to  the  old  Bohr  orbit  theory.  Since  many 
readers  may  not  be  familiar  with  these  assumptions  and  may  not  remember  many 
of  the  concepts  and  most  of  the  technical  jargon  used  in  the  field  of  atomic 
spectra,  we  review  some  of  these  briefly.  We  stick  strictly  to  those  concepts 
which  apply  to  transition-metal  ions  and  rare-earth  ions. 

For  the  transition-metal  ions  we  shall  consider  the  doubly,  triply, 
and  quadruply  ionized  states.  For  the  rare-earth  ions  we  consider  only  the 
triply  ionized  oase. 

The  electronic  structure  of  the  three  series  of  doubly  ionized 
transition-metal  ions  is  given  in  table  7.1.  The  triply  ionized  rare-earth 
ions  are  characterized  by  the  electronic  structure  shown  in  table  7.2.  In  the 
rare-earth  series,  it  is  assumed  that  the  atomic  interactions  are  very  strong; 
thus,  when  an  ion  is  placed  in  a  crystal,  the  crystalline  electric  field  acts 
as  a  perturbation  on  the  ion.  In  the  transition-metal  ions,  the  electronic 
interaction  with  the  crystal  can  in  some  oases  be  larger  than  the  free-ion 
interaction.  Nevertheless,  we  shall  continue  to  label  the  states  using  the 
free-ion  labels.  This  assunption  allows  the  notation  developed  for  the  free 
ion  to  be  used,  with  the  reservation  that  many  of  the  "good”  quantum  numbers 
of  the  free  ion  are  not  quite  good  when  the  ion  is  present  in  the  crystal.  It 
is  assumed  that  the  free  ions  have  the  zeroth-order  Hamiltonian 


where  p.  is  the  momentum  of  the  ith  electron  and  u(rj)  is  an  appropriate 
spherical  average  potential  of  the  remaining  electrons  in  the  ion  (other  than 
the  N,  ntN).  The  single-electron  solutions  to  equation  (7.1)  are  taken  in  the 
form 


♦  *  Rnt(r)W?> 


(7.2) 
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TABLE  7.1.  ELECTRONIC  STRUCTURE  OF 
TRANSITION-METAL  IONS 


First 

series* 

Second 

series8 

Third 

series0 

Z 

■s 

BUM 

m 

mm 

Z 

X2+ 

5dN 

21 

So 

3d1 

39 

Y 

Hd1 

71 

Lu 

5d’ 

22 

Ti 

3d2 

HO 

Zr 

Hd2 

72 

Hf 

5d2 

23 

V 

3d3 

HI 

Nb 

Hd3 

73 

Ta 

5d3 

2H 

Cr 

3d“ 

H2 

MO 

Hd4 

7H 

W 

5d4 

25 

Mn 

3d5 

H3 

To 

Hd5 

75 

Re 

5d8 

26 

Fe 

3d6 

HH 

Ru 

Hd6 

76 

Os 

5d6 

27 

Co 

3d7 

H5 

Rh 

Hd7 

77 

Ir 

5d7 

28 

N1 

3d8 

H6 

Pd 

Hd8 

78 

Pt 

5d8 

29 

Cu 

3d9 

H7 

Ag 

Hd9 

79 

Au 

5d9 

30 

Zn 

3d10 

H8 

Cd 

Hd10 

80 

Hg 

5d10 

31 

Ga 

3d10Hs 

H9 

In 

Hd105s 

80 

Tl 

5d106s 

32 

Ge 

3d10Hs2 

50 

Sn 

Hd105s2 

81 

Pb 

5d106s2 

a(ls22s22p63s23p6 )3<?  *  (Ar  co reJSc? 

b(Ar  core)(  3d104a24p64d?  )  -  (Kr  core) 4c? 

c(Kr  core)(4d104f14Sa25p6 )S<?  -  (Lu3+  core)5<? 


TABLE  7.2.  ELECTRONIC  STRUCTURE  OF 
TRIPLY  IONIZED  RARE-EARTH  IONS 


Number 

Element 

Symbol 

Outermost  electron  shell 

57 

Lanthanum 

La 

Hd10Hf°5s25p6 

58 

Cerium 

Ce 

Hd^Hf’ss^p6 

59 

Praseodymium 

Pr 

Hd10Hf25s25p6 

60 

Neodymium 

Nd 

Hd10Hf35s25p6 

61 

Promethium 

Pm 

Hd10Hf45s26p6 

62 

Samarium 

Sra 

Hd10Hf55s25p6 

63 

Europium 

Eu 

Hd10Hf65s25p6 

6H 

Gadolinium 

Gd 

Hd10Hf75s25p6 

65 

Terbium 

Tb 

Hd10Hf85s25p6 

66 

Dysprosium 

Dy 

Hd10Hf95s25p6 

67 

Holalum 

Ho 

Hd10Hf105s25p6 

68 

Erbium 

Er 

Hd10Hfn5s25p6 

69 

Thulium 

T» 

Hd10Hf125s25p6 

70 

Ytterbium 

Yb 

Hd10Hf135s25p6 

71 

Lutetium 

Lu 

Hd,0Hf 1 45s25p8 

(Sohiff,  1968),  where  the  Y «ra(r)  are  the  spherical  harraonios  with  2  -  2  for 
the  transition-metal  ions  and  1-3  for  the  rare-earth  ions.  (Remember  that 
1-0  for  s,  2  «  1  for  p,  and  1-2  for  d  eleotrons.)  ^he  radial  funotions  in 
equation  (7.2)  are  taken  to  be  the  same  for  all  the  1  in  the  ion,  while  the 
angular  funotions,  along  with  the  spin  of  each  eleotron,  must  form  a  determi¬ 
nants!  function  so  as  to  obey  the  exclusion  principle.  Depending  on  whatever 
determinantal  function  is  chosen,  the  radial  funotions  oan  be  found  by  some 
self-consistent  method.  These  radial  funotions  (Freeman  and  Watson,  1962; 
Fraga  et  al,  1976;  Cowan  ?Jid  Griffin,  1  976)  have  been  found  for  the  Hund 
ground  states  of  all  the  transition- metal  find  rare-earth  ions  from  cerium 
(Ce3+)  through  ytterbium  Orb3+). 

The  Hund  ground  state  for  the  transition-metal  ions  with  N  S  5  and 
the  rare-earth  ions  with  N  %  7  is  given  by  assuming  that  all  N  spins  are 
parallel  and  that  each  angular  momentum  projection  is  the  maximum  allowed  by 
the  exclusion  principle  (in  eq  (7.2),  2  is  the  angular  momentum  and  m  is  its 
projection).  Thus,  the  Hund  ground  state  for  two  electrons  is  the  determi¬ 
nantal  function  (unnormalized) 


a(1)Yft  1(l)a(2)Yl  1.1 (2)  -  a( 1 >Yfc > £_, ( 1 ) a(2)Ya ^ j(2) 


where  a  is  the  spin  "up”  wave  function  (8  -  spin  down).  A  convenient  notation 
for  such  a  determinant  is 

{it, 2-1}  ,  (7.3) 

where  the  upper  sign  is  the  spin  projection  (+  -  up,  -  -  down)  and  2  and  2  -  1 
are  the  z  projection  of  the  angular  momentum  (m  in  e<£  (7.2)).  Thus,  the  Hund 
ground  state  equation  (7.3)  for  *»f2  praseodymium  (Pr3  )  has  total  spin,  S,  and 
total  angular  momentum,  L,  given  by 

S  -  1/2  +  V2  -  1  , 


L  -  2  +  2-1  -  22-1-5  (for  f,  2  -  3)  . 


Hence,  the  ground  state  is  L  -  5,  with  multiplicity  of  2S  +  1  -3.  In  the  so- 
called  Russell-Saunders  notation,  this  state  is  referred  to  as  3H,  as  given  by 
the  following: 

Total  angular  momentum,  L,  of  ion:  0123^567... 

Russell-Saunders  label  for  the  state:  SPDFGHIK.  .  . 

(continues  alphabetically) 
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In  this  notation,  the  technical  reference  to  auoh  a  state  is  a  term;  other 
terms  for  Pr3+  are  'I,  3P,  and  S  (Condon  and  Shortley,  1959;  Nielson  and 

Koster,  1963).  For  the  ion  Ce3  +  ,  which  has  one  f  electron,  the  atom  notation 

becomes  identical  to  that  of  the  ion,  that  is,  i  -  L  -  3  and  s  -  S  -  1/2  with 
the  single  term  2F.  Those  ions  in  the  series  for  N  >  2£  +  1  have  the  same 
terms  as  for  N  <  21  +  1 ,  and  their  Hund  states  can  be  constructed  simply  as 

d,  1-1,  1-2,  ...  21+1,  I,  1-1,  1-2,  ...  !-(p+1)}  (7.4) 

N 

where  the  number  of  electrons  is  N  -  21  +  1  +  p.  The  1  shell  is  completely 
filled  when  N  -  2(21  +1),  which  for  f  electrons  is  triply  ionized  lutetium 
(Lu3+).  As  an  example. of  equation  (7.4),  consider  triply  ionized  terbium 
(Tb3+),  which  has  the  4f8  configuration.  The  determinantal  wave  function  is 

+  +  +  ■♦•  +  +  - 

(3  2  1  0  -1  -2  -3  3)  ,  (7.5) 

where  total  angular  momentum  L  -  3  and  total  spin  angular  momentum  S  -  6/2  - 
3.  Thus,  the  Hund  ground  state  is  'F.  In  all  oases,  the  Hund  term  has  been 
found  to  have  the  lowest  energy  in  atomic  systems.  In  general,  the  wave 
functions  for  the  higher  terms  are  very  difficult  to  construct,  but  sophisti¬ 
cated  techniques  have  been  devised  for  the  orderly  development  of  a  set  of 
wave  functions  for  each  ion  having  the  electronic  configuration  p  ,  d  ,  and  fN 
(Nielson  and  Koster,  1963).  The  Hamiltonian  gi^pn  in  equation  (7.1)  has  the 
same  value  for  all  terms  of  the  configuration  nl  ;  consequently,  we  ignore  Hq 
in  the  future  discussion. 

7. 2  Significant  Free-Ion  Interactions 

A  number  of  interactions  within  the  ion  do  not  depend  on  the  partic¬ 
ular  solid  or  are  modified  when  the  ion  enters  a  solid.  These  interactions 
are  termed  the  free-ion  interactions .  We  discuss  several  such  interactions 
here. 


7.2.1  Coulomb  Interaction 


The  largest  contribution  to  the  Hamiltonian  for  the  free  ion  is  the 
electrostatic  interaction  of  the  n£N  electrons,  which  may  be  written 


where 


N 

I 

i>J 


(7.6) 
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The  matrix  elements  of  this  interaotlon  for  the  state  {3  2}  (the  term)  of 
Pr’+  are 


<3H|H1|3H>  -  F(0)  “  9F(2)  -  (l7F(4)/363)  -  (25F(6)/1M57)  (7.7) 


(Judd,  1963).  where  the  F^k'  are  frequently  referred  to  as  the  Slater 
parameters. 

The  F^  are  radial  expectation  values  given  by 


,(k)  2 


6  /o  Jo  vk+1  ^Rnit^r1  ^Rnt^r2^  r1  drir2  dr2  •  (7,8) 


where 


lo  Rnl(r)r2  dr 


— r  -  —  if  r,  <  r,  and 


^  If  r  >  r  . 

P  ,  1  J 


For  the  dN  electrons  the  matrix  elements  of  the  Coulomb  interaction  are  given 
in  terms  of  F(k|  ,  while  the  same  interaction  for  the  f”  series  is  given  in 
terms  of  new  parameters  Ek  by  Nielson  and  Koster  (1963).  Nielson  and  Koster 
give  the  matrix  elements  of  the  Coulomb  interaction  in  the  form 


<a'L'S» |H1 |aLS>  -  «LLI5SS,  l  c^Ca * .a. L, S)F 


(7.9) 


and  they  tabulate  the  coefficients  o(<(c,,a,L,S)  for  each  of  the  states  of 
d.  Similarly,  Nielson  and  Koster  give  for  the  matrix  elements  of  the  Coulomb 
Interaction 


<o ' L ' S  ’  | H 1  | aL5>>  -  «LL.«SS.  I  gk(a' .a,  L,  S)  Ek  , 


(7.10) 
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and  the  coefficients  g^a'.a.L.S)  are  given  for  aaoh  of  th<>  states  of  fN.  In 
equation  (7.9)  the  values  of  k  are  0,  2,  and  4,  while  in  equation  (7.10)  the  k 
values  are  0,  1,  2,  and  3.  The  relation  of  E*  to  F'*'  is  given  in  a  number  of 
places  (for*  example,  Judd,  1963,  p  206), 


7.2.2  Spin-Orbit  Interaction 

The  second  interaction  of  reasonable  magnitude  in  the  free  ion  is 
the  spin-orbit  coupling,  which  is 


-  1  c(r. )£. *s. 
i-1  111 


(7.11) 


where 


'<ri) 


•n2  i 


2m2o2ri  *1 


This  interaction  was  derived  from  relativity  theory  in  the  Bohr  orbit  quantum 
mechanics,  but  it  is  also  a  natural  consequence  of  a  nonrelativistic  approxi¬ 
mation  to  the  Di-ac  equation.  Values  of  and  £  (where  (■  -  <nfc|5(r )  |nfc>) 
from  Hartrse-Fock  wave  functions  are  given  in  tables  7.3  to  7.6.  In  the  rare- 
earth  series,  the  interaction  H2  is  quite  strong  and  is  in  general  much  larger 
than  the  interaction  of  the  rare-earth  electrons  with  the  crystal  fields. 


TABLE  7.3  HAfITREE-FOCK  VALUES  FOR  FREE-ION  PARAMETERS 
FOR  DIVAL.7NT  IONS  WITH  3dM  ELECTRONIC 
CONFIGURATION  (Fraga  et  al ,  1976) 


z 

X2* 

ndN 

F(2)  (cm"1) 

F(i°  (cm"1) 

C3d  (cm"1) 

<r2>  (A2) 

CrS  (A*1) 

21 

Sc 

3d1 

— 

— 

85.95 

0/346 

1 . 4997 

22 

Ti 

3d2 

67,932 

42,357 

131.15 

0.6716 

0.9808 

23 

V 

3d3 

74,062 

46,171 

187.17 

0.5677 

0.7112 

24 

C 

3^ 

79,790 

49,726 

265.60 

0.4910 

0.5401 

25 

Mn 

3db 

35,637 

53,368 

342.85 

0.4277 

0.4145 

26 

Fe 

3d6 

89,877 

55,927 

441.38 

0.3893 

0.3527 

27 

Co 

3d7 

94,600 

58,817 

561  „  21 

0.3525 

0.2949 

28 

N1 

3d8 

99,392 

61,756 

703.19 

0.3203 

0.2478 

20 

Cu 

3d* 

... 

— 

869.65 

0.2923 

0.2097 

71 


TABLE  7.4,  HARTREE-FOCK  VALUES  FOR  FREE-ION 
PARAMETERS  FOR  DIVALENT  IONS  WITH  ELECTRONIC 
CONFIGURATION  4dN  (Fraga  at  al,  1976) 


z 

X2+ 

ndN 

p(2) 

(om”1 ) 

p(  4) 

( om" 1 ) 

«4e 
(om”1 ) 

<r2> 

(A2) 

<r^> 

(A14) 

39 

Y 

Ad1 

— 

— 

312.00 

1.5737 

4.4402 

40 

Zr 

4d2 

51,177 

33,321 

432.03 

1.2734 

2.89)4 

41 

Nb 

4d3 

55,682 

36,328 

566.11 

1.0769 

2.0761 

42 

Mo 

4d“ 

59,873 

39,117 

713.12 

0,9316 

1 .5580 

43 

Ta 

4d5 

64,052 

41,911 

891.17 

0.8145 

1 .1907 

44 

Fu 

4d6 

67,247 

43,978 

1081.70 

0.7365 

0.9869 

45 

Rh 

4d7 

70,673 

46,224 

1299.11 

0  .  6656 

0.8126 

46 

Pd 

4d8 

74,108 

48,480 

1544.04 

0.6045 

0.6744 

47 

Ag 

4d9 

-- 

— 

1820.08 

0.5516 

C.5644 

TABLE  7.5.  HARTREE-FOCK  VALUES  FOR  FREE-ION 
PARAMETERS  FOR  DIVALENT  TONS  WITH  ELECTRONIC 
CONFIGURATION  5<T  (Fraga  at  al,  1976) 


Z 

X2* 

ndN 

p(2) 

(cm-'1 ) 

p(  4) 

(cm"1 ) 

45d 
(cm-1 ) 

<r2> 

(A2) 

<r1*> 

(A1*) 

71 

Lu 

Sd1 

— 

— 

1390.74 

1.6197 

4.6324 

72 

Hf 

5d2 

50,350 

33,000 

1773.59 

1.3646 

3.2437 

73 

Ta 

5d3 

54,008 

35,526 

2170,38 

1.1926 

2.4612 

74 

W 

5d4 

57.369 

37,840 

2594.40 

1 .0610 

1.9385 

75 

Re 

5d* 

60,702 

40,149 

3052.92 

0.9510 

1.5467 

76 

08 

5d6 

63,123 

41 ,766 

3530.75 

0.8779 

1.3277 

77 

Ir 

5d7 

65,755 

43,550 

4056.26 

0.8087 

1 .1289 

78 

Pt 

5d® 

68,386 

45,344 

4626.25 

0.7474 

0.9649 

79 

AU 

5d9 

-- 

-- 

5247.58 

0.6930 

0.6646 

72 


a-hjlsils  ilh 


TABLE  7.6.  NONRELATIVISTIC  HARTREE-FOCK  INTEGRALS  FOR  TRIPLY  IONIZED 
RARE-EARTH  IONS  (Fraga  et  *1,  1976) 


N 

H3* 

r<2'  (em*1) 

r(l,)  to*'1) 

Fifi)  (ora'M) 

C  (om"1) 

<r2>  (A2) 

<rlt>  (A4) 

<r6>  (A6) 

1 

C« 

— 

— 

— 

778.14 

1 .1722 

3.0818 

15.549 

2 

Pr 

105,120 

66,213 

47,718 

919.16 

1.0632 

2.5217 

1 1 . 492 

3 

Nd 

109,731 

69,165 

49,860 

1069.87 

0.97822 

2.1317 

8.9525 

4 

Pm 

113,640 

71,641 

51 .647 

1228.24 

0.91401 

1.8701 

7.4224 

5 

Sm 

117,222 

73,893 

53,269 

1397.79 

0.86059 

1 . 6700 

6.3365 

6 

Eu 

120,885 

76,204 

54,937 

1583.54 

0.81064 

1.4989 

5.3886 

7 

Gd 

124,644 

78,585 

56,655 

1786.68 

0.76368 

1.3267 

4.5589 

8 

Tb 

127,137' 

80,091 

57,722 

1990.51 

0.73523 

1 . 2508 

4.2673 

9 

Dy 

129,960 

81,829 

58,962 

2214.87 

0.70484 

1.1644 

3.902 

10 

Ho 

132,929 

83,670 

60,281 

2458.58 

0.67481 

1.0788 

3.5296 

11 

Er 

135,859 

85,486 

61 , 580 

2719.76 

0.64714 

1.0028 

3.2111 

12 

Tm 

138,754 

87,276 

62,864 

2999.22 

0.62154 

0.93514 

2.9355 

13 

Yb 

-- 

3299.82 

0.59678 

0.87030 

2.6713 

Consequently,  it  is  convenient  to  perform  all  the  calculations  in  a 
set  of  basis  functions  in  which  Hj  and  H2  are  diagonal.  The  set  of  functions 
that  achieves  this  is  the  totaj  angular  momentum  function  |JMj>,  where  the 
total  angular  momentum  operator  J  -  L  +  S. 

The  spin-orbit  interaction  H2  given  in  equation  (7.11)  commutes 
with  _the  total  angular  momentum  and,  consequently,  since  H1  also  commutes 
with  J2,  the  wave  functions  can  be  characterized  by  the  eigenvalues 
of  J2  and  J  .  That  is,  we  can  write  i|>TU  or  I  JM>  for  the  wave  functions  with 

32|JM>  -  J(J*1)|JM>  and 


(7.12) 


JZ|JM>  -  MjJM>  . 


For  any  term  of  given  L  and  S  (eigenvalues  of  C2  and  §2),  the  values  of  J  are 
restricted  to 


|l-s|  sjs  |l  +  s|  . 
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Then  the  wave  functions  are  oustomarily  written  *t> 


or  | JMLS>i  and  we  have 


and 


£2|  JMI.S> 
$2| JMLS> 
3?i JMLS> 
Jz I JMLS> 


L( L+1 ) J  JMLS>  , 
S( S+1 ) | JMLS>  , 
J( J+1 ) | JMLS>  , 
M|JMLS>  , 


(7.13) 


CJ’M’L'S' |H1+H2|JMLS>  «  •  (7.14) 

As  Implied  in  equation  (7.14),  the  energy  H1  and  H2  is  independent  of  M,  or 
each  J  level  of  the  free  ion  is  2J+1-fold  degenerate.  The  matrix  elements  in 
equation  (7.14)  do  not  vanish  generally  for  L’  ■  L  ±  1  and  S'  -  S  ±  1;  thus,  L 
and  S  are  not  strictly  good  quantum  manbers.  Nevertheless,  the  energy  levels 
are  labeled  as  though  they  were,  as  in  the  Russell-Saunders  notation,  2S+1Lj. 
An  example  of  the  energy  levels  for  the  4f2  configuration  of  the  free  ion 
(Pr’+)  is  given  in  table  7.7.  Also  included  is  the  same  ion,  Pr^+,  in  the 
host  materials  lanthanum  trichloride  (LaClg)  and  lanthanum  trifluoride  (LaFg). 


The  results  in  table  7.7  are 
Interesting  in  that  they  show  that  most  of 
the  energy  levels  observed  in  the  free  ion 
are  lowered  when  the  ion  is  embedded  in  a 
solid.  This  shift  in  the  energy  levels  is  a 
general  effect  and  is  not  restricted  to 
Pr3+,  but  exists  in  all  the  rare-earth  ions 
where  a  comparison  with  energy  levels  of  the 
free  ion  oan  be  made.  In  fact,  this  shift 
has  been  observed  by  Low  (1958a,b)  in  ions 
with  an  unfilled  d  shell,  The  first  explan¬ 
ation  of  this  shift  in  energies  was  by 
Morrison  et  al  (1967).  who  showed  that  if 
the  ion  under  investigation  was  assumed  to 
be  embedded  in  a  solid  of  homogeneous  die¬ 
lectric  constant,  e,  then  a  decrease  in  the 
Slater  integrals  is  given  by 

AF(k)«  -e2(e-1)(k*1)<rk>2/{b2k+1[e+k(e+1)]}  , 

(7.15) 

where  b  is  the  radius  of  a  fictitious  cavity 
surrounding  the  rare-earth  ion.  The  result 


TABLE  7.7.  FREE-ION  ENERGY 
LEVELS  OF  TRIPLY  IONIZED 
PRASEODYMIUM  AND  CORRE¬ 
SPONDING  CENTROIDS  IN  TWO 
CRYSTALS  (all  in  cm"1) 


(Dleke,  1968) 


[LS]J 

Free 

LaCl3 

LaF3 

3Hj, 

0 

0 

0 

3h5 

2,152 

2,119 

2,163 

3«6 

4,389 

4,307 

4,287 

3f2 

4,997 

4,848 

5,015 

6,415 

6,248 

6,368 

3p4 

6,855 

6,684 

6,831 

’c* 

9,921 

9,704 

9,801 

'D2 

17,334 

16,640 

16,847 

3po 

21,390 

20,385 

20,727 

3pl 

22,007 

21 ,987 

21,314 

W 

22,212 

21,327 

— 

3p2 

23,160 

22,142 

22,546 

1so 

50,090 

48,710 

46,786 

given  in  equation  (7.15)  was  first  successfully  applied  to  Co++  in  MgAlgOi,. 
Later,  Newman  (1973)  showed  that  the  shift  In  F'k)  given  In  equation  (7.15) 
was  sufficiently  large  to  predlot  the  shifts  In  the  energy  levels  for  rare- 
earth  ions.  More  reoently,  Morrison  (1980)  derived  the  result 


AF 


(k) 


“iV 


<**»<rk>2 


(7.16) 


where  a.  Is  the  polarizability  of  the  Z,  ligands  at  Rj  and  <rk>  is  the  radial 
expectation  value  of  r  .  The  result  given  in  equation  (7.16)  is  believed  to 
be  more  fundamental  than  that  of  equation  (7.15)  because  the  latter  explicitly 
accounts  for  the  local  coordination  of  the  rare-earth  metal  ion.  Morrison 
(1980)  gives  a  predicted  shift  in  the  spin-orbit  parameter,  c,  but  beoause  of 
the  smallness  of  the  predicted  shift  and  the  errors  in  the  fitting  of  the 
experimental  data,  no  comparisons  were  made. 

Because  of  the  lack  of  experimental  data  on  the  free-ion  spectra  of 
rare-earth  ions,  measurement  of  the  shift  in  the  Slater  integrals  is  possible 
only  for  Pr^+.  The  experimental  F'k'  for  triply  ionized  rare  earths  in  LaClo 
have  been  obtained  by  Carnall  et  al  (1978),  and  these  results  are  given  in 
table  7.8,  and  the  corresponding  experimental  values  for  the  transition-metal 
ions  are  given  in  table  7.9.  These  data  can  be  used  in  conjunction  with 
equation  (7.16)  to  obtain  results  that  can  perhaps  be  applied  to  an  arbitrary 
host  material  to  predict  a  priori  the  energy  level  shift  of  that  host. 


TABLE  7.8.  FREE-ION  PARAMETERS  FOR  TRIPLY  IONIZED 
RARE-EARTH  IONS  IN  LaCl,  OBTAINED  FROM  FITTING 
EXPERIMENTAL  DATA  (all  in  dm  )  (Carnall  et  al,  1978) 
Values  in  parentheses  were  not  varied  in  the  fitting 


Ion 

F2 

F2* 

F6 

a 

6 

Y 

C 

Pr 

68,368 

50,008 

32,743 

22.9 

-674 

(1520) 

744 

Nd 

71 ,866 

52,132 

35,473 

22.1 

-650 

1586 

880 

Pm 

75,808 

54,348 

38,824 

21 .0 

-645 

1425 

1022 

Sm 

78,125 

56,809 

40,091 

21 .6 

-724 

(1700) 

1168 

Eu 

84,399 

60,343 

41,600 

16.8 

(-640) 

(1750) 

1331 

Gd 

85,200 

60,399 

44,847 

(19) 

(-643) 

1644 

(1513) 

Tb 

90,012 

64,327 

42,951 

17.5 

(-630) 

(1880) 

1707 

Dy 

92,750 

65,699 

45,549 

17.2 

-622 

1881 

1920 

Ho 

95,H66 

67,238 

46,724 

17.2 

-621 

2092 

2137 

Er 

98,203 

69,647 

49,087 

15.9 

-632 

(2017) 

2370 

75 


TABLE  7.9.  EXPERIMENTAL  VALUES  OF  FREE-ION 
PARAMETERS  (ors”1)  FOR  DIVALENT  3d"  ELECTRONIC 
CONFIGURATION  (Uylings  at  al,  1984) 


z 

ndN 

“177T- 

“Hn— 

*d 

a 

e 

22 

Ti 

3d2 

54,927 

32,206 

118 

20.52 

-466.2 

23 

V 

3d3 

59,924 

36,268 

170 

22.90 

-480.8 

24 

Cr 

3d4 

64,798 

40,288 

231 

25.83 

-509.2 

25 

Mn 

3d5 

69,485 

44,305 

316 

29.20 

-537.1 

26 

Fe 

3d6 

74,282 

48,241 

422 

33.21 

-533.2 

27 

Co 

3d7 

78,906 

52,227 

536 

37.48 

-532.48 

28 

Ni 

3d8 

83.514 

56,164 

668 

42.49 

-554.9 

7.2.3  Interoonf lguratlonal  Interaction 

An  interaction  that  has  been  frequently  used  in  fitting  the  ’'free” 
ion  levels  of  a  rare-earth  ion  or  a  transition-metal  ion  in  a  crystal  is  the 
so-called  interconf igurational  mixing  or  the  Trees  interaction.  For  the  rare- 
earth  ions  this  interaction  has  been  parametrized  by  WyboUrne  and  Rajnak 
(Wybourne,  1965)  and  is 


H1fJ  -  oL( L+1 )  +  BGtG,)  + 


tg(r?) 


(7.17) 


are  given  in 
TABLE  7.10.  EIGEN- 


where  a,  0,  and  Y  are  parameters  adjusted  to  fit  the  experimental  data.  The 
operator  g(g2)  is  the  Casimir  operator  for  the  group  G2.  and  c(Ry)  is  the 
similar  operator  for  Ry  (note  that  -  L(L+1)  is  the  Casimir’s  operator  for 
the  group  R^) .  The  values  for  these  operators  for  all  the  states  are  tabu¬ 
lated  by  Wybourne  ( 1 965,  p  73).  The  values  for  the  state  of  f2 
table  7.10.  The  values  of  a,  8,  and  Y  obtained  by 
fitting  experimental  data  for  the  rare-earth,  ions  are 
given  in  table  7.8.  To  my  knowledge,  no  successful 
attempts  to  derive  theoretical  values  of  a,  0,  and  Y 
have  been  published.  For  the  transition-metal  ions  the 
Trees  interaction  in  equation  (7.17)  uses  a,  but  8 
either  multiplies  g(Rc)  or  the  seniority  operator  Q 
(Wybourne,  1965). 

7.2.4.  Other  Interactions 

Many  other  Interactions  are  considered  in 
the  free  ion,  such  as  spin-other-orbit,  orbit-orbit, 
and  configuration  Interaction.  All  these,  to  a 
greater  or  lesser  extent,  improve  the  fit  of  theoretical  energy  levels  to  the 
experimental  data.  We  omit  these  interactions  from  further  discussion  since 


VALUES 

OF 

CASIMIR’ 

S 

OPERATORS 

OF 

FQjl  STATE 

State 

a 

136 

5Y 

3P 

lF 

2 

12 

5 

12 

6 

5 

?H 

S 

D 

0 

'I 

30 

12 

5 

0 

0 

0 

6 

14 

7 

20 

14 

7 

42 

14 

7 

76 


,  Hg.  and  H^q  give  a  sufficient  representation  of  the  free  ion  for  our  pur¬ 
poses  here.  However,  we  list  a  number  of  interactions  including  the  above 
which  have  been  considered  by  various  researoh  workers  (Wortman  et  al,  1973): 


H1 

h2 

H3 

H5 

h6 

h7 

h8 

H9 

H10 


the  Coulomb  interaction 
the  spin-orbit  interaction 
the  orystal-f ield  interaction 

the  interaction  with  a  magnetic  field  (Zeeman  effeot) 

the  hyperfine  interaction 

the  spin-spin  interaction 

the  nuclear  quadrupole  interaction 

the  spin-other-orbit  interaction 

the  orbit-orbit  interaction 

the  interoonf igurational  interaction 

tne  spin-orystal-f ield  interaction 


The  notation  listed  above  is  that  of  Judd  (1963),  with  a  few  obvious  addi¬ 
tions. 


7. 3  Summary 

We  have  considered  the  Coulomb  interaction,  H1 ,  and  the  spin-orbit 
interaction,  H2,  for  the  conf iguration  nlN  in  the  free  ion.  The  wave  func¬ 
tions  that  are  chosen  as  a  basis  for  diagonal ization  of  these  interactions  are 
| JMLS> ,  and  the  resulting  energy  levels  are  labeled  according  to  the  Russell- 
Saunders  notation  as  given  in  seotion  7.1.  This  same  notation  (plus  addi¬ 
tional  quantum  numbers)  is  used  for  describing  an  ion  in  a  crystal.  The 
values  of  <rk>  that  are  needed  in  equations  (7.15)  and  (7.16)  are  given  in 
table  7.6  for  the  rare-earth  ions,  and  tables  7.2,  7.3,  and  7.^  give  the 
corresponding  values  for  the  doubly  ionized  ndN  ions.  The  wave  functions  used 
for  the  calculation  of  the  energy  levels  of  rare-earth  and  transition-metal 
ions  in  a  solid  will  be  the  combination  that  simultaneously  diagonalizes  H1 
and  H2.  While  this  process  is  not  a  good  one  for  the  3dN  configuration,  it  is 
better  for  the  ^dN  and  5dN  configurations  and  is  excellent  for  the  triply 
ionized  rare-earth  ions. 


7.H  Problems 


1.  We  have  the  tensor  operator  Tkq  given  by 


N 

I 

1-1 


Ckq(1) 


I 
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(a)  For  N  •  2  evaluate  the  matrix  elements  or  TkQ  by  using  the  states 
given  by  equation  (7.3).  That  is,  show  that  q 


{£,t-i}nkq{i,a-n  -  [<u|ckq|u>  +  <u-i|ckq|u-i>]«q0  . 


(b)  By  application  of  the  Wigner-Eoksrt  theorem  to  the  problem  in  (a) 

we  have 


)*Tkq{t.t-1 )  -  <L( L)k (0) | L( L) >  <LS|Tk|LS> 


with  L  -  28.  -  1. 


Using  this  result  and  the  result  in  (a)  show  that 
<I-s|Tk|LS>  ■■  '<X( L)k(  0  j I L( L)  >  !*'£)>  +  <l(  i-1)k(0)  |  l(£-1 )  >] 


where 


and 


<L(  L)k( 0) |  L( L)> 
<2(Jl)k (0)  |  fc(fc)  > 


f  ( 2L+1 ) (2L) 1 (2L)  1  11/2 

[( 2L+k  +1  ) ( 2L+k ) ! ( 2L-k ) ! J 


as  above  with  L  -  l. 


Show  also  that 

au-i)k(0)|tu-i)>  -  } ( ( 2t-k ) \ 2l~ 1 ) "] 1  /2 [ 2&-k-k2 ]  . 


2.  By  extending  the  results  obtained  in  problem  1  to  the  operator 

Tkq  ’  ^  Ckq(i)  * 

show  that 


N 

<LLSS|Tk0|LLSS>  -  £  <lU~p)k (0)  |  tU-p)>  <i(0)k(0)  1 1(0)>  , 

P-0 

where  L  -  Nt  -  — ■ - ,  S  -  N  S  2.1  +  1. 
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3.  For  N  -  29,  *  1  what  ia  the  value  of  the  following  sum? 


24 

SM(k)  -  I  <4(4-p)k(0)|4(4-p)> 
P-0 


The  reoursion  relation 
[k(k+1 )  -  24(4+1)  +  2m2]  <4(m)k(0)  | 4(m) > 

-  -(4+m) (4-m+1 )  <4(m~1  )k (0)  |  4(ra-1 ) >  -  (4-m)  (J.+m+l )  <4(m+1  )k(0)  |  4(m+1  )> 


oan  be  used  to  reduce  the  number  cf  Clebsoh-Gordon  coefficients  in  problems  1 
and  2.  obtain  this  reoursion  formula  (see  Rose,  1957,  chapter  III). 


5.  For  <4(0)k(0) j 4(0) >  derive  the  recursion  formula  (eq  (2.10)) 


<4|Ck+2|4> 


k+1 

k+2 


[ 


( 24+k+2) (24“k) 
(24-k-i )(24+k+3) 


-,1/2 

-J  <Mck|a> 


for  even  k.  Thus  since 


<4  |CQ  |4>  -  1 


any  given  k<4|Ckj4>  can  be  evaluated  algebraicly. 
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8.  CRYSTAL-FIELD  INTERACTIONS— PHENOMENOLOGICAL  THEORY  OF  CRYSTAL  FIELDS 


it 


W 


8.1  Discussion 


In  the  presence  of  a  crystal  field  we  take  the  interaction  of  an  ion 
whose  electronic  configuration'  is  nAN  as 


H3  ‘  l 

kq  M 


(8.1) 


where  the  Ckq  are  unnormalized  spherical  harmonics  given  by 


Ckq (n )  -  [4u/(2k+1  )]1/2  Ykq(r)  . 

The  use  of  the  Ckq  in  expressions  for  electronic  interactions  (along  with 
other  shorthand  notation  that  we  will  not  use)  is  practically  universal.  The 
number  of  terms  in  equation  (8,1)  that  need  to  be  considered  is  limited  by  the 
symmetry  cf  the  site  occupied  by  the  particular  ion.  Also,  since  we  will  be 
discussing  only  the  n£N  conf iguration,  k  is  limited  to  values  of  4  for  2,-2 
and  6  for  2-3.  This  limitation  arises  because,  independent  of  the  basis 
chosen,  Individual  matrix  elements  of  Ck_  will  have  to  be  considered,  and 
these  are  such  that  <2,|ck|2>  »  0  if  k  >  21. 

For  our  purposes  in  the  discussion  of  the  nonvanishing  Bkq  of  equa¬ 
tion  (8.1),  it  is  sufficient  to  consider  a  single  electron;  thus,  the  Hamil¬ 
tonian  we  consider  is 


H 


3 


I 

kq 


BkVkq 


(8.2) 


Since  H^  must  be  hermitian,  we  have  the  property  that  Bk_q  -  (-1)%kq,  which 
is  the  same  as  the  Ckq  given  in  equation  (1.25).  The  basic  assumption  of 
group  theory  is  that  Ho  is  invariant  under  all  the  operations  o^  the  groups 
under  consideration.  That  is,  we  shall  assume  that 


°AH3 


81 


(8.3) 


whe^e  0.  is  any  operation  of  the  group.  Using  equations  (8.3)  and  (8.2)  we 
have 


°»H3  '  i  B8,°lCk,  • 
kq 


(8,.  4) 


where  we  assume  that  0.  operates  only  on  the  electron  coordinates.  The  Bkq 
contain  the  dependence  on  r,  but  since  r  is  invariant  for  all  the  32  poLnt 
groups,  there  is  no  loss  in  generality.  Thus,  we  need  only  consider  the  re¬ 
sults  of  OxCkq  for  all  0X  in  a  particular  group.  The  lowest  group  C1  contaLns 
the  identity  operator  only;  consequently  all  Bkq  are  allowed.  The  group 
contains  the  inversion  operation  I,  and  ICk~  -  (-1)  Ckq,  so  that  equations 
(8.3)  and  (8.4)  give  Bka  -  0  for  all  odd  k,  while  all  B^q  with  even  k  are 
allowed.  The  group  C2  contains  the  operation  C2,  a  rotation  about  the  z-e.xis 

by  180°;  thus  C2Ckq  -  (-1)qCkq.  Consequently  for  the  group  C2  all  k  values 
are  allowed  and  only  the  Bkq  with  even  q  exist. 

The  group  Cg  has  the  symmetry  operator  o^,  given  in  problem  2  of 
section  1 .  That  is, 

which,  for  even  k,  gives  the  same  nonvanir hing  Bkq  as  C2,  but  for  odd  k,  q 
must  be  odd,  giving  B^  1 ,  B^-j  t  B33  ('or  oc*d  k  <  5. 

We  next  consider  the  two  groups  and  S„.  The  symmetry  operators 
for  the  C14  group  are  (a  rotation  about  the  z-axis  by  90  ) ,  2*  an 

cj  01- c;'  (CJ  -  E). 


For  the  first  of  these  operations  and  using  table  1.1  we  have 


C4Ckq  -  e‘ 


''kq  ’ 


(8.5) 


and  no  new  restrictions  are  imposed  by  the  symmetry  operations  C2  and  cjj.  The 
nonvanishing  crystal-field  parameters  are 


Bkq  ■  0 


(jqj  <  k  and  k  is  any  positive  integer),  unless  q  “  0,  ±4,  ±8,  ±12, 


so  that 


—  ] 

For  the  group,  the  symmetry  operator  3^  can  be  written  *  IC^  , 


oC  m  (-i)k  e”il,q/2C 
S4Ckq  {  U  e  Lkq 


(8.6) 


Thus,  for  even  k,  S4  has  the  same  nonvanishing  Bkq  as  C^,  but  for  odd  k, 


Bkq  “  0 


unless  q  -  ±2,  ±6  ...  with  odd  k  and  |q|  S  k.  This  condition  gives  rise  to 
the  nonvanishing  Bkq  as  Bg2,  B^2,  By2,  B^§  ... 


As  a  final  example  we  consider  the  group  D2.  The  symmetry  operations 
for  this  group  are  C2(x),  C2(y),  and  C2(z).  From  problem  2  of  section  1  we 
have 


C2(x,Ckq 


(-nkc 


k-q 


C2(»)Ck, 


<-1>k4,,c 


k-q 


(8.7) 


°2<x)Ckq  *  '-'^kq 


The  last  relation  in  equation  (8.7)  requires  that  q  be  even;  from  this  result 
the  first  two  relations  are  Identical.  If  the  first  relation  in  equation 
(8.7)  is  substituted  into  equation  (8.4)  and  the  summing  index  q  is  changed  to 
-q,  we  have 


(-1  )kB, 


k-q 


(8.8) 


Since  B. 

k-q 


(  -1  )qB* 
kq 


(see  the  discussion  of  eq  (8.2)),  and  since  q  Is  even, 


<-DkB»q 


(8.9) 


The  result  given  in  equation  (8.9)  requires  that  Bkcl  be  real  for  even  k  and 
imaginary  for  odd  k;  also  Bk0  -  0  for  all  odd  k,  since  B^0  -  B#Q  in  general. 

Using  the  above  techniques  and  table  1.1  we  obtain  the  results  given 
in  tables  8.1  and  8.2,  which  are  sufficient  for  d  electrons;  for  f  electrons, 
however,  the  tables  must  be  extended  to  k  -  6.  For  the  cubic  groups  the 
operations  are  more  involved;  the  crystal-field  Hamiltonian  is  given  in  table 
8.3  for  k  S  6,  whicn  is  sufficient  for  the  rare  earths. 

Sometimes  in  the  0  group  the  (111)  cubic  axis  is  chosen  as  the  z- 
axis;  the  crystal  field  is  then  given  as 

HCEF  "  B40^C40  *  v^0/'7  ^C43  "  C4-3^ 

(8.10) 


4  B6o(C60  *  (c63  -  c6_3)  -  /ffTil?  (c66  ♦  C6_6)J  . 

It  may  be  of  interest  to  note  that  the  lowest  odd-k  term  in  the  crystal  field 
for  the  0  group  is  (from  Polo,  1961) 

B9)*[(C9ij  “  C9-iJ  "  y^_/1  ^  (Cgg  -  cg„g)] 
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TABLE  8.1.  ALLOWED  VALUES  OF  Bkq  FOR  POINT  GROUPS  3  THROUGH  15* 


No.b 

0roupc 

B10 

B 

11 

B20 

B 

22 

B30 

B 

31 

B 

32 

B 

33 

B40 

B 

42 

6 

44 

Re 

la 

Re 

In 

Re 

In 

Re 

ta 

Re 

Im 

R* 

Im 

Re 

In 

3 

c2 

X 

0 

0 

X 

X 

X 

X 

0 

0 

X 

X 

D 

0 

X 

X 

X 

X 

X 

i) 

C8 

0 

X 

X 

X 

X 

X 

0 

X 

X 

0 

0 

X 

X 

X 

X 

X 

X 

X 

5 

c2h 

0 

0 

0 

X 

X 

X 

0 

0 

0 

0 

0 

0 

X 

X 

X 

X 

X 

6 

°2 

0 

0 

0 

X 

X 

0 

0 

0 

0 

0 

X 

0 

D 

X 

X 

0 

X 

0 

7 

C2v 

X 

0 

0 

X 

X 

0 

X 

0 

0 

X 

0 

0 

0 

X 

X 

0 

X 

0 

8 

°2h 

0 

0 

0 

X 

X 

0 

0 

0 

0 

0 

0 

0 

0 

X 

X 

0 

X 

0 

9 

c* 

X 

0 

0 

X 

0 

0 

X 

0 

0 

0 

0 

0 

0 

X 

0 

0 

X 

X 

10 

Sit 

0 

0 

0 

X 

0 

0 

0 

0 

0 

X 

X 

0 

0 

X 

0 

0 

X 

X 

11 

C4h 

0 

0 

0 

X 

0 

0 

0 

0 

0 

0 

0 

0 

0 

X 

0 

0 

X 

X 

12 

d4 

0 

0 

0 

X 

0 

0 

0 

0 

0 

0 

0 

0 

0 

X 

0 

0 

X 

0 

13 

c«v 

X 

0 

0 

X 

0 

0 

X 

0 

0 

0 

0 

0 

0 

X 

0 

0 

X 

0 

14 

°3d 

0 

0 

0 

X 

0 

0 

0 

0 

0 

X 

0 

0 

0 

X 

0 

0 

X 

0 

15 

°4h 

0 

0 

0 

X 

0 

0 

0 

0 

0 

0 

0 

0 

0 

X 

0 

0 

X 

0 

*An  X  indicates  the  presence  of  Bk  and  a  0  its  absence .  teasing  B.  are  0. 

brhe  number  is  that  given  in  Koster  et  al  (1963).  ® 

cSchoenflies  notation.  For  the  relation  to  other  notations,  see  Kostnr  et  al  (1963). 


TABLE  8.2.  ALLOWED  VALUES  OF  Bk  FOR 
POINT  GROUPS  16  THROUGH  27^ 


No.b 

Group0 

B10 

B20 

b30 

B 

33 

B40 

b43 

Re 

Im 

Re 

Im 

16 

c3 

X 

X 

X 

X 

X 

X 

X 

X 

17 

c3l 

0 

X 

0 

0 

0 

X 

X 

X 

18 

»3 

0 

X 

0 

0 

X 

X 

X 

0 

19 

c3v 

X 

X 

X 

X 

0 

y 

X 

0 

20 

°3d 

0 

X 

0 

0 

0 

X 

X 

0 

21 

X 

X 

X 

0 

0 

X 

0 

0 

22 

c3h 

0 

X 

X 

X 

0 

X 

0 

0 

23 

c6h 

0 

X 

0 

0 

0 

X 

0 

0 

24 

D6 

0 

X 

0 

0 

0 

X 

0 

0 

25 

C6v 

X 

X 

X 

0 

0 

X 

0 

0 

26 

°3h 

0 

X 

X 

0 

0 

X 

0 

0 

27 

r6h 

0 

X 

0 

0 

0 

X 

0 

0 

aAn  X  indicates  the  presence  of  Bk  and  a  0  its 
absence,  teasing  B ^  are  0.  '  1 

bVhe  number  is  thit  given  in  Koster  et  al  (1963). 
cSchoenflies  notation.  For  the  relation  to  other 
notations ,  sea  Koster  et  al  (1963). 
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TABLE  8.3  ALLOWED  VALUES  OF  Bkq  (k  S  5) 
FOR  CUBIC  OROUPS  28  THROUGH  32*  *b 


No. 

Group 

hcef 

28 

T 

b32^c32  +  C3-2^  *  b40^c40 

±  /sTPT  (c*n  ♦  c^J] 

+  86otC50 

*  'fTTz  (c64  ♦  c6_„)] 

29 

Th 

Same  as  28  but  B^2  »  0 

30 

0* 

Same  as  29 

31 

Td 

Same  as  28 

32 

°h 

Same  as  30 

aIn  all  cases  the  signs  of  the  B4(J  and  Bg^  parameters  are 
correlated. 

°The  z -axis  is  chosen  along  the  (001)  cubic  axis  (see  dis¬ 
cussion  of  eq  (8.10)). 


A  method  of  checking  the  results  given  in  tables  8.1,  8.2,  and  8.3  is 
by  the  use  of  the  full  rotation  compatibility  tables  of  Koster  et  al  (1963). 
These  tables  can  be  used  to  determine  the  number  of  times  the  identity  repre¬ 
sentation  appears  in  a  given  Ckq.  Actually  the  method  can  be  used  to  find  the 
Ckq  that  form  basis  functions  for  the  identity  representation  and  can  be 
extended  to  other  representations  of  the  single  group.  For  a  given  k  in  Ck  , 
we  look  for  the  number  of  I^'s  corresponding  to  Dk ,  and  for  the  Ckq  we  use  the 
plus  sign  for  even  k  and  the  minus  sign  for  odd  k  (ICkq  -  (-l)kCkq). 

As  an  example  we  consider  the  C2  group.  Using  table  13  of  Koster  et 
al  (1963),  we  have  D*  -  +  2 r2,  or  one  Ti  which,  from  our  previous  work, 
corresponds  to  C^q.  Thus  we  have  the  parameter  B^q.  For  D2  there  are 
three  I^'s,  which  correspond  to  B2q  and  B22;  since  B?2  is  complex  there  are 
three  independent  constants. 

As  a  second  example  we  choose  the  group  C2v.  From  table  22  of  Koster 
et  al,  for  we  have  one  r1  which  corresponds  to  B^;  for  D2  we  have 
two  Ti's,  which  correspond  to  B™  and  ReB22,  Normally  there  would  be  three 
B2q,  but  from  our  previous  results  we  know  that  B22  is  real.  This  process  can 
be  repeated  for  all  the  point  groups;  the  results  are  summarized  in  table  8.4. 
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The  results  presented  In 
tables  8.1  through  8.3  were  given 
assuming  a  single  electron;  however, 
it  is  generally  assumed  that  the  ra¬ 
dial  dependence  whioh  is  contained  in 
is  the  same  for  all  N  electrons 
in  a  single  configuration  niN.  This 
assunptlor.  is  inherent  in  the 
crystal-field  interaction  given  in 
equation  (8.1). 

Despite  all  the  restrictions 
imposed  by  the  symmetry  operations 
used  above,  there  exists  still  one 
more  restriction  (possibly  more  in 
some  groups)  that  can  be  imposed  on 
the  B^q.  This  would  be  apparent  if 
we  were  to  consider  a  particular 
model  for  the  computation  of  the 
Bkq.  Any  such  model  would  be  based 
on^  a  coordinate  system  embedded  in 
the  crystal  and  could  be  used  to 
determine,  say,  Bk  for  a  particular 
k  and  q.  Now  if  we  wished  to  change 
to  a  second  coordinate  system  ob¬ 
tained  from  the  first  by  a  simple 
rotation  about  the  z-axis ,  we  would 
obtain,  say,  B£  ,  and  the  two  sets  of 
parameters  are  delated  by 

B’  -  elq(|,B.  (8.11) 

kq  kq 

where  <J>  is  the  angle  of  rotation. 

This  result,  equation 
(8.11),  shows  that  we  can  choose  <J> 
such  that  we  can  make  one  of  the 
real  and  positive.  To  show  this  we 
assume  that  B'  is  the  parameter  we 
would  like  to  cl  real  and  positive. 

First  assume 

BK0  -  a  -  iB  .  (8.12) 

Then,  from  equation  (4.10),  we  obtain 


TABLE  8.4.  NUMBER  OF  TIMES  IDENTITY 


REPRESENTATION  APPEARS  IN  Ckq 


No. 

Group 

The 

number  of 

r1 

•a 

for 

k  -  1 

2 

3 

4 

5 

6 

3 

C2 

1 

3 

3 

5 

5 

7 

4 

Cs 

2 

3 

4 

5 

6 

7 

5 

JZ 

CM 

O 

0 

3 

0 

5 

0 

7 

6 

°2 

0 

2 

1 

3 

2 

4 

7 

C2v 

1 

2 

2 

3 

3 

4 

8 

°2h 

0 

2 

0 

3 

0 

4 

9 

C4 

1 

1 

1 

3 

3 

3 

10 

S4 

0 

1 

2 

3 

2 

3 

11 

c4h 

0 

1 

0 

3 

0 

3 

12 

D4 

0 

1 

0 

2 

1 

2 

13 

c4v 

1 

1 

1 

2 

2 

2 

14 

D2d 

0 

1 

1 

2 

1 

2 

15 

D4h 

0 

1 

0 

2 

0 

2 

16 

C3 

1 

1 

3 

3 

3 

5 

17 

C31 

0 

1 

0 

3 

0 

5 

18 

°3 

0 

1 

1 

2 

1 

3 

19 

C3v 

1 

1 

2 

2 

2 

3 

20 

°3d 

0 

1 

0 

2 

0 

3 

21 

c6 

1 

1 

1 

1 

1 

3 

22 

c3h 

0 

1 

2 

1 

2 

3 

23 

c6h 

0 

0 

0 

0 

0 

3 

24 

d6 

0 

1 

0 

1 

0 

2 

25 

c6v 

1 

1 

1 

1 

2 

26 

D3h 

0 

1 

1 

1 

2 

27 

D6h 

0 

1 

0 

0 

2 

28 

T 

0 

0 

1 

0 

2 

29 

Th 

0 

0 

0 

0 

2 

30 

0 

0 

0 

0 

0 

1 

31 

To 

0 

0 

1 

0 

1 

32 

°h 

0 

0 

0 

1 

0 

1 

B'  -  (a  cos  Q4>  +  0  3in  Q$)  +  i(a  sin  Q<p  —  S  cos  Q$) 
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The  imaginary  part  of  vanishes  if 


tan  Q$ 


.  £ 
o  a 


or 


♦o  -5  *r 


where  p  -  0,  ±1 . 


(8.13) 


(8.110 


The  real  part  of  Bkq  oan 
of  p.  The  complete  set  of  bJq  is 
appropriate  $  to  give 


be  made  positive  for  an  appropriate  choice 
then  obtained  from  equation  (8.11)  with  the 


‘«o 

e 


(8.15) 


In  most  of  our  work,  the  lowest  B^_  for  which  k  is  even  and  q  #  0  has  been 
chosen  as  real  and  positive.  Thus,  for  C2  point  symmetry,  the  twofold 
crystal-field  interaction  is  written 


B20C20  +  B22(C22  +  C2-2) 


(8.16) 


with  B22  real  and  positive. 

This  reduction  of  the  number  of  phenomenological  crystal-field 
parameters  needed  to  fit  the  experimental  data  does  not  help  much  in  the  low- 
symmetry  point  groups  (point  groups  1  through  8)  and  in  fitting  the  rare-earth 
ions.  However,  for  the  ndN  ions  it  allows  point  groups  9  through  15  to  be 
fitted  with  the  same  set  of  phenomenological  Bkq  (B20,  B^q,  and  ReB^),  and 
point  groups  16  through  20  to  be  fitted  with  B2Q,  B^0,  and  ReB^.  Thus,  for 
the  computation  of  the  energy  levels,  point  groups  9  through  15  and  28  through 
32  have  B2Q,  B^q,  and  ReB^  (in  the  cubic  groups  B?0  -  0,  and  is  related 
to  Bjjq ;  see  table  8.3).  Similarly,  point  groups  1o  through  27  have  B20,  B^q, 
and  ReB^  (with  B^  ■  0  for  point  groups  21  to  27;  see  table  8.2).  These 
restrictions  greatly  reduce  the  computation. 

8.2  Problems 


1 . 


Using  the  Hamiltonian  in  equation  (8.2),  prove  the  statement  that 


k-q 


M>4B* 


kq 


by  assuming  that  H3  is  real. 

2.  In  the  point  group  S^,  the  tensors  C^i  and  form  a  basis 
for  r3  and  r^,  respectively  (table  25,  Koster  et  al,  1963).  From  table  26  of 
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Koster  et  al  we  have  Tj  *  rij  -  Tj  .  Show  that  the  tensors  formed  by  reooupling 

the  product  using  equation  (1.27)  form  a  basis  for  r1 . 

3.  Using  the  same  methods  as  in  problem  2»  determine  what  bases  are 

formed  by  C10C1_1  and  c10c11  (C10  is  a  basis  for  r2). 
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9.  MATRIX  ELEMENTS  OF  H3  IN  TOTAL  ANGULAR  MOMENTUM  STATES  FOR  THE  ELECTRONIC 
CONFIGURATION  ntN 

9.1  Dlaouaslon 

In  order  to  make  full  use  of  tabulated  data  in  our  calculations,  it 
is  necessary  to  make  some  modifications  in  equation  (8.1).  Nielson  and  Koster 
(1963)  nave  calculated  the  reduced  matrix  elements  of  the  unit  spherical 
tensors  introduced  by  Racah.  As  was  shown  in  equation  (4.3),  the  Ckq(i)  can 
be  written  in  terms  of  those  tensors  as 

C.  (i)  “  <*K|t>  uk(i)  (9.1) 

and  kq  k  q 

l  Ckq(i)  "  <»|Ck|t>  U<k)  ,  (9.2) 

where 

<a jck |£>  -  a(o)k(o)|A(o)>  . 

Thus  equation  (8.1)  may  be  written 

H3  -  l  B*q  <fc|Ck|Z>  U^k)  .  (9.3) 

The  matrix  elements  of  Hg  in  total  orbital  angular  momentum  states  can  now  De 
written 


(9.4) 


6SS'6MsM£ 


I  Bj$q  <l|Ck|ft>  <L(ML)k(q)|L’(M{j)>  <L’Sa’ |u (k)  |LSa> 

KQ 
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where 

1  /2 

<L,S«,|U(k>  |LSa>  -  [|£— -f]  (l’So*  |U(k)  |LSa)  (9.5) 


and  the  last  expression  in  parentheses  is  tabulated  in  Nielson  and  Koster 
(1963). 


Alternatively,  we  could  use  the  tables  of  Polo  (1961)  to  obtain  the 
matrix  elements  of  Hg  as  given  by 

<L'M^S,M^a,|H3|LMLS«> 


$S3»$ 


BkY‘ 


1)L,~ML(-m!  q  Ml)  (L'Sa’|C(k)|LSa)  (9.6) 


where 


is  a  3j  symbol  and 


c(k) 

q 


N 

l 


L'-1 


(i)  . 


The  result  given  in  equation  (9.6)  is  much  more  convenient  to  use  for  compu¬ 
tation  than  is  equation  ( 9. ^ ) . 


The  spin-orbit  energy,  H2,  can  be  written  as 


H2  -  c  /s(s+1  HU+1 )  l  (-1)X  vj  \x  . 

X 

The  matrix  elements  in  total  orbital  angular  momentum  states  are 

A’rML 


(9.7) 


CL'M^S'M^a'  |H2  |LMLSMsa>  -  c  /s(s+1  HU+1  )  (-1) 


<l(ml)ka)|  L'(M^)> 


*  <s(ms)1(-A)|S'(m^)>  <L’S’a*|V11  |LSa> 
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where 


1/2 

<L'S'a'  IV1 1  |LSq>  ■  -|  [^-!0(2S'.li]  (L'S ’"I*”  lLS“)  •  <9'8) 


and  the  reduced  matrix  element  in  parentheses  is  tabulated  by  Nielson  and 
Koster  (1963). 


The  matrix  elements  of  Hj  in  total  angular  momentum  states  J 
(J  -  t+S)  can  be  written  J 


<J'M'So,|H_|JMLS«> 

5  (9.9) 

-  I  B*  <l\C  |Z>  <J(M)k(q ) |  J '  (M*  )>  <J  •  L ' Sot '  |uk |  JLSot> 
kq  q  K' 

where 


<J  'L'Sot' |uk|JLSa> 


(-1 )L"L'+J,_J/2j+1  W(kLJ'S;L,J)(L,Sa' |uk|LSa)  (9.10) 


and  again  the  reduced  matrix  element  in  parentheses  is  tabulated  by  Nielson 
■and  Koster  (1963).  AI30,  the  quantity  <2. |Ck  | JL>  in  equation  (9.9)  is 


<2|Ck|Z>  -  <Jt( 0)k(0)  |  JL(0) >  . 


(9.11) 


The  matrix  elements  of  H2  in  total  angular  momentum  states,  J,  are  given  by 


CJ'M'L’S'ct'  |H2|  JMLSa> 


-  -c  /  2  (  ?  + 1  )  (21+1  )  W(S1J,L';S'L)(L’S'o,|V1  '  |LSa)6JJt(5MM, 


(9.12) 


where  the  reduced  matrix  elements  in  L.S  space  are  tabulated  by  Nielson  and 
Koster  (1963). 
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10.  GROUP  THEORETICAL  CONSIDERATIONS 
10.1  Discussion 


We  do  not  go  into  the  details  of  group  theory  here  but  discuss  the 
use  of  the  tables  presented  in  numerous  texts  on  group  theory.  An  excellent 
text  for  physicists  is  Tinkham  (1964).  For  our  discussion  here,  as  in  section 
8,  we  use  Koster  et  al  (1963).  This  reference  uses  the  Bethe  notation  for  the 
irreducible  representations  (1^);  the  relation  of  the  Bethe  notation  to  the 
Mulligan  notation  (Aj,  Bi ,  etc)  is  given  in  the  appendices  of  Griffith 
(1964).  In  our  discussion  we  use  the  Bethe  notation  for  all  the  single  groups 
and  both  notations  for  the  oubio  0  group.  In  the  double  groups  we  use  only 
the  Bethe  notation. 


As  our  first  example  we  consider  a  single  d  electron  in  a  crystal 
field  of  S4  symmetry.  The  crystal-field  interaction  oan  be  obtained  from 
table  8.1;  it  is 


H3  “  B20C20  + 


B52C32  + 


B32C3-2  +  B40C40  +  B44 ' 


(c44  +  C4-4^ 


4-4* 


(10.1) 


If  the  spin-orbit  interaction  is  small  and  the  other  configurations  remote,  we 
can  at  present  ignore  the  B term  and  assume  that  angular  wave  functions  are 
Yjh,  (|2m>),  with  -2  &  m  <  2.  With  these  assumptions  we  use  table  30  of  Koster 
et  al  (1963)  to  find  that 


d2  *  ri  +  2r2  +  r3  ♦  r„  ,  (10.2) 
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where  we  have  used  the  D+  table  since  l|2m>  -  (-1 )2| 2m> ;  that  is,  the  d  elec¬ 
trons  have  positive  parity.  -  The  wave  functions  for  ndN  also  have  positive 
parity  since  I | LSdN>  -  (- 1 ) 2N| LSdN> .  But  for  the  f  electrons,  l|LSfN>  - 
(-1  )™| LSf”>,  which  is  even  for  even  N  and  odd  for  odd  N.  Thus,  for  a  single 
f  electron  we  would  use  the  D3  entry  in  table  30  of  Koster  et  al.  The  result 
given  in  equation  (10.2)  shows  that  the  five  d  electron  states,  which  are 
degenerate  in  the  free  ion,  would  be  split  into  four  levels  in  the  presence  of 
the  crystal  field;  these  levels  are  shown  in  figure  10.1.  All  the  free-ion 
degeneracy  is  removed  except  for  the  rg,r4  degeneracy.  We  can  detect  this 
degeneracy  in  character  table  25  of  Koster  et  ai  (1563,  p  50)  by  observing 
that  the  product  of  the  characters  for  and  Tjj  gives  unity  for  all  the  group 
operations  of  S4.  We  can  also  check  this  product  by  directly  taking  the 
matrix  elements  of  Hg  using  wave  functions  that  transform  as  Fg  or  in  S4. 
To  find  the  wave  functions  which  transform  according  to  the  irreducible  repre¬ 
sentation,  we  use  the  operation  S4  on  the  states  |2m>.  That  is, 

S4|2m>  -  IC"1 | 2m> 

(10.3) 

=  e"iirm/2|2m> 


and  for  m  -  0 

S4|20>  -  | 20>  . 


Figure  10.1.  Hypothetical  splitting  of  a  single  d  electron  in  a  crystal 
field  of  S4  symmetry:  (a)  magnetic  dipole  and  (b)  eleotric  dipole. 


That  Is,  the  oharacter  is  1,  so  that  J 20>  transforms  as  r1  (see  table  25  of 
Koster  et  al).  Also,  since 


Sjj  j  2±2>  -  - 1 2±2>  , 


the  character  is  -1.  From  table  25  under  the  operation  S^,  we  find  that  the 
states  with  character  -1  are  Tg,  so  that  the  two  wave  functions  |22>  and  |2-2> 
transform  as  Similarly, 

Sjj  1 2±  1  >  -  Ti|2±1>  ; 


from  table  25  we  see  that  1 21  >  transforms  as  r_  and  1 2—  1  >  as  IV.  It  is  not 
difficult  to  show  that 


<21 |H3 | 21 >  -  <2-1 |H3|2-1> 


(10.4) 


which  shows  directly  that  the  energy  levels  for  and  are  degenerate. 

The  matrix  elements  <2m'|H3|2m>  of  the  crystal  field  are  given  by 
table  10.1,  where  the  results  given  for  r2  are  obtained  from  equation 
(3.18).  From  the  results  of  table  10.1,  we  obtain  the  energy  levels  given  in 
table  10.2  for  Tl3+  for  a  particular  choice  of  Bk_.  Table  10.2  also  gives  the 
energy  levels  when  the  spin-orbit  interaction  is  included. 


The  decomposition  of  the  free-ion  state  onto  the  symmetry  given 
in  equation  (10.2)  gives  the  dimensions  of  the  secular  determinant  which  has 
to  be  solved  to  determine  the  energy  levels  of  the  system.  In  the  case  con¬ 
sidered  in  table  10.1,  the  dimensions  of  the  secular  determinants  are  1  for 
r-| ,  1*3,  and  Tij,  and  2x2  for  Tg.  If  for  example  we  had  an  H  state  in  Sj, 
symmetry  (L  -  5)  of  d  ,  then,  from  table  30  of  Koster  et  al  (1963),  we  have 


D5  *  3ri  +  4r2  *  3r3  +  3r4 


(10.5) 


and  we  would  have  to  solve  a  3*3  determinant  for  T-j ,  a  4x4  determinant  for  r2, 
and  a  3*3  determinant  for  fg  or  r^. 
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TABLE  10.1.  MATRIX  ELEMENTS  OF 
CRYSTAL  FIELD  FOR  A  SINGLE  d 
ELECTRON  IN  S^  SYMMETRY 


Irreducible 

representation  Energy 


7  B20  +  7  BH0 

£  o  +  JL  B  +  /TO  R 
7  20  21  BH0  21  BHH 


2  D  +  _L  B  /70  n 
7  20  21  HO  21  BHH 


[3,H 


1  B  -  -i.  B 
7  20  21  HO 


TABLE  10.2.  ENERGY  LEVELS  FOR 
Ti3*  (3d)  IN  Sj,  SYMMETRY4 


O 

a 

k* 

c  - 

158  cm"1 

IR 

E  (cm-1 ) 

IR 

E  (om-1 ^ 

r1 

0 

r5,6 

0 

r2 

395 

r5,6 

0 

r2 

2931 

r7,8 

2938 

r3,H 

3720 

r5,6 

3662 

r7,8 

3827 

*b20  m  394.7  cm~*t 
cm~  ,  B44  -  3152  cm"'* 


“40 


-  -7932 


These  are 


approximate  crystal-field  param¬ 
eters  for  Ti3+  in  the  Ga  site  of 
Gd jSc 26*3® 12 •  Gadolinium  scandium 
gallium  garnet . 

IR  m  irreducible  representation 
E  -  energy 


If  we  wish  to  determine  the  energy  levels  experimentally,  it  is  in¬ 
formative  to  investigate  the  use  of  polarized  radiation.  The  magnetic  dipole 
operator  is 


Hmd  ‘  “65  •  <£  *  & 


(10.6) 


l  B*(L 
(3  £  a  a 


2S  ) 
a 


where  B  Is  the  strength  of  the  magnetic  field  and  pg  is  the  Bohr  magneton. 
From  group  theory,  we  can  deduce  the  transitions  induced  by  equation  (10.6)  by 
using  first  table  25  of  Koster  et  al .  This  table  shows  that  Lz  and  Sz  are 
basis  functions  for  T1  ;  L+1  and  S+1  are  basis  functions  for  Tg ;  and  L_1  and 
S_ 1  are  basis  functions  for  Tn.  Then  from  table  26  of  Koster  et  al,  we  have, 
for  B|Z  (Lz  -  1*1 ), 


ri  x  r,  -  r. 


r,  x  r„  -  r. 


ri  *  r3 


-  r. 


ri  *  rH 
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and  for  B|X  or  B|Y  (L±  ♦  iy  r^)  , 


r3  *  F1  "  r3  ’ 


r4  "  ri 


r3  x  r2  "  r4  • 


r4  «  r2  “  r3  * 


r3  *  r3  "  r2  * 


r4  *  r3 


r3  *  r4 


r4  *  r4 


The  allowed  transitions  (In  absorption)  are  shown  in  figure  10.1a, 
where  the  labeling  is  according  to  the  orientation  of  the  light  polarizer. 
That  is,  the  polarization  is  determined  by  the  orientation  of  the  eleotrio 
vector;  for  u  polarization,  Ejz  (and  BjZ) ,  and  for  a  polarization,  EjjZ  (and 
B  |  Z ) . 


The  eleotrio  dipole  operator  is 


•  r 


•  I 


<10. 7) 


From  table  25  (Koster  et  al ,  1963),  we  find  that  C10(z)  is  a  basis  for  r2 
while  C1±1  (or  x  ±  iy)  has  the  same  representations  as  L^.  For  C10,  from 
table  26  (Koster  et  al),  we  have 


r  x  r  -  r 
->  *1  l2  ’ 


r2  *  r2  *  F1  ’ 


r2  x  f3  “  r4 


r2  x  r4 


and  the  allowed  transitions  (in  absorption)  are  shown  in  figure  10.1b.  In  the 
electric  dipole  case  we  have  to  assume  that  the  odd  terms  in  the  crystal  field 
given  in  equation  (10.1)  mix  either  the  p  or  f  configuration;  otherwise,  the 
electric  dipole  matrix  elements  vanish. 


As  a  second  example  we  consider  the  d^  configuration  in  the  crystal 
field  given  in  equation  (10.1).  The  Hund  ground  state  is  ^F  (S  «  3/2,  L  »  3) 
and  the  only  other  state  with  S  -  3/2  is  the  (all  the  states  of  pN,  dN,  and 
fN  are  given  on  pages  1  through  3  in  Nielson  and  Koster,  1  963,  and  pages  15- 
14  through  15-20  in  Polo,  961).  In  the  absence  of  spin-orbit  coupling,  the 


r.i  .MM  sa&iX-iMXiLlCVJCiiCii.  ifWXV 


state  is  the  only  state  that  couples  to  the  J*F  Hund  ground  state.  From 
table  30  of  Koster  et  al  we  have 

°3  *  ri  +  2r2  +  2r3.4  ’ 

(10.3) 


and  we  have  a  2x2  secular  equation  for  ^  and  r2,  and  a  3x3  secular  equation 
for  the  ^4  states.  Operating  on  the  states  | LM>  with  the  operator  Sj,  pro¬ 
duces  the  values  of  M  belonging  to  the  different  irreducible  representations. 
The  resulting  states  are 


|30>  , 

|10>  , 

for 

ri  ’ 

3-2>  , 

!  32  >  , 

for 

r2  ■ 

(10.9) 

3-1  >  . 

1 33>  , 

1 1-1  >  , 

for 

r4  • 

The  states  for  T3  are  obtained  by  changing  the  projections  (M  -*•  -M)  of  the 
states.  The  matrix  elements  of  the  crystal  field  for  the  states  given  in 
equation  (10.9)  are  presented  in  table  10.3.  7  The  results  given  in  table  10.3 
are  also  applicable  to  the  configuration  d  ,  provided  that  all  values  are 
multiplied  by  minus  one.  This  latter  result  is  obtained  by  the  rules  for  the 
states  conjugate  to  the  configuration,  which  has  4<l+2-N  electrons;  these 
rules  are  given  in  Nielson  and  Koster  (1963).  The  results  given  in  table  10.3 
were  taken  from  the  Polo  tables,  and  equation  (9.6)  was  used.  The  3-j  symbols 
were  found  in  Rotenberg  et  al  (1  969). 

Before  the  energy  levels  can  be  found  using  the  results  of  table 
10.3,  the  Coulomb  interaction  must  be  added  to  the  diagonal  elements.  From 
Nielson  and  Koster,  these  are 

w  f1**)  .  .  11  p(2)  _  72  _(  4) 


(10.10) 


The  energy  levels  of  d^  in  symmetry  is  given  in  table  10.4  for  a  represent¬ 
ative  set  of  Bkq  and  F^k'.  Included  in  the  table  is  the  effect  of  spin-orbit 
coupling,  which  we  consider  next. 

The  above  examples  have  been  restricted  to  the  orbital  angular 
momentum  only.  This  restricts  the  use  of  group  theory  to  the  single  group. 
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H,(4p) 


147  _(4) 

W  F 
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TABLE  10.3. 
OF  CRYSTAL 

Note:  For 
for  and 

,  MATRIX  ELEMENTS.  IN  Sj,  SYMMETRY, 
FIELD  FOR  *>F  AND  *>P  STATES  OF  nd3 

the  corresponding  matrix  elements 
states  of  nd?,  multiply  all 
entries  by  -1 .a 

TABLE  10.4 
Cr3  + 

-  IN 

C 

Free 

i0"  ID 

1.  ENERGY  LEVELS  FOR 

4F  AND  '‘P  LEVELS 

Sjj  SYMMETRY3 

-  0  C  -  275  cm-1 

Energy  Energy 

(cm"M  IR  (cm-’) 

L'M 

'  L  M 

IR* 

B20 

b40 

b44 

3  o 

3  0 

r1 

-4/35 

2/7 

0 

V 

r1 

0 

r5.6 

0 

3  0 

1  0 

r1 

12/35 

-4/21 

0 

r3,4 

352 

r7,8 

94 

1  0 

1  0 

r1 

2/5 

0 

0 

r3,4 

1  .574 

r7,8 

249 

3-2 

3-2 

r2 

0 

-1/3 

0 

^2 

2,203 

r7,8 

496 

3-2 

3  2 

r2 

0 

0 

/70/21 

r2 

3,672 

r5,6 

622 

3  2 

3  2 

r2 

0 

-1/3 

0 

r5,6 

788 

3-1 

3-1 

-3/35 

1/21 

0 

r7,8 

1,713 

3-1 

3  3 

r4 

0 

0 

A2/21 

r7, 8 

1,748 

2-1 

1-1 

^4 

2/6/35 

/ft/21 

0 

r5,6 

1 ,783 

3  3 

3  3 

^4 

1/7 

1/7 

0 

r5,6 

1  ,901 

3  3 

1-1 

r4 

0 

0 

2/7/21 

r5,6 

2,381 

1-1 

1-1 

^4 

-1/5 

0 

0 

r7,8 

2,460 

arhe  matrix  elements  of  f 

?  are  not  given. 

r7 ,8 

3,918 

as  \ 

they  are  equal 

to  those  for  f . , 

Before 

r5  6 

3.930 

the 

matrix 

is  diagonalized , 

the  Coulomb 

energies  from  Nielson  and  Foster  (1963)  for 

the 

qF  and 

*P  states  should 

be  added  to  the 

P'S  Jj 

16,294 

16,286 

dial 

Tonal  elements 

• 

5,  ** 

0,0 

3 

°IR  m  irreducible  representation 

• 

rl 

1 6 , 504 

r5,6 

16,373 

r7,8  ’M42 

Vrj  g  16,621 
r7f’8  16,7^1 

_ r5,6  16>771 

_«B20  «  296.9  c->B40-4597 
cm  .  B44  -  2544  cm  \  - 

74,202  cm";  F(4'  -  45,922  cm-2. 

The  F(2) ,  F'*  '  ,  and  C  (275  cm-2; 
vaiues  are  for  the  free  ion. 

X/?  “  irreducible  representation. 
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If  we  oonsider  the  case  of  the  d3  as  above  but  assume  the  spin-orbit  energy  to 
be  strong,  we  use  total  angular  momentum  states,  J,  with  J  -  L+*.  We  further 
restrict  the  disoussion  to  the  states  ^Fj  with  3/2  £  J  £  9/2.  Thus,  we  have 
the  levels  ^3/2*  ^5/2,  ^7/2*  and  ^9/2  to  oonsider.  From  table  30  (Koster 
et  al)  we  obtain 

D3/2  *  r5,6  +  r7,8  • 


°5/2  +  r5,6  +  2r7,8 


°7/2  *  2r5,6  +  2r7, 8 


(10.11) 


D9/2  "  3r5,6  *  2r7,8 


and  we  see  that  the  secular  equation  for  or  Tg  is  7x7,  and  for  or  Tg  the 
secular  equation  is  also  7«7.  We  see  that  the  inclusion  of  the  spin-orbit 
energy  greatly  increases  the  difficulty  of  the  problem.  The  wave  functions 
belonging  to  the  different  irreducible  representations  can  be  obtained  from 
table  25  of  Koster  et  al  and  are 

|J  1  +  4p>  ,  where  p  -  0,  ±1,  ...,  |g-  +  4p  j  Si  J,  for  , 

(10.12) 

|J  -  ^  +  4p>  ,  where  p-0,  ±1,  ...,  |-  ^  +  4p|  £  J,  for  . 

The  resulting  energy  levels  of  the  Cr^*  ion  with  the  Slater  parameters  and  the 
spin-orbit  constant  from  the  free  ion  are  given  in  table  10.4.  The  crystal- 
field  parameters  are  rough  estimates  for  the  crystal  field  for  Crd+  in  the  Ga 
site  (S^)  in  the  material  GdgS^GagOi  2. 

The  previous  two  examples  were  for  the  group  and  were  simple  to 
manipulate  in  that  the  group  was  cyclic  (all  the  operations  can  be  expressed 
in  terms  of  a  single  generator,  S^). 

We  consider  a  single  d  electron  in  a  crystal  field  of  D2  symmetry. 
The  crystal-field  interaction,  given  by  use  of  table  8.1,  is 


H3  “  B20C20  +  B22^C22  +  C2-2^  *  ImB32^C32  “  C3-2^ 

(10.13) 

+  B40C40  +  B42^C42  +  C4-2^  +  B44^C44  +  C4-4^  * 
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am  Km  am  am  am  *  nttn  MAM  hmmmauama 


jr**i,*  ■  *  w  a  wb  .  e  u  wiiki*  win. 


The  operations  of  D2  are  given  Jn  table  17  of  Koster  et  al  and  are  C2  -  C2(z), 
C2  “  c2(y)»  311(1  "  Cg(x)  .  The  result  of  these  operations  on  |  tm>  is  given 

in  problem  2  of  ofiapter  1.  These  are 

C2<y)|S,m>  -  (-1  )2,+m| n-m>  , 

C2(z)|im>  -  (-1  )m| Jlm>  ,  (10.14) 

C2(x)|lm>  -  (-1  >*'| H-m>  . 

Now  from  table  21  of  Koster  et  al ,  the  |2m>  states  deoompose  to 

°2  2ri  *  r2  +  r3  *  r4  ’  (10.15) 

and  the  T.  states  are 

<|n  -  j  20 >  ,  -  [  1 22>  +  1 2-2 >]  (10.16) 

/2 

which  are  essentially  the  crystal  potential  for  the  C?_  terms  in  equation 
(10.13).  Q 

For  the  other  representation  we  take  the  combination 

%  -  —  [  1 21  >  -  (2-1  >]  for  rp  , 

/2 

-  —  C|22>  -  |2-2>]  for  r,  , 

*u  “  —  C | 21 >  +  | 2-1 >]  for  r,  , 

/2 

which  can  be  verified  by  using  equation  (10.14)  and  the  character  table  for 
the  D2  group  in  Koster  et  al . 

10,2  Annotated  Bibliography  and  References 

Bullhausen,  C.  J.  (1962),  Introduction  to  Ligand  Field  Theory,  McGraw-Hill, 
New  York,  NY. 

Griffith,  J.  S.  (1961),  The  Theory  of  Transition-Metal  Ions,  Cambridge  Univer¬ 
sity  Press,  Cambridge,  U.  K. 

101 


(10.17) 

(10.18) 
(10.19) 


Konig,  E. ,  and  S.  Kremer  (1977),  Ligand  Field  Energy  Diagrams,  Plenum  Press, 
New  York,  NY.  This  book  contains  a  large  number  of  energy  levels  of 
3dN  ions  in  different  crystalline  environments.  Pages  18  through  22 
describe  the  different  notation  used  to  specify  the  constants  (crystal- 
field  parameters)  of  the  system. 

Koster,  0.  F. ,  J.  0.  Dimmook,  R.  G.  Wheeler,  and  H.  Statz  (1963),  Properties 
of  the  Thirty-Two  Point  Groups,  MIT  Press,  Cambridge,  MA. 

McClure,  D.  S.  (1959),  Electronic  Spectra  of  Molecules  and  Ions  in  Crystals, 
Part  II. — Spectra  of  Ions  in  Crystals,  Solid  State  Phys.  9,,  399. 

i 

Nielson,  C.  W.,  and  G.  F.  Koster  (1963),  Spectroscopic  Coefficients  for  the 
pn,  dn,  and  fn  Configurations,  MIT  Press,  Cambridge,  MA. 

Polo,  S.  R.  (1961,  June  1),  Studies  on  Crystal  Field  Theory,  Volume  I — Text, 
Volume  II — Tables,  RCA  Laboratories,  under  contract  to  Electronics  Re¬ 
search  Directorate,  Air  Force  Cambridge  Research  Laboratories,  Office  of 
Aerospace  Research,  contract  No.  AF  1 9( 60M ) -55^1 .  [Volume  II  gives  date 
as  June  1 ,  1 961 . ] 

Rotenberg,  M.,  R.  Bevins,  N.  Metropolis,  and  J.  K.  Wooten,  Jr.  (1969),  The  3~j 
and  6-j  Symbols,  MIT  Press,  Cambridge,  MA. 

Tinkham,  M,  (1964),  Group  Theory  and  Quantum  Mechanics,  McGraw-Hill,  New  York, 
NY. 


102 


11.  NUMERICAL  EXAMPLE:  4Fj  STATES  OF  Nd3+  (4f3) 

As  a  numerical  example  of  the  calculation  of  the  crystal-field  splitting, 
we  calculate  the  splitting  for  Nd3+  (4f3)  in  a  field  of  Sjj  point  symmetry.  We 
assume  that  the  levels  are  pure  4Fj.  We  assume  that  L,  S,  and  J  are  all  good 
quantum  numbers;  then  we  consider  matrix  elements  of  in  equation  (9.9)  with 
J'  -  J  and  L'  «  L.  Thus, 


<JM*LS|H.|  JMLS>  -  l  B*  <J(M)k(q)| J(M')>  <j|Uk|j>  <£(  0  )k(  0)  |  )l(0)  >  .  (11.1) 

J  Isn  K9 


The  values  of  the  reduced  matrix  elements  <j|uk|j>  <fc|Ck|Jl>  for  the  V  state 
of  Nd3+  are  as  in  table  11.1.  In  obtaining  these  values  we  have  used  Nielson 
and  Foster's  results  (1963)  in  equations  (9.9)  and  (9.10)  for  the  reduced 
matrix  elements,  (LSjUk|LS),  for  L  -  3  and  S  ■  3/2;  the  Racah  coefficients 
(6-j  symbols)  are  found  in  Rotenberg  et  al  (1969). 


The  calculation  of  the  energy  levels  is  made  somewhat  simpler  by  using 
wave  functions  that  transform  according  to  some  irreducible  representation  of 
the  Sjj  group  (Foster  et  al ,  1963).  The  irreducible  representations  of  the  Sjj 
group  are  all  one-dimensional ,  but,  since  the  ion  under  investigation  has  an 
odd  number  of  electrons,  the  energy  levels  will  be  at  least  doubly  degener¬ 
ate.  Thus,  of  the  four  irreducible  representations — Tg,  Tg,  Ty,  and  Tg — only 
two  need  be  chosen;  we  chose  Tg  and  Ty.  The  energy  corresponding  to  Tg  is 
degenerate  with  T5,  and  that  corresponding  to  Tg  is  degenerate  with  Ty.  The 
wave  functions  belonging  to  Ty  with  a  particular  J  value  are 


,  8k  +  1 

J  -  > 

2 


2  J  +  1 


9 
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TABLE  11.1  VALUES  OF  <j|Uk|j>  /2fc*1  <«,|Ck|lt>a 


J 

k-2 

k-M 

k-6 

3 

-  lu-'t 

J1/2 

0 

0 

2 

5 

11  1 

‘  1  1 1 /2 

1  f 

’  1  11/2 

n 

2 

2-5  1 

.2-5 -7J 

2  L 

„2*3*7J 

7 

1  r  5 11/2 

1  [■  1  11/2 

2*5  r  1  11/2 

2 

3  L3  *  7j 

2-3  L7-11J 

3  L3  •  1 1  *  1 3 J 

9  _1_  f  Hll/2  1  f  13  11/2  5  f  1  11/2 

2  2‘3  1.2 •  3 J  "  2/3  L2 •  1 1 J  3  * 3 *  1 1  •  1  3 J 


al Vote:  2*5  ”  10,  etc. 


and  those  for  Te;  are 


8k  -  3 


2 


> 


1  k  S 


2J  +  3 

8 


9 


where  k  is  an  integer,  and  the  number  of  k  values  occurring  for  a  given  J  is 
the  number  of  times  a  representation  will  occur.  The  number  of  Bkq  for  the 
calcium  site  (S^)  in  calcium  tungstate  is  five:  B20,  B^q,  Bg0,  B^,  and 
B61).  Of  these  parameters  only  Bgi,  is  complex.  The  matrix  elements  of  the 
crystal  field  given  in  the  above  equation  are  presented  explicitly  under 
particular  states  in  the  following  paragraphs. 


11.1 


This  level  of  the  free  ion  is  split  into  two  doublets  by  the  crys- 

talline  field.  The  wave  functions  corresponding  to  r„  and  rc  are  lx  ■x  > 
3  3  (  o  '  d  d 

and  If  "  f  >  respectively.  From  equation  (11.1)  we  have 


and 


2  1  lu  I  1  1  „  .  2_  B 

2  2  1  31  2  2  25  20 


3  _  3  ,  |  3  _  3  >  „  2_ 

2  2  |H3 '  2  2  >  25  B20 


(11.2) 


(11.3) 


fi  WkV  n  ,r.A'dx^Kd  f  mn  mim  a*  ^  n  ^  - - 


where  the  appropriate  values  of  the  reduced  matrix  elements  in  equations 
(11.2)  and  (11.3)  were  taken  from  table  11.1.  The  total  splitting  of  the 

U 

F3/2  state  is  then 


R 

25  20 


(11.4) 


11.2  > 


5/2 


Unlike  the  previous  case,  this  state  contains  two  I'-'s,  and  their 

C  C  C  o  ^ 

wave  functions  are  -  >  and  ~  >.  The  wave  function  for  the  state 


.  5  1 

is  |-  >.  The  energy  for  is 


<  5  1  j  ,51 
2  2  1  31  2  2 


1 1 
700 


(",B20  *  §  ■  E(r7  I)  •  ("-5) 


The  necessary  matrix  elements  for  the  energy  in  r,.  are 

25 


<  |  |  Ih3I  II  > 


11 

700 


t5B20  -  *>11 


(11.6) 


<  2  “  I  lH3^  "2  “  2  >  "  700  [”B20  "  11  B40^  "  b22  ’ 


(11.7) 


<1-1  |H31  II  >  -  —  B,.,.  -b 


6/nr 


44  1 2 


(11.8) 


The  two  energy  levels  corresponding  to  are 


,  /„  5\  b11  +  b22  +  ^b11  "  b22^  +  ^b12b?2^ 

V  5  2) - 2 - 


(r5  I) 


b11  +  b22  "  ^b11  "  b22^  +  4b12b12^ 


(11.9) 


2 


(11.10) 


11.3  F 


7/2 


fur  Ty  are 


This  state  contains  two  rc'3  and  two  r_’s.  The  matrix  elements 

o  f 


<  H  l«,|  y  y  > 


1  r  55  9  50  -i 

2  2  '“3 1  2  2'  9"9  *■  7  B20  +  T¥  B40  +  TT  B60-l  "  a11 

<  2  "  2  I H3 I  2  "  2  >  99  ^ 1 1 B20  +  2  B40  "  f3  B60^  *  a22  ’ 

<1-1  |h3|  11  >  .  ±  ^  B,,  -  If  /To  B, 


64 


»1 2  . 


The  two  energy  levels  are 
Mr7  l>  ■ 


7  a„  *  .22  ♦[(»„-  *  "a,  j.^]1 


E2(r7  {)  - 


7i  al1  +  a22  '  ^311  "  a22^  +  i<al2a?2^ 


1/2 


(11.11) 

(11.12) 

(11.13) 


(11.14) 


(11.15) 


The  matrix  elements  for  r 5  are* 


<75||75>,_L  fli  B  _  13  o  +  50  i 
2  2  lH3 '  2  2  99  7  B20  iff  B40  13  B60-*  b 


1 1 


(11.16) 


<  1  -1  \U  I  7  .  3  ,  1  r_  33  n  _  3  p  90  _  i  h 

<  2  2  lH3'  2  2  >  99  L  7  B20  Tff  B40  13  B60J  “  b22 


(11.17) 


<  |  -  I  |H3I  \  |  >  -  TO  [/iTTTT  BjlU  +  fl  /2i?2  BfiJ  -  b 


99 


44  13 


64  ■ 


12 


(11.18) 


The  corresponding  energies  are  given  by  substituting  the  above  values  of  b^ . 
into  equations  (11,9)  and  (11.10).  J 


*The  symbol  a^  will  be  used  for  the  matrix  elements  In  T7  and  b^  for  those  In  r5  to  avoid 
introducing  new  symbols  for  each  new  value  of  J. 
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11.4 


9/2 


The  number  of  r^'s  in  this  state  is  three,  with  two  iv's.  The 

«  J  _  1  J  A _  n  _ _  *' 


matrix  elements  of  the  crystal  field  for  are 


<  1 1  l"3l  I!  >  ■  ih  t6B20  -  T  %  *  fr  b6o)  ■  •„ 

<  I  i  l"3l  I  r  >  -  ffe  l-“B2o  -  ^  -  It  b601  - a 


22 


(11.19) 

(11.20) 


|H3I  jf-J>  -ps  [2B20  *  ^  B,0  -  B60]  -  a 


9  7 


22 


110 


33 


<21  lu  I  2  9  s  ,  7  r_  6  yp 
<?o  I H  o  I  TH?r  L  7  ®i 


2  2  '“31  2  2'  395" 

9  7  i„  ,  9  1 


150 


44  91  64 


12 


<2_2  lH3^  22>  “  395"  7  ^  B44  +  9l"  B64^  ”  a23  * 


(11.21) 

(11.22) 

(11.23) 


The  three  energies  are  given  by  the  solutions  of 


e3  '  (a1l  *  a22  *  a33)fi2  *  ta;ia22  *  31 1 a33  *  a22a33  -  a23a23  '  a12a*2)E 


+  a1ia23a23  +  a33a12a12  "  ana22a33  “  °  * 


The  matrix  elements  for  rc  are 

0 


<22  lH3 1  2  2  >  39?  ^B20  +  "T  B40  +  IT  B60-I  “  b11  *  (11.25) 

<  2  ~  2  2  '  I  >  jfe  ["3B20  T  B^0  +  9f  B60^  "  b22  ’  (H-26) 

<  2  '  2  lH3 1  2  2  >  39?  7  ^ 0  B44  “  9T  ^  B6iJ  “  b1 2  *  (H*27) 


The  energies  Ei and  £2^5  %)  are  given  by  equations  (11.9)  and  (11.10), 


respectively,  using  the  b^j  given  above. 
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11.5  Cal oul at  Ions 


Three  of  the  orystal-f ield  parameters  oan  be  obtained  quite  simply 
from  the  experimental  data.  These  are  E^q,  B4q,  and  BgQ.  If  we  express  the 
sums,  in  terms  of  the  Ei  where  the  E^  are  experimental  data,  then  we 
obtain 

SJ-1/2  “  ^  Ei^r7J)  * 
and 


-  -1  B 

25  20  ' 

(11.28) 

-  11  R  +  1  R 

1  75  20  142  40  * 

(11.29) 

2  R  *  8  R  +.  d 

53  B20  593  B40  T2&7  B60  ’ 

(11.30) 

iB _ Z_  n  -  J°  B 

9  °20  T5ff  B40  lT5?  “60  ’ 

(11.3D 

where  S4  is  for  the  J  <*  9/2 

level . 

These  equations  can 

be  inverted  to  give 

CD 

rv> 

o 

B 

.  25 
“2  S1 

P 

(11.32) 

B40  = 

-33S1  + 

42S2  , 

(11.33) 

B60  “ 

1001  - 
“W  si 

78  .  1287  - 

"  “5  S2  +  ~W  S3  * 

(11.34) 

The  other  crystal-field  parameters  are  slightly  more  involved.  f 

From  equations  (11.9)  and  ( 1 1 .1 0)  we  have 

b4J|  -  3/nr  [w2  -  n2]1/2  , 

where 

"t  --Ts,  “’d  "l  -  E,  (r5  I)  -  e2  (f5  !)  • 

To  determine  Bg4,  we  use  equations  (11.14)  and  (11.15)  to  give 

a,2a?2  -  |  [wf  - 


108 


I 

traniwun  wv  kv kv .  -  in «  . »  v  Km  m  uhum mj«  m »  n mjt  ns  HJin*  .  n  *_*  *  i  jji m* mji *  n  a/i 


(11.35) 


I 

« 

t. 

i 

(11.36)  ; 


k 

c 


where 


r  S1  - 1  s2  *  i  h  and  w2  ■  E1  (r7 I)  -  Ea(r7  i)  ■ 


A  similar  expression  can  be  obtained  using  equations  (11.16)  and 
(11.17)  in  equations  (11.9)  and  (11.10),  yielding 

b12b1*2  "  ¥  [W3  -<]  *  (1137) 

where 

h3  -T3!  -  2S2  *ls3  and  W3  '  E,M)  -  ^M)  • 


Substituting  equations  (11.13)  and  (11, 14)  into  the  left  side  of 
equations  (11.36)  and  (11.37),  we  obtain  two  equations  for  Bg^.  These  two 
equations  can  be  solved  simultaneously  for  both  real  and  imaginary  parts  of 
Bg4  to  give 


11  CL 

20  L 


99 


(" 2  *  w3 


-  n|)  -  6 


K  -  Ni)] 


COS  0  - 


13 


15*99  r,,2 
~ S  (W3 


N|)  -  8(W‘  -  N‘) 


2  >  7-99 


1/2 


(WJ 


N  l) 


20/70 


r6(w,  -  **) 


2n  1/2 


(11.38) 


(11.39) 


where  Bg^  =  Rgei9. 

All  the  crystal-field  parameters  can  be  determined  once  the  experi¬ 
mental  data  are  taken  on  the  Tg/j*  ^5/2*  311(1  ^F7/2  levels* 

As  tedious  as  the  above  procedures  may  have  seemed,  the  crystal- 
field  parameters  we  obtain  are  only  approximate,  since  we  have  ignored  L-S 
mixing  by  the  spin-orbit  coupling  in  the  free  ion  and  J  mixing  caused  by  the 
crystal  field.  Nevertheless,  the  crystal-field  parameters  obtained  by  the 
above  procedure  can  serve  as  very  good  starting  values  in  the  fitting  of  a 
more  sop'  isticated  calculation  to  experimental  crystal-field  levels. 


The  crystal-field  parameters  B^  obtained  by  the  above  procedure 
for  Nd^+  in  CaWOjj  are  given  as  follows,  along  with  crystal-field  parameters 
for  the  same  ion  but  with  full  diagonalization,  that  is,  L-S  mixing  and  J 
mixing  (Wortman  et  al,  1977). 
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B*q 

B20 

b40 

b44 

B60 

B64 

Re  Im 

Calculated 

above 

403 

-635 

711 

-219 

885  0 

Calculated  with 
full  diagon- 

509 

-866 

1042 

-509 

903  243 

alization 
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12.  CLASSICAL  POINT-CHARGE  MODEL 
1 2.1  Discussion 

In  the  simplest  model  of  the  crystal  field,  the  point-oharge  model 
introduced  by  Bethe  (1929),  the  lattice  is  replaoed  by  an  array  of  point 
charges  placed  at  the  nuclei  of  the  constituent  ions.  A  multipole  expansion 
is  made  of  the  point-charge  potential  energy  at  the  rare-earth  site.  Thus, 
if  Rg,mn(j)  is  the  vector  position  of  constituent  J  at  site  j  in  the  l,  m,  nth 
cell ,  we  have 


-e2Z 


L 


3  *»"  J  '  fl 


(12.1) 


where  ^j_mn  -  ta  +  mb  +  nc  +  pj ,  and  a,  b,  and  c  are  lattice  vectors.  The 
charge  at  site  j  is  eZ. ,  and  r  is  the  position  of  an  electron  on  the  ion  being 
discussed.  The  multipolar  expansion  of  equation  (12.1)  is 
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(12.2) 


The  multipolar  crystal-field  components  Akq  are 

y  y  zjckq[^tmn^  )  3 

:mn  j  [Rlmn^J  J  l1'1'*’1 


A.  -  -e 
kq 


Thus  the  point-charge  Hamiltonian  is 
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(12.3) 


(12.4) 


N 

where  we  have  summed  over  all  the  N  electrons  In  the  48,  configuration.  If 
all  the  lengths  are  measured  in  angstroms  (10~8  om),  then 
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,  r  Z.1clA.n(J)] 
inn  J 


(12.5) 


where  the  fine-structure  oonstant,  cu, 
and  the  units  of  Akq  are  onf'/A  .  If 
cm  • 


is  /fio  (thus,  a  /2ir  x  10®  -  116,1110), 
<rk>  is  in  angstroms,  then  Akq<rk>  is  in 


The  sum  in  equation  (12,5)  always  converges — even  for  the  lowest  k 
value  (k  -  0) — if  taken  in  the  order  indicated.  That  is,  the  sum  over  j  is 
performed  with  £,  m,  and  n  fixed.  The  unit  cell  is  neutral,  that  i3, 


» 


* 


l  Z  -  0  .  (12.6) 

J  J 

In  many  (but  not  all)  space  groups,  it  is  possible  to  choose  an  origin  for  the 
lattice  coordinates  such  that  the  dipole  moment  of  the  unit  cell  vanishes; 
that  is, 


l  PjZj  -  0  ,  (12.7) 

J 


where  is  the  position  of  the  jth  charge  in  the  unit  cell.  The  result  in 
equation  (12.7)  can  be  anticipated  by  observing  the  point  symmetry  of  the  ions 
in  a  specific  solid.  If  the  ions  occupy  C1  ,  C2,  Cg ,  C2v,  Cji,  C^y,  C^,  C^y, 
Cg,  or  Cgv  point  symmetry  (Schoenflies  notation),  then  it  is  impossible  to 
satisfy  equation  (12.7)  with  these  sites  as  the  origin  in  a  unit  cell.  I1"  it 
is  possible  to  satisfy  equation  (12.7),  then  the  sum  given  in  equation  (12.5) 
converges  very  rapidly.  This  can  be  shown  from  the  expansion 


Ckq(^"x) 
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<a(a)k(q)  |a+k(ot+q)>  xaC  (x) 

cl  01 
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(12.8) 


(Carlson  and  Rushbrooke,  1  950).  With  -  x  and  ^mn(0) 
”  ^Hmn^  +  Pj)»  f‘or  the  sum  in  equation  (12.5)  we  have 
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Now  if  equation  (12.7)  is  satisfied,  then 


(12.9) 


J  ZjpJC1>V  "  0 


(12.10) 


Thus,  the  sum  in  equation  (12.9)  is  for  a  >  1 ;  we  see  that  even  for  the  lowest 
terra,  k  -  0,  for  large  A,  m,  and  n  the  individual  terms  decrease  as 
1/r2  (o).  While  the  expansions  in  equations  (12.8)  and  (12.9)  are  good  for 
demonstrating  the  rate  of  convergence,  the  computation  of  Ak_  by  equation 
(12.5)  is  more  practical.  However,  in  equation  (12.5),  the  sum  over  j  should 
be  done  for  each  cell  first,  with  fixed  values  of  A,  in,  and  n.  In  programming 
language,  this  is  expressed  by  stating  that  the  j  loop  is  the  innermost  of  the 
nested  A,  m,  n,  and  j  loops.  In  some  lattices,  the  condition  in  equation 
(12.7)  may  place  some  of  the  point  charges  on  the  cell  faces.  In  these  oases 
it  Is  a  simple  matter  to  oalance  these  charges  by  an  adjustment  to  fractions 
of  equal  charges  on  opposite  faces. 


The  convention  we  use  for  our  lattice  sums  Is  that  given  in  the 
International  Tables  for  Crystallography  (1952);  table  12.1  is  reproduced  from 
volume  I  (the  other  two  volumes  give  data  strictly  for  x-ray  orystallog- 
raphers).  The  data  used  in  the  lattice  suras  are  generally  those  reported  in 
Acta  Crystallographies,  section  B  (now  predominately  section  C).  Care  should 
be  taken  to  make  certain  the  correct  setting  is  used. 


Typical  data  used  in  the  calculation  of  the  lattice  sums  are  given 
in  table  12.2  (LIYF^,  calcium  tungstate  space  group  88,  scheelite  structure). 
All  the  data  given  in  table  12.2  are  given  in  the  International  Tables,  except 
that  the  x,  y,  and  z  coordinates  are  determined  by  x-ray  diffraction.  The 
lattice  constants  a,  b,  and  c  are  also  determined  by  x-ray  diffraction  and,  as 
customary,  the  true  positions  of  the  ions  are  xa,  yb,  and  zc  (these  relations 
hold  for  all  the  ions  in  a  unit  cell).  The  polarizability  (from  Kittel,  1956, 
p  165)  of  each  ion  is  given  at  the  bottom  of  table  12.2.  For  this  particular 
solid,  only  the  fluorine  ions  can  have  dipole  moments  that  contribute  to  the 
crystal  field  (we  discuss  the  role  of  the  dipole  moments  in  sect.  I1*).  Not 
all  the  data  for  space  group  88  are  contained  in  table  12.2  because  the 
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TABLE  12.1.  CRYSTALLOGRAPHIC  AXIAL  AND  ANGULAR  RELATIONSHIPS 
AND  CHARACTERISTIC  SYMMETRY  OF  CRYSTAL  SYSTEMS 


Spaoe 

group 

System 

Axial  and  angular 
relationships 

X-ray 

for 

data  needed 
unit  cell 

1.  2 

Trlolinic 

a  *  b  *  o 
a  #  8  •*  Y  *  90° 

a,  b, 

c,  a,  8,  Y 

3  to  5 

Monoclinic 

First  setting 
a  *  b  *  o 
a  -  8  »  90°  *  Y 

a,  b, 

c,  Y 

Second  setting 
a  *  b  •  o 
a  -  Y  -  90®  *  8 

a,  b, 

c,  8 

16  to  7*1 

Orthorhombic 

a  #  b  «*  o 
a  -  8  -  Y  *  90® 

a,  b, 

c 

75  to  142 

Tetragonal 

a  -  b  *  c 
a  -  8  -  Y  -  90® 

a,  c 

143  to  167 

Trigonal 
(may  be  taken 
aa  subdivision 
hexagonal) 

(Rhombohedral  axes) 
a  *  b  •  c 

a  ■  8  “  Y  <  120®  «*  90® 
a  -  b  *  c 
a-8-90® 

Y  -  120® 

a,  a 

168  to  194 

Hexagonal 

a  -  b  -  c 
a  -  8  -  90® 

Y  -  120® 

a,  c 

195  to  230 

Cubic 

a  -  b  -  c 
o  -  8  -  Y  -  90® 

a 

Sourcei  International  Tables,  1952,  Vol.  I,  p  11,  table  2,3.1, 


TABLE  12.2.  CRYSTALLOGRAPHIC  DATA  FOR  L1YF4 
(SCHEELITE,  CaW04),  TETRAGONAL  SPACE 
GROUP  88  (FIRST  SETTING),  Z  -  4 


Ion 

Position 

Symmetry 

x 

y 

z 

Y 

4b 

s4 

0 

0 

1/2 

Li 

4a 

S4 

0 

0 

0 

F 

I6f 

C1 

0.2820 

0.1642 

0.0815 

Note:  a  *  5 

.1668  A ,  b 

-3,0” 

10.733  A, 

a  m 

90°  ,  i  -  90°  ,  Y  ■  90°  ,  <xy  -  0.55  A3,  aLi  -  0.05  A3, 
aF  m  1.04  A 3  (reduced  to  0.104  in  the  lattice  sum). 


Zy  -  +3,  ZLi  -  +1,  ZF  -  -I. 


m 


equivalent  positions  given  In  the  International  Tables  are  generated  Inside 
the  program.  The  centering  position  in  the  cell  can  be  taken  as  either  the  Y 
or  Li  site,  since  these  positions  have  Sh  symmetry,  and  their  lowest  crystal- 
field  component  is  ^20*  Equation  (12.10)  is  therefore  automatically 
satisfied.  The  resulting  lattice  sum  for  the  Y  site  in  Li YF ^  for  the  param¬ 
eters  in  table  12.2  is  given  in  table  12.3.  The  sum  covers  all  the  complete 
cells  in  a  sphere  of  30-A  radius  and  should  be  an  accurate  result.  Also 
included  in  table  12,3  are  the  results  for  the  dipole  contributions  due  to  the 
presence  of  dipoles  at  the  fluorine  sites  (see  sect.  14). 

Aj  a  seoond  example,  we  choose  a  very  low  symmetry  solid,  YCI3, 
which  is  character! zed  by  monoolinio  space  group  12.  As  can  be  seen  in  table 
12.4,  all  the  ions  are  in  very  low  point-symmetry  positions,  and  each  position 
can  have  a  dipole  moment  (another  way  of  saying  this  is  that  each  position  has 
a  onefold  field,  A1m).  We  then  have  to  consult  the  International  Tables  for  a 
higher  symmetry  position  in  order  to  satisfy  equation  (12.7),  which  in  this 
case  i 3  the  site  4e  with  Ci  symmetry.  The  Ci  point  group  ha3  only  the  inver¬ 
sion  operation,  and  all  the  odd-n  A__  vanish  in  this  symmetry.  Thus  if  the 
position  4e  is  used,  equation  (12.7)  will  automatically  be  satisfied.  The 
lattice  sum  for  YClo  was  also  run  over  a  lattice  30  *  30  x  30,  and  only  the 
even-n  A^  are  given  in  table  12.5.  The  dipole  contributions  were  also  calcu¬ 
lated;  these  calculations  were  more  complicated  in  this  solid  because  of  the 
three  types  of  sites  (Y,  Clf,  Clg);  all  have  a  dipole  moment.  For  many  of  the 
km,  the  dipole  contributions  are  much  larger  than  the  monopole  terms.  This 
frequently  happens  when  the  handbook  values  for  the  dipole  polarizabilities 
are  used.  We  have  had  more  believable  results  when  we  reduce  the  polarizabil¬ 
ity  to  one  tenth  of  the  handbook  value. 

The  lattice  sums  given  in  tables  12.3  and  12.5  are  incomplete  in 
that  the  results  are  not  in  a  usable  form  for  many  of  our  computer  programs. 
Before  we  can  use  these  results,  the  A^  should  be  rotated  using  the  standard 


TABLE  12. 

3.  LATTICE 

SUMS  FOR 

Y  SITE  AT  (0, 

0,  1/2) 

FOR 

liyf4  with  zy  -  3, 

ZLi  ■  -1* 

ZF  -  -1,  ap  ■ 

•  0.104  a 3 

Lattice 

Monopole  Anm 

Dipole  Anm 

Monopole 

and  dipole 

sum 

Real 

Imaginary 

Real 

Imaginary 

Real 

Imaginary 

A20 

1074 

0 

340 

0 

1414 

0 

a32 

373 

859 

-358 

74.0 

15 

933 

-1957 

0 

-98.1 

0 

-2055 

0 

A44 

-2469 

-2362 

-3.83 

-80.3 

-2473 

-2442 

A52 

1050 

-2456 

1.28 

-74.7 

1051 

-2531 

A60 

-17.2 

0 

7.96 

0 

-9.24 

0 

a6U 

-615 

-420 

-29.03 

-9.37 

-644 

-429 

A72 

-15.7 

0.90 

1.55 

-9.94 

-14.2 

-9.04 

A76 

250 

-63.8 

17.8 

7.96 

268 

-55.9 
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TABLE  12.4.  CRYSTALLOGRAPHIC  DATA 
FOR  YCI3,  MCNOCLINIC  SPACE  GROUP  12 
(C2/m)  (SECOND  SETTING),  Z  -  4 


Ion 

Position 

Symmetry 

X 

y 

z 

Y 

4g 

Cp 

0 

0.166 

0 

O  O 
H*  l-» 
\>  -* 

41 

ca 

0.211 

0 

0.247 

8J 

c1 

ci 

0.229 

0.179 

0.760 

“  ”■ 

4e 

1/4 

1/4 

0 

Notet  a  ■  6.92  A,  b  »  22.94  A,  c  ■  6.44  A, 
a  ■  90°,  B  “  11.0°,  Y  •  90 °t  Charges t  qy  m  3 , 


Polarizability  1  ay  -  0. 55  k3  ,  acl  -  3.66  k3. 


TABLE  12.5.  LATTICE  SUMS  FOR  Y  SITE  AT 
(0,  0.166,  0)  FOR  YClo,  EVEN-n  A 
ONLY,  ALL  A  REAL 


Lattice  sum 

Monopole 

Dipole 

Total 

a20 

1738 

3227 

4965 

a2i 

-913 

2916 

2003 

A22 

245 

2574 

2819 

A40 

-73.9 

246 

172 

A41 

85.8 

-398 

-312 

A42 

-41  .3 

47.7 

6.4 

A^ 

10.4 

-791 

-781 

a44 

-3.64 

516 

512 

a60 

-0.06 

-80.2 

-80.3 

a61 

-3.76 

21  .3 

17.5 

A62 

3.35 

-27.4 

-24.0 

a63 

-0.58 

60.5 

59.9 

a64 

2.73 

39.2 

41 .9 

A65 

5.07 

13.7 

18.8 

A66 

8.14 

-65.9 

-57.6 

Euler  angle-rotation  matrix,  so  that  the  lattice  sums,  A^,  rotated  from  A^ 
by  the  angles  a,  6 ,  and  Y,  are 


A' 

nm 


m 


D"  (a.Sj)!  . 

rn'm  nm' 


Explicit  forms  for  the 


D"tm(a,6,Y) 

rn'm 


are  given  in  Rose  (1957, 


oh. 
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13.  POINT-CHARGE  MODEL  DEVELOPED  AT  HDL 


In  the  classical  point-charge  model,  the  crystal-field  parameters,  B^, 
for  the  crystal- field  interaction  of  the  form 


H,  -  l  B*  l  C  (i)  (13.1  ) 

3  nm  f  nm 

nm  i 

were  calculated  as 

Bn«,  •  «•">*, m  ■ 

where  the  <rn>  are  the  expectation  values  of  r11  of  the  rare-earth  ion,  and  the 
Anm  are  the  multiP°le  components  of  the  energy  at  the  site  occupied  by  the 
rare-earth  ion.  In  the  earlier  models,  the  radial  integrals  used  in  the 
evaluation  of  rn  were  taken  from  Hartree-Fock  calculations  (Freeman  and 
Watson,  1962),  and  the  Am  were  calculated  using  the  point  charges  at  the 
valence  values  for  the  constituent  ions.  These  calculations  generally  gave 
the  twofold  field  10  times  too  large,  the  fourfold  fields  approximately  in 
good  agreement,  and  the  sixfold  fields  10  times  too  small. 

1  3*1  Screening  and  Wave  Function  Spread 

Several  errors  in  the  classical  theory  were  immediately  obvious. 
If  the  radial  wave  functions  (Hartree-Fock )  for  the  free  ion  were  correct, 
then  these  wave  functions  should  give  the  correct  values  for  the  Slater  inte¬ 
grals  F2,  F1*,  and  F^.  They  did  not  for  Pr^+.  A  simple  procedure  was  then 
applied.  The  radial  wave  functions  were  assumed  to  be  of  the  form 

♦(r)  -  CRHF(ir)  (13.3) 

where  x  is  a  parameter,  C  is  a  normal ization  factor,  and  RHF(r)  are  the  Har“ 
tree-Fock  radial  wave  functions.  With  the  radial  function  given  by  equation 
(13*3),  the  Slater  integrals  become 

Fk  -  xFkF  ,  (13.4) 
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and  it  was  found  that  a  value  of  t  of  approximately  0.75  was  needed  to  fit  the 
Fk  that  a^e  found  by  fitting  the  experimental  spectra  of  Pr’+.  Thus,  the 
Hartree-Fook  radial  wave  functions  had  their  maxima  too  near  the  origin  and 
needed  to  be  spread  out  even  in  the  free  ion,  and  perhaps  more  spreading  was 
neoessary  for  an  ion  in  a  solid. 

From  radial  wave  functions  given  in  equation  (13.3),  it  is  not 
difficult  to  snow  tnat 


<f(r)>  -  <<|» | f (r)|<J>>/C<t>|4>> 


*  ^f  (r/ t ) > 


(13.5) 


so  that  any  quantity  that  has  been  calculated  using  Hartree-Fock  functions  is 
immediately  obtained,  particularly 

<rk>  -  <rk>HF/xk  .  (13.6) 

A  second  error  of  the  classical  method  was  the  omission  of  the 
Sternheimer  shielding  factors  (Sternheimer ,  1951,  1966;  Sternheimer  et  al, 
1968).  In  1952  Sternheimer  showed  that,  in  a  multipolar  expansion  of  the 
energy  of  a  point  charge  embedded  in  a  solid,  the  r11  should  be  replaced  by 
rn(l  -  on) ,  where  the  on  are  the  shielding  factors.  He  further  showed  that 
these  factors  were  independent  of  azimuthal  angle;  that  is,  if  the  angular 
variation  in  the  multipolar  expansion  was  given  by  Ynm,  the  on  were  independ¬ 
ent  of  m.  The  values  of  crn  have  been  calculated  for  Pr3  +  and  TnH+  and  are 

u2  ■  0.666  ,  on  -  0.09  ,  og  -  0.04  for  Pr^+  , 

(13.7) 

o2  ■  0.545  ,  04  -  0.09  ,  05  -  0.04  for  Tm^*  , 

where  the  replacement  is 

rn  -  rn(l  -  on)  .  (13.8) 

More  recent  calculations  of  the  shielding  factors  have  been  done  (Sengupta  and 
Artman,  1970,  and  perhaps  others)  which  we  shall  need  if  further  refinements 
of  the  theory  are  undertaken. 

13.2  Effective  Charge  and  Position 

The  crystal-field  components,  km,  are  a  function  of  the  position 
of  the  ions  in  a  solid;  in  solids  such  as  CaWOjj  the  (WO4)”2  complex  is  known 
to  be  covalent.  That  is,  the  charges  on  the  tungsten  and  the  oxygen  ions  are 
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not  necessarily  at  their  valence  values.  If  we  let  the  charge  on  the  tungsten 
ion  be  (fy,  then  we  require  that 


%  *  ^Q0  “  "2  (13.9) 

with  the  charge  on  the  o/^gen  being  qQ.  The  result  given  in  equation  (13.9) 
then  assumes  that  the  Ca^  site  is  purely  ionic  with  charge  2.  This  assump¬ 
tion  is  consistent  with  many  of  the  experimental  results  on  compounds  such  as 
Ca»rUh  or  Vso,}.  We  introduce  a  second  parameter,  the  effective  position  of  the 
oxygen  ion  relative  to  the  tyj-igsten  site  that  would  reproduce  the  effective 
dipole  moment  seen  from  the  Ca2+  site.  This  parameter,  n,  is  introduced  by 

Ro-wCeffective)  -  nRo-vf(measur’e<;0  ,  (13.10) 

where  Rq_w  is  the  distance  from  the  oxygen  nucleus  to  the  tungsten  nucleus. 
Thus  there  are  two  parameters  in  the  A™:  qQ,  the  effective  charge  (qw  is 
eliminated  by  eq  (13.9)),  and  n,  the  effective  distance  of  the  oxygen  site 
from  the  tungsten  site.  The  calculated  crystal-field  parameters  Bnm  then  are 


®nm(T’^0,r*}  ’  ~  °n^nm^cb,ri^Tn  *  (13.11) 


with  the  three  parameters  t,  qQ,  and  n. 

The  experimental  data  that  were  taken  at  HDL  on  the  rare-earth  ions 
in  CaWO^  were  analyzed  using  the  effective  spin-orbit  Hamiltonian  (Karayianis, 
1970),  and  a  set  of  phenomenological  Bnm  was  obtained.  These,  given  in  table 
13.1,  are  the  that  the  theory  has  to  fit. 

The  fitting  was  done  by  minimizing  the  square  quantity  given  by 


Q 


2  CBniJT;q0’n) 

nm 


B  ]' 
nmJ 


(13.12) 


where  Bnm(T;q0,n)  is  given  by  equation 
for  each  ion.  The  minimization  was  done 
ion.  Since  the  qQ  and  n  are 
assumed  to  be  ion  independent 
and  t  is  assumed  to  be  host  in¬ 
dependent,  the  average  qQ  and  n 
were  taken  and  fixed.  The  proc¬ 
ess  was  then  repeated  with  mini¬ 
mization  with  respect  to  t 
only.  This  process  yielded  the 
following: 


% 


-1.09 


n  -  0.977 


(13.13) 


(13.11),  and 
with  respect 


Bnm 
to  T 


is  from  table  13.1, 
,  qQ,  and  n  for  each 
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IONS  IN  CaWO, 


(all  Tn 


cm 


-1 


Ion 

B20 

b40 

b44 

o 

\o 

CD 

B64 

Re 

Im 

Nd 

549 

-942 

1005 

-17 

947 

1 

Tb 

468 

-825 

872 

-290 

595 

160 

Dy 

428 

-825 

978 

-7 

H48 

2.5 

Ho 

436 

-664 

779 

-33 

558 

196 

Er 

404 

-685 

728 

12 

452 

164 

Tm 

417 

-688 

754 

17 

506 

359 
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The  t  values  were  well  fitted  by 

t  -  0.767  -  0.00896N  ,  (13*1^) 

where  N  is  the  number  of  f  electrons  In  the  configuration  4f^.  The  values  of 
on  used  in  the  above  were  not  varied  in  the  minimizing  process  but  were  inter¬ 
polated  from  the  values  given  in  equation  (13.7)s  that  is, 

o2  -  0.6902  -  0. 01  21  N  , 

o4  -  0.09  (all  N)  ,  (13.15) 

o6  -  -0.04  (all  N)  . 

The  predicted  values  of  the  Bnm ^ T * ^ ’ n ^  for  the  entire  rare-earth 
series  are  given  in  table  13.2.  The  results  given  in  table  13.2  when  oompared 
to  table  13.1  show  that  the  difference  between  the  derived  0,^(1;  qo,n)  and  the 
phenomenological  Bnm  is  greater  for  the  low-N  ions  in  the  4fN  configuration. 
This  may  be  a  defect  in  the  theory,  but  not  enough  data  on  the  low-N  ions  are 
available  for  analysis.  One  of  the  significant  results  of  the  analysis  was 
that  It  led  to  the  reanalysis  of  the  3peotrum  of  Tb:CaW04  with  a  different 
interpretation  of  the  experimental  data  (Leavitt  et  al,  1974). 

More  recent  work  on  CaW04  (Morrison  et  al,  1977)  obtains  the  fol¬ 
lowing  values: 

a2  -  0.6846  -  0.00854N, 
a4  -  0.02356  +  0.001  82N, 
t  -  0.75( 1.0387  “  0.01 29N) , 
o6  “  -0.04238  +  0.00014N, 
qQ  -  -1.150,  and 
n  ■  0.962. 

The  on  values  are  Interpolated  from  the  calculations  of  Erdos  and  Kang  (1972) 
for  Pr^+  ar.d  Tm  .  The  factors  in  equation  (13.11)  were  combined  so  that 


Pn  -  <rn>(l 


(13.16) 


and  the  pn  along  with  the  t  are  given  in  table  13.3.  Thus  we  have 


Bnm(T;qO,ri)  “  pnAnm^q0,n^ 


(13.17) 
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At  present  we  use  the  results  given 
crystal-field  parameters  given  by  equation  (13.17) 
starting  values  to  best  fit  experimental  data.  We 
n  •  1  in  the  prooess  (q  here  is  the  effeotive 
neoessar’ily  oxygen).  After  obtaining  the  best- 
oaloulation  of  Anm(q)  and  vary  q  to  obtain  the 
quantity 


Q  -  l  [  B 


nm 


nm 


in  table  13.3  to  oaloulate 
and  use  these  parameters  as 
generally  use  A^q^)  with 
charge  on  the  ligands,  not 
fit  we  return  to  the 

best  fit  by  minimizing  the 


03.18) 


Following  this,  we  obtain  the  An_(q)  for  odd  n  and  use  them  to  oaloulate  the 
intensities  using  the  Judd-Ofelt  theory. 


TABLE  13-2.  DERIVED  CRYSTAL-FIELD  PARAMETERS, 
B  (T;q0.n)  FOR  MfN  CONFIGURATION  OF  TRIPLY 
IONIZED  RARE-EARTH  IONS  (all  in  onf1 ) 


N 

Ion 

B20 

bmo 

bmm 

b60 

B6M 

Re 

Im 

1 

Ce 

MM1 

-1M29 

1M62 

16 

1251 

52 

2 

Pr 

M2M 

-122M 

1253 

13 

996 

M2 

3 

Nd 

J408 

-1059 

1083 

11 

805 

3M 

M 

Pm 

Mil 

-1017 

10M1 

10 

76M 

32 

5 

Sm 

M08 

-938 

960 

9 

676 

28 

6 

Eu 

M08 

-887 

908 

8 

626 

26 

7 

Gd 

M06 

-82M 

8M  3 

7 

559 

23 

8 

Tb 

M2M 

-856 

876 

8 

591 

25 

9 

Dy 

M28 

-831 

851 

7 

563 

2M 

10 

Ho 

Ml  7 

-756 

77  M 

6 

M88 

20 

11 

Er 

Ml  5 

-707 

72M 

6 

M39 

18 

12 

Tm 

M  35 

-729 

7M6 

6 

M5M 

19 

13 

Yb 

M3M 

-701 

717 

6 

M29 

18 

TABLE  13. 

3.  VALUES  FOR  t,  <rn>HP,  on,  AND  p 
OF  TRIPLY  IONIZED  RARE-EARTH 

FOR  MfN 
I0NSa 

CONFIGURATION 

Ion 

N 

T 

<r2>HF 

<r4>HF 

<r,6>HF 

°2 

°9 

°5 

02 

P9 

p6 

Ce 

1 

0.7693 

0.3360 

0.2709 

0.9659 

0.6757 

0.0259 

-0.0922 

0.1891 

0.7536 

2.3917 

Pr 

2 

0.7597 

0. 30*11 

0.2213 

0.3959 

0.6667 

0.0272 

-0.0921 

0.1756 

0.6969 

1.8759 

Nd 

3 

0.7500 

0.2803 

0.1882 

0.2715 

0.6578 

0.0290 

-0.0920 

0.1706 

0.5776 

1 .5897 

Pm 

14 

0.7903 

0.2621 

0.1655 

0.2297 

0.6988 

0.0308 

-0.0918 

0.1679 

0.5339 

1.9213 

Sm 

5 

0.7306 

0.2972 

0.1988 

0.1 929 

0.6298 

0.0327 

-0.0917 

0.1668 

0.5099 

1 .3210 

Eu 

6 

0.7210 

0.2397 

0.1353 

0.1686 

0.6309 

0.03M5 

-0.0915 

0.1666 

0.9836 

1.2503 

Gd 

7 

0.7113 

0.2232 

0.1737 

0.1977 

0.6220 

0.0363 

-0.0919 

0.1668 

0.9656 

1  .1873 

Tb 

8 

0.7016 

0.2129 

0.1131 

0.1287 

0.6130 

0.0381 

-0.0913 

0.1673 

0.9990 

1.1232 

Dy 

9 

0.6919 

0.2033 

0.1037 

0.1119 

0.6091 

0.0399 

-0.0911 

0,1681 

0.9391 

1  .0619 

Ho 

10 

0.6823 

0.1995 

0.0995 

0.0981 

0.5951 

0.0918 

-0.0910 

0.1692 

0.9217 

1 .0119 

Er 

1 1 

0.6726 

0,1865 

O.O883 

0.0879 

0.5861 

0.0936 

-0.0908 

0,1706 

0.9126 

0.9826 

Tm 

12 

0.6629 

0.1790 

0.0820 

0.0787 

0.5772 

0.0959 

-0.0907 

0.1722 

0.9053 

0.9699 

Yb 

13 

0.6532 

0.1717 

0.0753 

0.0681 

0.5b83 

0.0972 

-0.0906 

0.1737 

0.3938 

0.9120 

a 

The  units  of  <rn>HF  end 

P„  ere  An 

• 

At  present  we  have  not  included  the  dipole  contribution  to  the 
Anm(q)  but  intend  to  do  so  as  soon  as  possible.  The  old  n  in  the  three- 
parameter  theory  will  be  replaoed  by  a,  the  polarizability  of  the  constituent 
ions  in  low-symmetry  sites.  We  believe  that  this  procedure  (including  'inding 
new  pn  values)  will  give  muoh  better  results  than  obtained  by  the  older 
theory.  In  our  projected  reanalysis  we  will  have  the  good  phenomenological 
Bnm  for  R^jLaFg,  R3  :LaClg,  and  R3  +  iLiXF\  (these  are  reported  by  Morrison  and 
Leavitt,  1982),  and  will D soon  have  R3  s 2°3 *  in  ad(JitiGn  t0  the  B  for 
R^  +  :CaW0i1  used  in  the  older  theory.  These  data  should  be  sufficient  to  form  a 
very  stringent  test  of  a  newer  three- parameter  theory. 

For  the  ndN  ions  (Xq  for  q  -  +2,  +3,  and  +4,)  we  have  taken  the 
values  of  F^k'  obtained  by  fitting  the  free- ion  data  (Uylings  et  al,  1984) 
along  with  the  Hartree-Fock  values  of  F'k^  (Fraga  et  al,  1976)  and  obtained 
values  of  t  using  equation  (13.4).  Using  these  values  of  t  ,  we  obtained  the 
estimated  values  of  <rk>  from  equation  (13*6).  The  resulting  <rk>  are  given 
in  table  13.4. 

The  results  given  in  table  13.4  along  with  reported  values  of  Ak(}  (Morri¬ 
son  and  Schmalbaoh,  1985)  can  be  used  in  equation  (13.2)  to  obtain Morystal- 
field  parameters,  Bk  ,  whioh  can  be  used  as  starting  values  ln+ fitting  the 
experimental  data.  AM  similar  process  can  be  performed  for  the  X+q  (q  -  2,  3, 
and  4)  for  the  4dN  and  5dN  series.  However,  because  of  the  lack  of  free-ion 
parameters  F^k'  for  these  ions,  we  shall  have  to  interpolate  from  the  3dN 
series . 


TABLE  13.4.  ESTIMATED  VALUES  OF  <rk>  (Ak)  DIVALENT,  TRIVALENT,  AND 
QUADRIVALENT  IONS  WITH  3dN  ELECTRONIC  CONFIGURATION 


ndN 

X2  + 

<r2>a 

<r4>a 

X3+ 

<r2>b 

<rV 

X4  + 

<r2>c 

<rV 

3d2 

Sc 

1  .372 

4.053 

Ti 

0.7558 

1  .281 

V 

0.6217 

1  .298 

3d3 

3d3 

Ti 

1  .073 

2.505 

V 

0.6689 

0.9145 

Cr 

0.5172 

0.8177 

V 

0.8822 

1  .718 

Cr 

0.5776 

0.6911 

Mn 

0.4911 

0.7761 

3d 

Cr 

0.7423 

1  .234 

Mn 

0.5054 

0.5363 

Fe 

0.3958 

0.4955 

3d6 

Mn 

0.6293 

0.8973 

Fe 

0.4436 

0.4177 

Co 

0.3648 

0.4081 

3d7 

Fe 

0.5576 

0.7236 

Co 

0.4020 

0.3506 

Ni 

0.3304 

0.3282 

3d 

8d? 

Co 

0.4917 

0.5738 

Ni 

0.3627 

0.2903 

Cu 

0.2982 

0.2600 

Ni 

0.4353 

0.4577 

Cu 

0.3280 

0.2413 

Zn 

0.2708 

0.1995 

3d9 

Cu 

0.3871 

0.3678 

Zn 

0.2977 

0.2016 

Ga 

0.2291 

0.1296 

a<rJc>Wf./Tk  calculated  using  t  *  0.76878  +  0.11128N.  ! 

b<r*>Wf./T?  calculated  using  t  -  0.311184  +  0.0073953N. 

c<rk>WF./T  calculated  using  t  »  0.833540  +  0 .0056609N .  j 
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14.  CRYSTAL-FIELD  EFFECTS  NOT  YET  FULLY  INCORPORATED 


14.1  Self-Conslafcent  Point  Dipole  and  Point  Multipole 

In  seotlon  12  we  disoussed  the  point-charge  contribution  to  the 
multipolar  field  oomponenta  A^.  It  waa  early  reoognlzed  by  Hutohlnga  and  Ray 
(1963)  that  the  multipolar  oomponenta  of  the  constituent  Iona  oontribute  to 
the  A^,,  at  the  site  oooupied  by  the  unfilled  ahell  nd  .  For  a  point  oharge 
eZj  located  at  Ri  from  the  origin  ion  site,  we  have  the  eleotrio  potential 


♦(?) 


(14.1) 


The  potential  energy  of  one  of  the 


electrons  at  r 


is 


u(r)  -  -e$(r ) 


-  -e^Z, 


nm 


nm 


(14.2) 


where  we  have  expanded  the  denominator  of  equation  (14.1)  in  the  spherical 
tensors  discussed  in  chapter  1.  If  we  write  equation  (14.2)  as 


then 


U(r) 


A*  rnC  (r) 
nm  nm 


zl^nm^l ) 
Rn+1 


(14.3) 


(14.4) 
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where  the  sum  on  i  oovers  ail  the  iona  of  oharge  eZ1  in  the  solid.  This 
result  we  derived  In  seotion  11,  expressed  in  slightly  different  form.  It 
seems  natural  to  extend  equation  (1*4.3)  to  the  form 


u(r)  -  l  A<>>*rnc  fr)  ,  (14.5) 

n,m,k  ™  "" 

and  relate  the  to  the  various  k-pole  moments  of  ligands  at  ft..  To  relate 
the  to  the  multipole  moments,  eQkq,  we  need  first  to  express1the  eleotrlo 
potential  at  the  rare-earth  eleotron  due  to  the  multipole  moment  eQkg(i) 
at  ft^ 

The  eleotrlo  potential  due  to  a  multi pole  distribution  at  ft^  Is 

given  by 


♦(?) 


I 

kq 

am 


(”1  )kQkq(i)(2k2Sn) 


1/2 


<n(m)k(q) |n+k(m+q)> 


cft+k,m+q^l) 

Rfk  +  1 


r°C  (r) 
nnr  ' 


(14.6) 


where 


c2k+2n\  (2k+2n)l 

1  2n  >  “  T2nTlT2kTT 


(the  details  of  the  derivation  of  this  result  are  given  in  sect.  15).  Thus, 
since  u(r)  -  -e<j>(r),  we  find,  using  equation  (14.5),  that 


A(k)  - 

run 


-•* 


l  ("1)l<Qkq(i)(2n2nk)1/2  <n^m )k Cq ) | n+k (m+q ) >  - 

i 


q.i 


(14.7) 


If  we  let  k  -  0  in  equation  (14.7),  we  obtain 


-.2  ZQJ0(D 


R>J+1 


(14.8) 
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(if  k  -  0,  q  -  0),  which  is  identical  to  the  result  given  in  equation  (14.4) 
if  we  identify  Qg0(i)  with  ZA  (the  number  of  oharges)  givon  there.  The  result 


for  k  -  1  is 


Anra)  "  ®2  ^  qiq(1H2^2)1/2  <n(m)1  (q) |n+1  (m+q)>  • 

q,i  H  R1 


04.9) 


Sinoe 


<n(m) 1 (q) |n+1 (m+q)>  -  (-1)1  q  (§?t£)1/2  <1 (-q)n+1 (m+q) |n(m) > 


and 


a?q  ■  • 


we  oan  use  these  results  in  equation  (14.9),  to  obtain 


Anm^  "  “®2  E  /(n+1 ) TZn+Tf  Q1q(i)  <1  (q)n+1  (m-q) |n(ra)> 

Q » 1  i 


(*i) 


(14.10) 


which  is  identical  to  the  result  published  by  Morrison  (1976),  if  we  identify 
eQiq(i)  -  PiqO)  (p1q  is  the  dipole  moment  component).  Thus,  if  we  knew  the 
Qk  (i),  we  could  easily  calculate  the  A^*'  by  using  equation  (14.7).  Unfor¬ 
tunately,  the  real  difficulty  is  determining  the  Qkq(i).  In  what  follows  we 
restriot  our  discussion  to  the  dipole  oase,  k  -  1  ,  and  let  eQ1q  -  pq  and 
express  the  results  in  Cartesian  veotors. 


At  sites  of  low  symmetry,  an  electric  field  can  exist  whose  value 
is  determined  by  the  various  point  oharges  of  the  solid.  The  eleotrio  field 
due  to  the  point  oharges  of  the  solid  at  a  site  of  low  symmetry  is  given  by 


« -  - 1  ^ 

3  1  R?j 


(1 4.11 ) 


and  the  field  generated  by  the  point  dipoles  is 


if  -  I 

J  i 


3^ij(^ij*Pj)  Pj 


(14.12) 
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The  dipole  moment  at  site  j  is  then  given  by 

h  '  “Sj  '  “[Ej  *  f?]  •  t1"-’3) 

where  a  is  the  polarizability  of  the  ion  at  site  J.  (If  more  than  one  speoies 
is  considered,  the  polarizability  of  each  type  must  be  used.)  The  sum  in 
equation  (14.11)  presents  no  problem  and  oan  be  done  quite  simply.  To  perform 
the  sum  in  equation  (14.12),  it  is  convenient  to  assume  a  fixed  coordinate 
system  in  the  unit  oell  ar^d  an  associated  reference  point  (say  position  1); 
then  each  dipole  jnoment,  pj ,  oan  be  related  to  the  dipole  located  at  the 
reference  moment,  p1 ,  by  the  symmetry  operation  of  the  orystal.  Similarly, 
the  field  at  each  point,  £.,  can  be  related  to  E..  Having  done  this,  we  oan 
write  J 

-  0(1 )-pt  ,  (14.14) 

and  from  equation  (14.13) 

p,  -  a[i£®  +  0(1  )«p1j  .  (14.15) 

The  result  given  in  equation  (14.15)  oan  then  be  solved  for  the  dipole  moment 
to  give 

p,  »  aB(  1  )-'£°  ,  (14.16) 


where 

B  -  (1  -  oG) 


The  result  obtained  in  equation  (14.16)  is  rather  interesting;  if 
the  polarizability,  a,  is  near  the  reciprocal  of  one  of  the  eigenvalues  of 
the  G  matrix,  then  the  dipole  moment  becomes  excessively  large.  This  is 
suggestive  of  the  type  of  catastrophe  that  occurs  in  the  onset  of  a  ferroelec¬ 
tric  transition.  Such  a  situation  would,  perhaps,  be  modified  by  the  inclu¬ 
sion  of  the  higher  multi  pole  moments  in  the  calculation.  It  should  be  pointed 
out  that  the  G  matrix  defined  In  equations  (14.12)  and  (14.15)  is  dependent 
only  on  the  lattice  constants  and  the  symmetry  of  the  crystal.  The  results 


here  were  expressed  in  terras  of  Cartesian  coordinates  but  can  equally  well  be 
done  in  3pherioal  tensors.  If  higher  moments  were  considered,  the  spherioal 
tensor  form  would  be  much  more  convenient.  (This  statement  has  been  confirmed 
by  M.  Fauoher  in  private  oommuniaation,  1982.  She  has  extended  the  self- 
oonslstent  moments  Uirough  quadrupole  moments.) 

14.2  Self-Consistent  Results  for  Soheellte  Structure 

The  procedure  given  above  is  rather  involved,  so  we  shall  go  into 
the  derivation  of  the  G  tensor  for  the  soheellte  structure  (CaWOip  LiYF^ , 
etc).  The  spaoe  group  for  soheellte  is  88  in  the  International  Tables;  the 
position  of  all  the  constituents  is  given  in  table  14.1.  To  be  specific, 
LiYF^  has  been  chosen;  the  fluorine  is  in  site  1;  x,  y,  and  z  have  been  chosen 
as  the  reference  points  for  the  dipoles  (u,  v,  w);  and  all  other  dipoles  in 
the  unit  cell  are  related  to  u,  v,  and  w.  No  dipoles  can  exist  at  the  Y  or  Li 
sites  since  the  lowest  fields  at  these  sites  are  quadrupole  (k  -  2) . 

To  evaluate  G  for  the  soheellte  structure,  we  choose  the  ion  at 
site  1  in  table  14.1  as  j  in  equation  (14.12).  The  R  in  equation  (14.5), 
including  the  translational  vectors  U  in  x,  m  in  y,  n  in  z),  is 


^i,1  "  U+xi-x)aex  +  (m+y^yJaSy  +  (n+zi-z)cej 


(14.17) 


where  we  shall,  during  this  discussion,  suppress  the  explicit  dependence 
of  1$  on  l,  m,  and  n.  We  write  equation  (14.12)  as 


where 


and 


^  -  f?  *  !?  ■ 


-  - 1  Jp- 

’  i  R?,i 


M  ¥■-  ■ 

1  i ,  1 


(14.18) 


(14.19) 


(14.20) 


where  sums  over  2,,  m,  and  n  are  implicit.  Then  using  table  14.1,  we  write 
F^  explicitly  as 


u 


v  +  _Ji_  +  v 


U  _  V  u 

X  R3  R3  R3  R3  R3  R3  R3 

R1,1  R2,1  R3, 1  R4,1  R5,1  R6,1  R7,1 


(1 4.21) 


v  +  _H _  +  v 


R8,1  R9,1  R?0 , 1  R?1 , 1  R?2,1  R?3.1  "14,1  "15,1  "16,1 


_1_  +  JL_  + 

R?..  .  R?r-  ,  R? 


1  32 


TABLE  14.1.  SPACE  GROUP  88  (FIRST  SETTING):  COORDINATES  OF  ALL 
IONS  IN  A  UNIT  CELL  OF  YLiF^  AND  DIPOLE  MOMENTS  OF  EACH  ION 
(px»  Pyi  and  p2  of  site  1  are  chosen  as  u,  v,  and  w,  respectively) 


Site 

Ion 

X 

y 

z 

Px 

py 

Pz 

<0 

c r 

or 

1 

F 

X 

y 

z 

u 

V 

w 

1 

2 

F 

y 

-X 

-z 

V 

-u 

-w 

(-1)k(i)f‘ 

3 

F 

-x 

-y 

z 

-u 

“V 

V# 

(•  1)q 

4 

F 

-y 

X 

-z 

-v 

u 

“W 

(-1)k(-i)q 

5 

F 

1/2  ♦  x 

1/2  +  y 

1/2  +  z 

u 

V 

w 

1 

6 

F 

•1/2  ♦  y 

1/2  -  x 

1/2  -  z 

V 

-u 

-w 

(-1)k(i)q 

7 

F 

1/2  -  x 

1/2  -  x 

1/2  -  z 

-u 

-v 

w 

(-1)q 

8 

F 

1/2  -  y 

1/2  -  y 

1/2  +  z 

“V 

u 

-w 

(-1)k(~i)q 

9 

F 

X 

1/2  +  y 

1/4  -  z 

u 

V 

-w 

(-1)k+q 

10 

F 

y 

1/2  -  x 

1/4  ♦  z 

V 

-u 

w 

(-i)q 

11 

F 

“X 

1/2  -  y 

1/4  -  z 

-u 

-v 

-w 

(-1)k 

12 

F 

-y 

1/2  +  x 

1/4  ♦  z 

-v 

u 

w 

(i)q 

13 

F 

1/2  ♦  x 

y 

3/4  -  z 

u 

V 

-w 

(-1>k+q 

14 

F 

1/2  ♦  y 

-x 

3/4  +  z 

V 

-u 

w 

(-1)q 

15 

F 

1/2  -  x 

-y 

3/4  -  z 

“U 

-v 

-w 

(-1)k 

16 

F 

1/2  -  y 

X 

3/4  +  z 

-V 

u 

w 

(l)q 

17 

Li 

0 

0 

0 

— 

— 

— 

— 

18 

Li 

0 

1/2 

1/4 

— 

— 

— 

— 

19 

Li 

1/2 

1/2 

1/2 

— 

— 

— 

— 

20 

Li 

1/2 

0 

3/4 

— 

— 

— 

— 

21 

Y 

0 

0 

1/2 

— 

— 

— 

— 

22 

Y 

1/2 

0 

1/4 

— 

— 

-- 

— 

23 

Y 

1/2 

1/2 

0 

-- 

— 

— 

-- 

2k 

Y 

0 

0 

3/4 

— 

— 

— 

— 

aThe  last  column  relates  those  Q^q  to  the  reference  point  Q^g- 
Thus  the  Q^g  for  fluorine  are  all  related  to  site  1. 
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Then  we  can  write 


8 

Fd  -  I 

v  u 


(-1) 


i 


X  1-1  R? 


8  f-i)1 

u  +  l  v  +  <0>w  • 

i-1  R; 


'21+1,1  *  "21,1 

If  we  let  G  -  G’  +  3G"  and  restore  the  i,  m,  n  sum,  we  have 

8  r-u1 

C-X  -  £  I  4^-  • 


t.m.n  1-1  R2i+i,i 


(14.22) 


(14.23) 


G’ 

xy 


8  /  .  \  1 

i  l  411- 

&»tn*n  i*1  Rgj  ‘i 


(14.24) 


Gxz  "  0  • 


(14.25) 


By  similar  methods  we  obtain 


8 


G’ 

zz 


l  l  M  )J 

i,m,n  1-1 


l3 
1+8J 


(14.26) 


and  the  G'  tensor  Is  symmetrical. 


To  evaluate  the  G"  term,  the  procedure  Is  precisely  the  same  as  to 
evaluate  the  G'  term,  except  that  we  relate  G"  to  equation  (14.20).  It  is 
convenient  to  express  ft.  •?  explicitly  In  tabular  form,  as  given  in  table 
14.2,  for  easy  reference  when  writing  out  each  term  of  G" .  We  shall  not  write 
out  the  detailed  expression  as  in  equation  (14.21),  but  this  procedure  gives 


^  (-’>‘[4.1  •  X21Y2il  •  «8.27) 

8  1  ? 

°;y  ■  -  <-')  [x21  -  X21I2l.,]  ,  (U.28) 

8  , 

°xz  •  -  <-’>  [Vl  -  Wl-sl  •  (,"-29> 

8  , 

°?z  "  -  J,  (-,)  [I1Z1  '  Yi.8z1.8l  -  n',-30) 

8 

Gzz  -  "  ^  (_1)  tzi  -  zi+8^  •  <14.31) 
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where 


Vi 


a  x(l )y(l ) 
R5 

i,1 

2  ,, \2 
a  x(l) 

5  ’ 

Vi 

a2y(i)‘i' 

R5 

Ri,1 


XiZi 


acx(i ) z ( i ) 


i*1 


and  all  the  sums  in  equations 
(14.27)  through  (14.31)  have 
the  sum  over  l,  m,  and  n 
implied.  The  G"  is  symmetric 
(this  can  be  shown  directly 
from  evaluating, _  for  example, 
G"  and  G"  independently). 


TABLE  14.2.  VALUES  OF  R*P  FOR  DIFFERENT 
SITES  IN  SCHEELITE3 


Site 

Px 

py 

Pz 

s.? 

1 

u 

V 

w 

x(1 )u  + 

y(1 )v  ♦ 

Z(1  )w 

2 

V 

-u 

-w 

x(2)v  - 

y(2)u  - 

z(2)w 

3 

-u 

-V 

w 

-x(3)u  - 

y(3)v  + 

z(3)w 

Ji 

-v 

u 

"W 

-x(4)v  ♦ 

y(4)u  - 

z(4)w 

5 

u 

V 

w 

x(5)u  + 

y(5)v  + 

z(5)w 

6 

V 

-u 

-w 

x(6)v  - 

y(6)u  - 

z(6)w 

7 

u 

-v 

w 

-x(7)u  - 

y(7)v  + 

z(7)w 

8 

-v 

u 

-w 

-x(8)v  + 

y(8)u  - 

z(8)w 

9 

u 

V 

-w 

x(9)y  + 

y(9)v  - 

z(9)w 

10 

V 

-u 

w 

x(10)v  ■ 

-  y(10)u 

*  z(10)w 

11 

-u 

-v 

-w 

-x(ll)u  ■ 

-  y(1l )v 

-  z( 1 1 )w 

12 

-V 

u 

w 
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r 

The  equations  for  G’  and  G"  were  calculated  for  several  lattices, 
and  the  results  are  given  In  table  14.3.  The  crystal  axial  field  components, 
A_0,  were  computed  for  CaWOjj  using  o  -  2.4  A  and  oxygen  charge  of  -2e,  and 
using  a  -  0.24  A'  and  oxygen  charge  of  -e.  The  results  are  shown  in  table 
14.4  (Morrison,  1976). 


After  the  above  work  hai‘  been  done,  the  dipole  terms  in  the 
were  programmed  for  a  computer  for  a.ll  the  230  space  groups.  In  the 


TABLE 

14.3.  0 

TENSOR  AND  X- 

RAY  DATA  FOR  SEVERAL 

COMPOUNDS  (1/A* 

}a,b 

Compound 

a 

c 

X 

y 

z 

CaWOjj 

5.248 

11.376 

0.2413 

0.1511 

0.0861 

PbMoOj, 

5.4312 

12.0165 

0.2353 

0.1  3660 

0.08110 

tLlF|i 

5.1668 

10.733 

0.2820 

0.1642 

0.0815 

yvo„ 

7.120 

6.289 

0.1852 

0 

0.1749 

Gxx 

Gxy 

Gxz  Gyy 

V 

Gzz 

CaWOj, 

-0.252608 

-0.0731076 

-0.0979719  0.197101 

-0.040229 

0.123419 

PbMoO^ 

-0.224152 

-0.0758124 

-0.0976969  0.168517 

-0.0426228 

0.0768659 

YL1F,, 

-0.252131 

-0.0826652 

-0.125165  0.192571 

-0.0446824 

0.168842 

yvo„ 

-0.734174 

0 

-0.241965  0.127677 

C 

-0.173784 

^^The^eferenceaitelTinallthecalculatlona ,  ia  the  ligand  at  x,  y,  and  z. 

ia  not  a  s?heelite  structure  ( YVO *  ia  the  zircon  atructure,  apace  group 
141,  in  the  International  Tablea)  but  can  be  done  in  the  acheelite  atructure  by 


t ranalating  the  oxygen  poaitiona  to  the  above  (aea  Karayiania  and  Moxriaon,  1973). 
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TABLE  14. 4.  AXIAL  COMPONENTS  OF  CRYSTAL  FIELD  FOR  TWO 
VALUES  OF  POLARIZABILITY  OF  OXYGEN  AND  OXYGEN  CHARGE 

(om“  /A* ) 


program  any  number  of  inequivalent  3ites  can  have  an  associated  dipole  moment 
(we  only  considered  one  type  of  site  above).  Recently  the  members  of  Caro's 
group  in  France  and  de  Sa's  group  in  Brazil  (Fauoher  and  Malta,  1981)  have 
included  the  dipole  and  quadrupole  moments  in  a  self-consistent  manner  for 
LaCl^;  they  have  found  that  with  the  reported  values  of  the  dipole  and  quadru¬ 
pole'  polarizabilities  the  resultant  A9  is  much  larger  than  A0  or  Ad  . 

nm  nm  nm 

1 4.3  Self-Induced  Effects 

When  a  rare-earth  ion  is  immersed  in  a  solid  it  is  possible  for  its 
electrons  to  experience  a  field  due  to  the  reaction  of  the  medium  baok  on  the 
electrons.  Both  this  type  of  field  and  the  external  fields  due  to  the  point 
charges  of  the  medium  can  exist.  Thi3  reaction  is  identical  to  the  classical 
problem  of  a  charged  particle  interacting  with  its  induced  image  in  a  conduct¬ 
ing  plate  or  sphere.  The  interaction  was  recognized  by  Judd  (1977),  and  it 
was  he  who  suggested  the  polarization  of  the  ligands  as  a  possible  source  of  a 
two-electron  crystal-field  interaction.  In  this  section  we  consider  the 
development  of  this  interaction  as  derived  earlier  (Morrison,  1980),  using  the 
same  technique  used  in  the  earlier  work.  In  later  sections  this  interaction 
is  developed  in  a  more  general  way,  deriving  the  multipolar  interaction. 

Wo  consider  an  electron  at  r  on  a  rare-earth  ion  and  a  ligand  at  ft 
with  polarizability  a.  The  electric  potential  created  by  the  electron  is 


*  -  —  .  (14.32) 

I*  -?l 

The  electric  field  at  the  ligand  is 

ft  »  -VR4>  , 


where  VR 
Then, 


indicates  that  the  derivative  should  be  taken  with  respect  to  ft. 


ft 


-e(ft  -  r? 

I*  -  ?i3 


(14.33) 


i 

> 
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The  dipole  moment  induced  on  the  ligand  is  given  by 


aft  , 


(14. 3*0 


where  a  is  the  polarizability  of  the  ligand. 

Now  if  we  consider  a  dipole  from  some  arbitrary  origin,  the  eleo- 
trlc  potential  at  point  from  that  origin  is 

p-S, 


(14.35) 


To  find  this  potential  at  the  electron  it3eif,  we  let  ^  -  -(ft  -  r).  Then 
equation  (14.35)  becomes 


-p  •  (ft  -  r) 


|5-?|3 


The  energy  of  the  electron  interacting  with  this  potential  is  given  by 


*  *1  -G  i  /_ \ 

T  *i(r) 


U(r ,ft) 


e  p  «  (ft  -  r) 

2  i*  -  ?|3 


where  the  1/2  is  due  to  a  self-interaction.  We  can  write 


ft  -  r  .  _v  _ 1_ 

|ft  -  r|3  "  R  |ft  -  r  ’ 


Then  equation  (14.37)  becomes 


(14.37) 


(14.38) 


U<M>  -=£?  ■  VR— V  , 

I*  -  fl 


(14.39) 


and  similarly 


2  =  eV  - - - 

R  |ft  -  ?| 


(14. HO) 
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Using  the  result  of  equation  ( 1  *1. 40)  in  equation  (I1!. 34)  and  substituting  the 
result  Into  equation  (14.39);  we  have 


(14.41) 


whers  Vo  operates  only  on  the  function  on  its  immediate  right.  To  further 
reduce  the  result  given  in  equation  (14.41),  we  consider  the  operation 


+  ^V2^  +  2(7^)  •  (v<i>2) 


If  tjii  and  \l>2  satisfy  Laplace's  equation  (whioh  they  do),  then 


^2)  “  2(^i  )  *  (^2)  • 


(14.42) 


(14.43) 


If  we  identify  ^  and  i|>2  with  1  / 1 S  -  r|  in  equations  (14.41)  and  (14.43),  we 
can  write  U(r,R)  as 


u(M) 


“ae  „2  1 

T  7R  ,'i  *,2  ■ 


is  -  ?r 


(14.44) 


To  proceed  further  we  must  expand  the  factors  on  the  right  side  of 
equation  (14.44).  First  we  notice  that 


If  we  let 


I* 


-  r|2  -  R2  +  r2  -  2r-R 


-  2Rr 


|~r2  +  r2  ^  1 

L  2Rr  "  r‘RJ 


_2  2 
R  t-  r 

2rR 


(14.45) 


(14.46) 


then 


The  expansion 


-  l  (2n  ♦  1)Q(t)P(z) 
t  -  z  **  n  n 

n 

is  given  by  Rainville  (I960);  the  leading  terra  for  larga  t  is 

an(t; 

From  equation  (14.46)  we  have 


a"(ni)2 
tn*’(2n  *1)1 


«2n+1 .2  n+1 
Qn(t)  “(2n  ♦  ijl  Rn+1 

for  large  R. 

Substituting  the  result  of  equation  (14.50)  into 


gives 


1 


02n,  , ,2  n 

y  2  r  P  (Z) 

1  ,n+2  nU'  * 


|R  -  r|‘  (2n)!R‘ 

From  the  Legendre  addition  theorem  (see  seot.  1),  we  have 


P  (z)  -  P  (r*R) 
n  nv  ' 


ana 


R  -  r  |: 


l  C  (?)C»  (R) 
u  nnr  ‘  nnr  J 
nfm 


n,m 


The  remaining  necessary  operation  is  VR,  which  can 


-  4  a-  (f.2  5-) 

R  R2  dR  \  dR/ 


t2 

R2  ‘ 


The  only  term  in  equation  (14.53)  that  this  operates  on  is 


~  C*  (R) 

4-  v-  Mill  v  * 


pn+«i  nm 


(14.48) 

(14.49) 

(14.50) 

equation  (14.47) 

(14.51) 

(14.52) 

(14.53) 

be  written 

(1 4.54) 

(14.55) 
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Using  equations  (14.54)  and  (14.55),  we  have 


Finally, 


,g[R-"-^(R)]  -  .  at-^cyR) 


where  we  have  uoed 


(tPcjR)  -n(n*  1)cjR)  . 


(14.56) 


(14.57) 


(14.58) 


a  result  we  discussed  in  section  1.  The  result  in  equation  (14.57)  substi 
tuted  into  equation  (14.53)  gives 


VI 


|R  -  r |: 


22n+1 (ni ) (n 

l  72nT T 

nm 


+  D! 


Rn 


crJ?>c;u(5) 


(14.59) 


The  result  given  in  equation  (14.59)  is  substituted  into  equation  (14.44). 
This  result,  when  summed  over  all  ligands  at  R  with  polarizability  ,  pro¬ 
duces  J 


U  (r ) 


nt (n  ♦  1 )  t 
(2n) ! 


C*  (r.) 


rnC  (r 
run'* 


) 


(14.60) 


If  we  write  equation  (14.60)  as  we  have  previously  done  with  the  point- charge 
model , 

U(r)  -  l  (ASI)*  r"C  f?)  ,  (14.61) 

nm  nm 

nm 

we  have 


*£  •  (4)  1  •  »»•«> 

^  J 

which  are  the  seif-induced  crystal-'ield  components  due  to  induced  dipoles 
only.  Higher  order  multipole  moments  can  be  induced  on  the  ligands,  and  these 
multipoles  will  contribute  a  correction.  From  previous  experience,  we  should 
anticipate  the  total  self-induced  multipole  fields  to  be  of  the  form 


k-1 


ASI(X) 

nm 


SI 

with  the  result  above  being  A  (1). 

nm 


(14.63) 


1  40 


•J«iAIWyUMMMJUUMUMJw;yuuuWiAI'dUMil 


As  in  the  point-oharge  model,  if  we  express  all  lengths  in 

angstroms  and  otj  in  angstroms  oubed,  then  equation  (1*4.62)  becomes 

•  -(&  *  ’°8)  -~-ri?7r—  I  •  <«•«> 

To  express  in  units  of  om“1/An,  use  the  conversion  faotor  a_/8Tr  x  10®  - 

29,035.  ""  0 
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15.  MISCELLANEOUS  CRYSTAL-FIELD  EFFECTS 


15.1  Judd's  Interaction  for  Two  Electrons 


The  interaction  considered  here  is  a  development  of  a  suggestion  by 
Judd  (1977)  concerning  a  possible  origin  of  two-electron  crystal-field  ef¬ 
fects.  Specifically,  Judd  suggested  that  such  term3  would  arise  if  one  of  the 
electrons  in  the  configuration  nS,N  polarized  a  nearby  ion,  and  the  remaining 
N  -  1  electrons  interacted  with  the  induced  multipolar  moments.  The  investi¬ 
gation  of  this  interaction  was  performed  later  (Morrison,  1980),  assuming  only 
a  dipole  polarizability.  The  interaction  for  two  electrons  that  resulted  is 


V(  1 , 2, R )  ■  l  F(abk)raCaa(P1  )rbcbjq_a(P2)  <a(o)b(q-a)|k(q)>  (15.1) 

b.k.q 

where 


F(abk)  -  -  (a  <a(0)b(0)  | k(0) >  [ (a+b+1 ) (a+b+2)  -  k(k+1)]  , 

and  a  is  the  dipole  polarizability  of  the  ion  at  R. 

The  development  of  the  result  given  in  equation  (15,1)  was  similar 
to  that  given  in  the  derivation  of  the  self-induced  field  in  section  14.3. 
For  the  full  multipolar  result  we  shall  use  more  general  methods. 

The  electric  potential  of  an  electron  at  r. ,  as  seen  at  a  ligand 

at  is 


(15.2) 


where 


The  multipolar  inducing  field  Elnm  at  ^  oan  be  defined  by 


*(S’  *  *!  ‘  'J  • 

nm 


(15.3) 


By  expanding  equation  (15.2),  we  obtain 


C»  (Ri)xn 

c»<:>  • 

nm 


(15. n) 


then  oomparing  equation  (15. *0  with  equation  (15.3)  gives 


E„m  “  e(-DnC  (R.,)/Rn+1  . 

nm  run  I  l 


(15.5) 


The  multipole  moment,  Q^,  is  given  by 


*nm 


°nE 


nm 


(15.6) 


where  the  multipole  polarizability  is  o^. 

The  electric  potential  at  an  arbitrary  point  from  a  multipole 
distribution  is  given  by 


,  rfi  )  y  QnmCnm^3^ 

W  1  „n.1 


(15.7) 


and  the  energy  of  an  electron  at  r2  interacting  with  the  multipoles  is 


U  -  -ed»Q  ($3  -  -R2) 


(15.8) 


with  R2  “  R  ”  r2' 


From  equation  (15.7),  we  obtain 


•e  I 
nm 


Q*  (-1  )nC  (R_) 
nm  nm  2 

Rj+1 


(15.9) 


Now  from  equations  (15.6)  and  (15.5)  we  have 

,nCnm(S1> 


Q  =  ea(-1 ) 
nm  n  pn^l 


(15.10) 


HU 


whioh,  when  substituted  into  equation  (15.9),  gives 


If  we  write 


■W>  •  -  „K 

nm 


u(f,?2.5)  -  1  u<n)(?,?2,S)  . 


nm 


(15.11) 


(15.12) 


we  have 


•  u(")(?,?2.*)  -  -.2a„  l 


Cm,f52) 


"  .  «?*'  RT' 


(15.13) 


f  we  were  considering  the  self-interaction,  at  this  point  we  would  let  * 
fti  and  take  half  the  results.  The  sum  on  m  would  then  collapse  to  unity  and 
-  -e2(an/2)/R2n+2. 

However,  the  two-eleotron  interaction  is  more  complicated.  We  use 
the  two-center  expansions  (Carlson  and  Rushbrooke,  1950;  Judd,  1975)  to  obtain 

i  <a(a)n(m) |a+n(a+m)>  • 

(15.1*0 


where  ^  ?  -  r1 ,  and 


CniJ*y  r  /2b+2 n\  1 /2  roc£fl(rJ 


R£+1 


j  CT)  <»<6)n(.)|».n(B^)>  , 


(15.15) 


where  -  $  -  r^, 


As  Indicated  in  equation  (15.11),  equations  (15.1*0  and  (15.15)  are 
to  be  multiplied  together.  When  these  two  equations  are  multiplied,  the  two 
spherical  tensors  in  R  can  be  recoupled  as 


CI.n,a„»(S>Cb.nWS)  '  (-,)C“ra  J  «*n(0)..n(0)|k(0)> 

K 


(15,16) 


<b+n(8+m)a+n(-ot-m)  |k(6-o)>  C  (fi)  , 

k  f  m  "a 
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X 


where  we  have  used 


c;.n,a.rJS’ 


(-1  )o»+mc 


a+n ,-a-m 


(R) 


(15.17) 


It  should  be  noted  that  the  resultant  projection  J.n  equation  (15.16), 
[Ck,0-a(R)  ] .  is  independent  of  m.  Thus  with  a  proper  recoupling  of  the  C-G 
coefficients  in  equations  (15. 11*)  and  (15.15),  the  sum  over  m  can  be  per¬ 
formed.  Selecting  the  independent  terms  from  the  product  of  equations  (15. 1*0 
and  (15.15)  and  the  result  of  equation  (15.1 6),  we  have 


S  -  £  (-1)m+a  <a(a)n(m)  |a+n(cfm)>  <n(m)b(B)  |b+n(B+m)> 
m 


(15.18) 


x  <b+n(6+m)a+n(-a-m) |k( B-a)> 


which,  when  further  reduced,  gives 


U(n)(r1?2,^)  -  -e2an  I  l  <b*n(0)a+n(0)  | k ( 0) >  [(2a^n)  (2b^n)] 


aa  k 

be 


(15.19) 


.....  .  A  . 

x  r1Caa^r1 ^r2CbB^r2^  Ra+b+2n+2  * 


Thus  the  final  desired  result  is  obtained  if  we  knew  S.  In  equation  (15.18) 
we  rearrange  the  C-G  coefficients  as  follows: 


<a(ot)n(m)  ja+n(a+m)> 
<n(m)b( B) j  b+n(B+m)> 


(-1)3  a  <a(-a)a+n(a+m)  |n(m)>  , 

(15.20) 

1  /2 

(-1)n  m  - gg-  ;-n1  ^  <n(m)b+n(-B-m) |b(-6)>  . 


We  then  recouple  (eq  (5.8))  the  two  C-G  coefficients  on  the  right  to  give 


<a(-a)a+n(o-*m)  |n(m)>  <n(m)b+n(-B-m)  |  b(-B)> 


-  £  /(2f +1 ) ( 2n+1 )  W(a,a+n,b,b+n;nf )  <a+n(o+m)b+n(-B_m) |f (a-B) > 
f 

X  <a(-a)f (a-B)|b(-B)>  . 


(15.21) 


1  H6 


The  sum  on  m  oan  now  be  performed;  note  that  the  phase,  (-1)m,  in  equation 
(15.20)  oanoels  the  (-1)ra  in  equation  (15.19)  if  we  ohange  the  phase  in  the 
first  C-G  ooeffioient  on  the  right  side  of  equation  (15.21).  This  then  fixes 
the  sums  on  f  at  k.  Thus, 

„  ,  <sk-b+nr(2a+2n+1)(2b+2n+1) 
s  -  (-1)  |_ - - 


(15.22) 


*  W(a,a+n,b,bn;nk)  <a(-a)k(a-S) ( b ( —0 ) >  . 


The  C-G  in  equation  (15.22)  can  be  rearranged  to  give 

S  -  (-1)°[(2a*2n+1)(2b+2n+1 ) j1 /2W(a,a+n,b,b+n;nk)  <a(a)b(p) j k(a-0)>  . 

(15.23) 


If  we  let 


F  (abk)  -  -fa  e2)  <a+n(0)b+n(0) |k(0) >  /(2a+2n H )(2b+2n+ I ) 
n  n  1 

(15.24) 

*  W(a,a+n,b,b+n;nk)^aJ)(2b2^)]  , 

then,  substituting  into  equation  (15,19),  we  have 

U(n)(^r2.R)  -  l  Fn(abk)rar£  I  <a(a)b(q~a)|k(q)>  )cb  (r2) 

a,b  a 

k,q 

(15.25) 


Rq+b+2n+2  ’ 

which  is  the  final  form  of  the  two-electron  multipolar  Interaction.  To  obtain 
the  result  given  in  (15.1),  we  would  have  to  relate  <a+1 (0)b+1 (0) |k(0)>  to 
<a(0)b(0) ] k ( 0) >  and  evaluate  W(a,a+1 ,b,b+1 ;1k),  both  of  which  procedures  can 
be  found  in  Rose  (1957,  pp  47,  277).  If  this  is  done,  then  equation  (15.25) 
will  reduoe  to  equation  (15.1).  In  a  solid  the  ligands  at  R  are  suoh  that, 
when  the  sum  is  performed  over  the  ligands,  only  certain  k  and  q  survive. 
Much  of  the  above  derivation  has  been  given  by  Judd  (1976)  In  a  different 
context,  and  many  of  his  elegant  techniques  could  be  used  to  simplify  the 
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resulting  expressions.  For  example,  using  Judd’s  notation  (1976),  equation 
(15.25)  beoomes 


-  I  Fn(abk)r^r|[Ca(ri  )c&(r2)]k  •  Ck(R)/Ra+b+2n+2 


ab 

k 


(15.26) 


where 


[S.(?1)SbMk(1  -  I  <a(.)b(q-.)|k«i)>  Cj?,)^  q.0(^)  . 

The  tensor  in  orbital  space,  given  in  equation  (15.25), 

Tk  (a,b)  -  l  <a(a)b(q-a)  jk(q)>  C  (r2)  .  (15.27) 

q  a  *q 


should  be  considered  carefully.  For  a  fixed  value  of  k  the  number  of  terms  In 
the  sum  over  a  and  b  is  restricted  by  a  +  b  <  k;  for  equivalent  electrons  a 
and  b  are  restricted  to  even  Integers;  and  for  0  <  (a,b)  <  6,  the  total  number 
of  terms  is  not  excessive.  But  since  a  and  b  can  reach  the  maximum  value  of  6 
for  the  configuration  4fn,  the  value  of  k  in  the  k  sum  (similar  to  the  lattics 
sum)  must  go  up  to  12,  that  is,  k  <  12. 

If  as  in  previous  work  (Morrison,  1980)  the  sum  over  all  the  elec¬ 
trons  is  performed  ir  equation  (15.25)  along  with  the  sum  over  the  ligands, 
the  results  are 


u(n)(8)  . 


1/2  l  U 
ij 
R 


(n)fi  t 


,$) 


(15.28) 


where  the  factor  1/2  accounts  for  the  self- intar act  ion  terms  that  are  present 
when  an  electron  interacts  with  its  own  induced  multipole,  as  well  as  for  the 
interactions  that  occur  twice  when  i  *  J.  This  interaction  contains  a  large 
number  of  corrections  to  the  free-ion  parameters,  a  few  o^  which  are  discussed 
in  the  following. 

15.2  Slater  Integral  Shifts 


The  Slater  integrals  for  the  free-ion  interactions  are  given  by  the 
Coulomb  interaction  as 


H 


1 


l  — 5 - 

1>J  |ri  -  r j | 


(15.29) 
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which  for  equivalent  electrons  can  be  written 


where 


(15.30) 


.(k) 


-  e 


I  1 


r< 


0  0  r> 


k  +  1 


[Rnji^r2^2  dridr 


1  2 


Since  the  interaction  represented  by  equation  (15.30)  is  spheri¬ 
cally  symmetric  in  the  space  of  all  the  electrons,  corrections  to  the  F'*)  can 
only  arise  from  terms,  in  an  interaction  that  are  spherically  symmetric  in  the 
space  of  the  electrons.  Thus,  in  equation  (15.25)  if  we  let  k  -  0,  we  have 
such  an  interaction,  and  the  following  results  are  achieved: 


F  (abO)  -  -a  e2  <a*n(0)a+n( 0) I 0( 0) >  /(2a+2n+1)  ( 2b+2n+1 ) 
n  n  1 

*  W(a,a*n,a,a+n;n0)  [(^a")  (15.31) 

If  in  equation  (15.31)  we  use  the  relations 

<a(a)b(-a)|0(0)>  -  (-1  )a"a//2bTT  6^  (15.32) 


and 


W(a,a+n,a,a>n;n0)  -  (-1 )n/[ ( 2a+1 ) (2a+2n+1 ) j1 


(Rose,  1957),  then  equation  (15.25)  becomes 

U<n>(iy2.Sj  -  -V2  l  (2>*2n)  I  rfr-a 

which  is  the  same  form  as  equation  (15.30). 


C*  (r.  )C  (r  J 
aa  1  aa  2 

p2a+2n+2 
Thus , 


(15.33) 


- 1 


Vi,e 


.  kv2 
<r  > 

,2k+2n+2 


(15.34) 


for  the  Slater  integral  shifts,  because  of  the  electron  multipolar  interaction 
with  the  ligands  of  multipolar  polarizabilities  a^. 


15.3  Problems 


i.  If  we  ijave  a  oharge  distribution  p(x)  at  the  origin,  then  the 
eleotrio  potential  at  r  is 

_>  p(x)dt„ 
d$(r)  ■ 


|r  -  x 


where  dix  is  the  volume  element  at  x.  By  expanding  - - -  for  Ixl  <  Irl,  show 

that  |r  -  x| 


♦  (r) 


1 

nra 


Q#  C  (r) 
nm  run 

n+1 


where 

Qnm  “  /  P(>:)xnCnnj(x)  dtx 

Is  the  multipolar  moment  of  the  oharge  distribution.  What  is  the  common  name 
for  Q00  and  QlB? 

+2.  If  we  consider  the  charge  distribution  in  problem  1  to  be  at 
the  point  r,  show  that  the  electrio  potential  is 

Q*  C  (r) 

V  /  . \ n  nm  nm' 

♦(r)  •  i  (  °  —  • 

nm  1 

3.  The  electric  potential  of  a  uniform  electric  field  ft  is 


$(x)  -  -x*f 
or 

♦<*>  •  - z  <*> 

m 

when  £  and  x  are  written  in  spherical  tensor  form. 

If  we  generalize  equation  (a)  above  to 


♦(ft  ♦  x)  -  -  J  B*  (ft)x"C  (x)  , 

__  nm  nm 

nm 


(b) 


we  can  define  the  induced  multipoie  moment  caused  by  a  field  E  (ft)  as 

nm 

Q__ (ft )  -  a  E  (ft)  . 
nm  n  nm 


(o) 
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3y  taking  a  point  charge  q0  at  R,  show  that  the  multipolar  field  at  the  origin 
is 


CJR) 

_  nm 

run  "  “qo  Rn+1 


and  that 


Qnm(R) 

nm 


-q  a 
Mo  n 


cn_( R ) 

nm 

Rn+1 


(d) 


Using  the 
potential  at  r  due  to 
polarizability  an  is 


results  of  problems  1,  2,  and 
the  above  multi  pole  moments 


*  ■  l  °nVR> 

nm 


C  (r) 
nm 

rn+1 


3,  show  that  the 
of  an  ion  with 


with  Q^CR)  given  in  the  last  part  of  problem  3* 

5.  By  considering  a  dielectric  sphere  of  radius  a  and  dieleotric 
constant  e  in  the  presence  of  a  point  charge  qQ  at  |r|  >  a,  show  that  the 
multipole  polarizability  of  the  sphere  is 


c<n  -  [n(e  -  1)a2n+1]/(£n  *  n  +  1) 


This  result  can  be  obtained  by  solving  the  electrostatic  problem  of  a  dieleo 
trio  sphere  in  the  presence  of  a  point  charge  (Jackson,  1975). 

^  6.  The^  energy  of  interaction  of  two  rig^d  charge  distribu' 
tions  pA(x)  and  PB(y)  separated  by  a  distance  ft  (|ft|  >>  |x|,  |y|)  is 

pk(h  dixPB(y)  <lty 

^ab  "  *  J  ,  ±  *  ■*  i 

jft  -  X  +  yj 

By  using  the  two-center  expansion  (Judd,  1975)  on  the  integrand,  show  that 

UAB  ‘  £  <a(a)b(6)|a+b(a+B)>  Qaa(*)a„s<B) 

ad 
bB 

where 

Qaa(A)  "  /  pA(^x\a(^  dTx  •  etc‘ 


15! 


The  total  energy  for  a  solid  consisting  of  suoh  multipoles  would  then  be  tne 
sum  over  all  A  and  B  suoh  that  Ft  -  R.  -  Ra  . 


7.  The  energy  of  Interaction  of  two  ions,  with  the  first  having  the 
electronic  configuration  r>A®*AA  an<*  the  aeoond  ion  having  the  electronic 
configuration  ngig8,  can  be  obtained  from  the  result  .1.r.  problem  6  by  letting 


QaoU)  ‘  J,  1Caa(?l) 


and 


N, 


abB(B) 


The  resulting  interaction  can  be  used  to  calculate  energy  transfer  from  ion  A 
to  B. 
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